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£Sj | This thesis is a theoretical analysis of sample two- and one- dimensional systems. The two- 

■ dimensional examples are the quantum Hall liquid and anomalous paired states. The most 
widely accepted effective theory of the quantum Hall liquid is based on the so-called Chern 

■ Simons Lagrangian, but it is not entirely satisfactory. We obtain the first derivation of an 



alternative effective theory from microscopic principles. Our formulation allows for a first 

principles derivation of physical quantities such as the effective mass and compressibility 

and contains the first analytical observation of the magnetoroton. The formalism developed 

along the way is also applied to paired states in anomalous supeconductors, a topic of 

much recent interest. The one-dimenisonal system is the carbon nanotube. Gas uptake 

in nanotube bundles is currently attracting a wealth of research with both applied and 

fundamental implications. We propose adsorption of gases on the surface of a single tube, 

finding strong correlations and symmetries that have not been observed yet. The properties 

■ of these states are directly relevant to other one-dimensional structures such as spin ladders 

IT) , . 

t^j- and stripes and raise interesting and open questions. 
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Chapter 1 



Introduction 



1.1 Composite Particles in the Fractional Quantum Hall Effect 

The technological innovation that made discovery of the quantum Hall effect (QHE) possible 
is called a MOSFET, or metal oxide semiconductor field effect transistor. Under suitable 
conditions, electrons are effectively confined to two dimensions by an inversion layer Q. 
Inversion layers are formed at an interface of a semiconductor and an insulator or between 
two semiconductors with one of them acting as an insulator. The original system in which 
the QHE was discovered was between Si (semiconductor) and SiC>2 (insulator). More re- 
cently the semiconductor-semiconductor system GaAs-ALjGai-^As (with GaAs acting as 
the semiconductor) has been used. The parameter is approximately x ~ 0.2. Typically the 
layers of GaAs and AlGaAs are grown with atomic precision using molecular beam epitaxy. 
The necessary donors that are required for the inversion layer to form are implanted away 
from the interface allowing very high electron mobility within the inversion layer. For some 
samples, an external voltage may be used to control the density of electrons. 

The original discovery by von Klitzing, Dorda, and Pepper in 1980 has come to be 
known as the integer quantum Hall effect (IQHE) . In brief, they found that when a magnetic 
field is applied perpendicular to the electron layer, the electron current response is purely 
transverse and quantized. More precisely, the current density responds to an electric field 
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by ji = Ylj a ijEj, where the conductivity tensor is 







ve 2 /h 



a = 



ve 2 /h 







h is Planck's constant, e is the electron charge, and v is a small integer. Thus the con- 
ductivity is quantized in fundamental units, independently of specific material parameters. 
The diagonal conductivity vanishes, so the state is dissipationless. The actual measured 
quantities are the conductance and/or resistance. Theoretical understanding of this effect is 
based on models of independent electrons in the presence of disorder, and is quite developed 
by now. 

In this thesis we will be concerned with an effect that was discovered shortly after 
the IQHE, which is known as the fractional quantum Hall effect (FQHE). In 1982, Tsui, 
Stormer, and Gossard found that, in extremely pure samples at very low temperature, 
the integer v above can be replaced by a hierarchy of rational numbers, v = p/q. In the 
quantum limit uj c t 3> 1, where oj c = eB/m and r is the electronic scattering time, v 
shows plateaus at these fractions as the chemical potential is varied. The leading fractions 
were found to have q =odd, forming an incompressible quantum liquid. Since then, q = 
even states are understood to have their own complementary set of phenomena [||. In 
contrast to the odd denominator fillings, they are compressible Fermi liquid-like states and 
are not characterized by plateaus in the conductivity (at least for p < 3; v = 5/2 may be 
an exception Q). Theoretical understanding of the FQHE has progressed rapidly, but is 
not yet complete. For instance, problems which depend in detail on the interaction and 
disorder, such as transitions between the plateaus, are not well understood. Questions have 
also arisen recently on the nature of the effective theory for v < 1 even in pure samples. 
This is the issue we will tackle in this part of the thesis. 

Let us briefly summarize the salient ingredients of the theory. The fundamental length 
scale for electrons in a magnetic field B is 



particle states are parameterized by Landau levels (LL) of energy E n = hco c (n+ 1/2). Each 




It is independent of material parameters and is in the range of 50 — lOOA. The independent 



9 



LL is highly degenerate; the number of states per unit area of one full LL is given by 

po = 1/2tv£% = eB/hc 

The last equality can be rewritten as po = B/<&q, where <E>o is the unit flux quantum, so 
that the degeneracy of one full LL is counted by the number of flux quanta in the external 
field. In this simplified model, electrons successively occupy the Landau levels, and the 
proportion of occupied states is denoted by the filling fraction, v. If the electron density is 
p then 

v = p/p . (1.1) 

The IQHE occurs at integral u, when an integral number of levels are full. 

However, the FQHE occurs when certain rational fractions of LL's are filled. The 
theoretical understanding of this phenomenon starts with Laughlin's approach in 1983 Q. 
He proposed a variational wavefunction for the ground state at v = 1/p, p =odd, which is 
a fractionally filled lowest Landau level (LLL). Assuming that co c is large compared to the 
electron-electron interaction, only the LLL should describe the physics at low energies. For 
N electrons with coordinates z% = X{ + iyi, the Laughlin wavefunction is, 

N N 
i<j i 

This function is composed of single particle states in the LLL and is properly antisymmetric 
in keeping with fermionic statistics. The two main features of are (i) there is a zero on each 
electron, and (ii) each electron sees other electrons as magnetic flux due to the accumulated 
phase in dragging one coordinate around another. The basic excitations are quasiholes 
and quasielectrons with fractional charge ±e/p. Clearly this kind of effect is due to strong 
correlations in the fluid. The quasiparticles also obey fractional statistics, as articulated by 
Halperin || and by Arovas et al. || . Experimental data is consistent with fractional charge 
at v = 1/3 and that the basic excitations at v = 1/2 are neutral M, p|, M. Already we see 
a collective behavior that is very different from Fermi liquids. In fact, the Laughlin state 
describes a strongly correlated quantum liquid that cannot be reached perturbatively from 
the Fermi liquid. As originally noted by Laughlin, it is accepted that the FQHE liquid is 
isotropic and incompressible. 
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Much of the current understanding of quantum liquids relies on effective field theories. 
In the case of the FQHE the most successful has been the Chern-Simons (CS) theory 
Fortuitously, CS actions were being formally developed concurrently in a purely field 
theoretic context [11 1. One starts with a transformation that represents each electron as 
a boson plus an odd number of <ft of <5-flux tubes. Fractional statistics emerges as a Berry 
phase when particles or vortices are dragged around flux tubes. After the transformation, 
the action contains a U(l) Chern-Simons term that couples to the fermion density. At the 
fractions v = 1/q (q , =odd), the statistical gauge potential a is determined by the relation 

C = ~Vxa (1.2) 

(where Ti = c = 1). The full Lagrangian includes a by minimal coupling and the statistical 
gauge transformation allows us to replace the fermions by a bosonic field, <j>. The complete 
gauge invariant Lagrangian is 

Cos = <t>Hidt -Ao- a Q )<f> - -L |(iV + A + a) 0| 2 + ^Le^\d v a x 

2m 47r 



d 2 yp(x)V(x - y)p{y) , (1.3) 



where Aq, A is the external potential. The third term is the U(l) CS term (e is the Levi- 
Civita symbol) and the last term is the interaction. The defining equation ( |1.2| ) for a can 
be viewed as an equation of motion. To completely determine a, one can use conservation 
of charge, dtp + diJi = 0, to obtain the dynamics equation, Jj = —-^e^dtaj. Together with 
eqn. ( |1.2|) this completely determines the statistical gauge field a. 

The net field is the sum of the external field and the 5-flux tubes. On the average, for 
a uniform density of particles, the two fields cancel when v = l/<j) (that is, (A + a) = 0) 
and we are left with bosons in zero net field. The Laughlin state can thus be interpreted 
as a bose condensate of the electron-flux tube composite. In principle, one can substitute 
bosonic particles for electrons and use an even <fi. Some understanding of the states with 
even denominator, e.g. v = 1/2, has been achieved in this framework In this case, the 
composite fermions form a Fermi sea and the resultant phase is Fermi-liquid like. Paired 
states of fermions or bosons are envisioned as pairing of particle-flux tube composites in 
zero field, since the external field vanishes at mean field, as discussed in the introductory 
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Section [E^ and in Chapters || and |5[ Any analysis beyond mean field must proceed with 
caution in any case; standard techniques are exact only in the limit <j> — * 0, but the FQHE 
states require (ft > 2. 

Other, perhaps more fundamental, difficulties with the CS approach have been appreci- 
ated by several workers since the beginning (see e.g. [jT^] ). At zero temperature, the usual 
assumption is that the inter-electron interactions ~ v l l 2 e 2 j 's£b are weak compared to the 
cyclotron frequency oj c so the physics should be dominated by LLL states when v < 1. The 
kinetic energy is just a constant in any given LL so we are faced with a macroscopically 
degenerate perturbation theory with a purely dynamical Hamiltonian. In particular, con- 
sider the Fermi liquid-like state at v = 1/2. In the mean field approximation, the effective 
mass of excitations close to the Fermi surface is the bare mass m, however the low-energy 
excitations should have an effective mass m* determined solely by the interactions. The 
problem is partially resolved by Fermi liquid fluctuations, which introduce a Landau inter- 
action parameter F\ that renormalizes the bare mass by to -1 = m* -1 (l + F\). The most 
serious problem is that interactions do not play a role in the compressibility; whether or 
not interactions are included in the fluctuations, the result is a finite compressibility for the 
Fermi liquid due to the mass m. On the other hand, a partially filled Landau level of non- 
interacting particles should have an infinite compressibility. The same difficulty is present 
in the odd denominator filling fractions; regardless of the interaction, the compressibility 
vanishes flOfl . There has been a renewed interest in these puzzles recently, stemming from 
new developments in composite fermion theory. 

It has become clear that a more physical way of looking at the composite particles 
is as bound states of one particle and <j> vortices. The composites are literally dipoles. 
Since the quantum Hall liquid cannot be reached perturbatively from a normal electron 
fluid, Landau's Fermi liquid theory cannot be applied, and the new quasiparticle must be 
built up from scratch. This notion developed steadily starting from Laughlin's observation 
that particles see other particles as flux Gfl. Jain used it to obtain the basic sequence of 



fractional quantum Hall plateaus [13|, and Haldane ]14|] and Halperin [gj constructed a 
complete hierarchy of states. Their approach relied solely on the analytic properties of 
ground state wavefunctions. The corresponding field theoretic implementation is the class 
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of CS models outlined above. In spite of the successes of these descriptions, neither does 
justice to the particle-vortex composites as bound states (dipoles) in their own right. 

This thesis will clarify that perhaps this is the essential ingredient required for a self- 
consistent understanding of the effective mass and compressibility issues. Several authors 
have applied the dipole scenario with some success [|TJ, [l5|, [l6|, 17, 18, 19, 20]. The works 
of Lee [20 1 and of Shankar and Murthy |l5|, [0], |l^] reconsider the Chern-Simons action and 
recover some of the dipole physics. Here, we will approach the problem from the opposite 
direction by working in the LLL at the outset without any singular flux attachment. The 
language no longer includes (5-flux tubes or <j>, rather particle-vortex dipoles are the basic 
building blocks. Our guide is a formalism introduced by Haldane and Pasquier [19| and 



applied by Read [18]. These authors considered bosons at v = 1, which is presumably qual- 
itatively identical to fermions at even denominator filling fractions. Both the effective mass 
and compressibility puzzles inherent in the original Chern-Simons approach are resolved in 
this way (at least for bosons at v = 1). 

Our first task, in the following chapter, is to extend the Haldane-Pasquier formalism to 
arbitrary v for either fermions or bosons. In the case of bosons, Chapter ||, we will obtain an 
effective theory microscopically by following Read's analysis at v = 1. The effective action 
does lead to a consistent picture of the mass and compressibility, but differs fundamentally 
from the CS action; there is no a priori reason the two actions ought to look similar since 
we project to the LLL from the beginning. However, it is gratifying that some of our results 
overlap with those of Shankar and Murthy, who do start with CS. 

In Chapter ^, we construct a phenomenological Landau-Ginzburg field theory for v = 
1/p, i.e. the underlying particles must be fermions for p =odd and bosons for p =even. The 
spectrum in this case contains a magnetoroton dip, which is the first analytical observation 
of this phenomenon. 

A central theme running through both of these chapters is the internal structure of the 
composite particles. To compensate for the extra degrees of freedom of the vortices, a set 
of constraints is introduced and appears thorughout our models. 
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1.2 Pairing in Two Dimensions 



The standard theory for superconductivity was introduced by Bardeen, Cooper, and Schri- 



effer (BCS) almost fifty years ago [21]. The original ground state was a paired state of 
fermions in the s-wave channel, or relative angular momentum I = 0. Since then, there have 
been numerous generalizations to non-zero angular momentum and to other more compli- 



cated order parameters. For a review, of the rich phase diagrams in He 3 , see the book by 



Vollhardt and Wolfle [22]. We refer to the non-zero angular momentum states as anoma- 
lous in the sense that they violate both parity and time reversal symmetry. Typically, BCS 
theory is applied to fermionic particles, which can be thought of as pairing into bosons 



which then condense. A less widely appreciated body of work treats pairing of bosons [23]. 
We explore instances of both statistics in Chapters |3] and |5[ There is growing evidence of 
anisotropic superconductivity in the perovskite oxides. In Sr2Ru04, both experiment 



and theory [25] support p-wave (I = —1) pairing. Similarly, d-wave pairing (I = —2) is by 



now well established in the high temperature superconductors [26]. p-wave pairing has also 
been observed in superfluid He 3 in the so-called "A-phase" |22j . 



In the context of the fractional quantum Hall effect, Halperin [27] proposed that under 
certain conditions, electrons can form pairs that condense into a Laughlin state of charge-2 
bosons. Since then, various alternatives for p- and d-wave pairing have been explored by 



several groups [g8|, [3C], 31], |33|, gj, |35[. An intriguing possibility is that the observed 



plateau at v = 5/2 is the Pfaffian state of Moore and Read |28|, 33, pq ]. 

Recently it has been proposed that the spin conductivity of the class of p- and d-wave 
states is transverse and quantized f37j . A remarkable series of earlier papers by Volovik 
contained some of these predictions, as well. However, the precise proof of this proposal has 
not been shown. Section [5] contains our derivation of this effect with a conserving approxi- 
mation. In the context of the original BCS theory, a conserving approximation is required 



for a correct description of the collective mode 21, |3^] since it respects charge conservation. 
Analogously, in our case, the conserving approximation will respect the conservation of spin 
current, leading to the correct result. 

A less familiar application of BCS theory is to paired states of bosons (see for example 
Typically, there is a competition between the usual single-particle condensate and a 
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pure paired state, or a phase of coexistence is possible. 

Pure p-wave pairing of spin-1/2 bosons is characterized by a BCS wavefunction in the 
form of a permanent. We will show that p-wave pairing can be interpreted as a condensate 
of spin waves. As parameters in the Hamiltonian are tuned, one reaches another single 
particle condensate with helical spin order. The permanent sits right on the transition and 
contains a single anti-Skyrmion, which is yet another single particle condensate. 

In the context of the fractional quantum Hall effect (FQHE), the permanent describes 



singlet pairs of spin-1/2 composite bosons with filling factor v = 1/p, p odd [28, 30|. It is 
the unique ground state wavefunction of a Hamiltonian which penalizes the closest approach 
of three spin-1/2 fermions for a fixed number of flux quanta piercing the bulk. The exact 
number of flux threading the bulk depends upon the geometry, e.g. on whether the system 
is on the torus, on the sphere, or on the plane. Although the permanent itself is difficult 
to treat analytically, its Hamiltonian contains an infinite set of degenerate zero-energy 
eigenstates when flux is added, among which is the polarized Laughlin state, and the rest 
are interpreted in one of two equivalent ways as either quasiholes or as spin wave excitations. 
The p = 1 Laughlin state is most amenable to analysis as it is a Slater determinant of single 
particle wavefunctions (Section |3.1.c ). Accordingly, it will serve as the prototype for our 
exact statements. 

We also consider d-wave pairing of spinless bosons, which is known as the Haffnian 



in FQHE literature [32]. A rich phase diagram of competing single particle and paired 
condensates emerges, with the Haffnian sitting on a phase boundary. 

In the case of pairing in the FQHE, the bosons (fermions) are to be thought of as 
composite bosons (fermions), as described in Chapter [| We can use CS mean field theory 
such that the CS gauge field cancels the external field on the average, and we are left with 
composite particles in zero net field. Then we can treat the composite particles within a 
BCS approximation. The case of composite fermions, which was addressed recently [35], 



serves as a point of reference for composite boson pairing. The central theme in ref. [35] 
was the existence of two regimes, weak- and strong-pairing. The weak-pairing phase can be 
characterized by a nontrivial topological winding of the BCS order parameter in momentum 
space; we will see an example in the discussion of the quantum Hall effect for spin in Chapter 
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||. On the other hand, the strong-pairing phase is topologically trivial, and, in the simplest 
case, the two phases are separated by a transition at zero chemical potential. For example, 



the Pfaffian state [28] is a weak-pairing phase, while the Haldane-Rezayi state f29|] , which 
is an I = —2 state, was found to lie at the weak-strong transition point. It is useful to keep 
these results in mind as we consider composite boson pairing in this chapter. 

1.3 Adsorption on Carbon Nanotubes 

Monolayer adsorption of noble gases onto graphite sheets has proven to be an interesting 
problem both theoretically and experimentally [40, 41, 42]. Many of the observed features 



can be understood within a lattice gas model, where the underlying hexagonal substrate 
layer forms a triangular lattice of preferred adsorption sites. An equivalent formulation 
is in the language of spin models on a triangular lattice, where the repulsion between 
adsorbed atoms in neighboring sites translates into an antiferromagnetic Ising coupling. 
The frustration of the couplings by the triangular lattice leads to the rich phase diagram of 
the monolayer adsorption problem pi]] . Introducing hopping adds quantum fluctuations, 
further enriching the phase diagram f42] , [43|| . 

In this Chapter we address what happens if, in addition to the triangular lattice frustra- 
tion, one has an extra geometric frustration due to periodic boundary conditions. In fact, 
such a system is physically realized by a single walled carbon nanotube [[44|], which may be 
viewed as a rolled graphite sheet. In this context, adsorption has been the subject of grow- 
ing experimental and theoretical interest [^5], |4(| spurred by potential applications. Stan 
and Cole jS| have considered the limit of non-interacting adatoms at low density, finding 
that they are localized radially near a nanotube's surface at a distance comparable to that 
in flat graphite (~ 3A). In that work, it was sufficient to omit the hexagonal structure of 
the substrate. However, the corrugation potential selects the hexagon centers as additional 



commensurate localization points [4C[. In view of the similarity to flat graphite, we include 



both the substrate lattice and adatom interactions and consider a wider range of densities. 



In fact, very recently, it has been shown [47] that the adsorbate stays within a cylindrical 
shell for fillings less than ~ 0.1/A 2 (or ~ 0.5 adatom/hexagon), justifying the densities 
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studied here. 

Our adsorption model is equivalent to a new sort of XXZ Heisenberg quantum spin 
tube, which is type of spin ladder with periodic boundary conditions. A simple example 
with highly anisotropic couplings was considered recently in references Hfl. We find density 
plateau structures for armchair, zig-zag and chiral nanotubes. In the language of spin 
systems, the density plateaus correspond to magnetization plateaus. The zig-zag tubes 
turn out to be special, and have extensive zero temperature entropy plateaus in the classical 
limit. Quantum effects lift the degeneracy, leaving gapless excitations described by a c = 1 
conformal field theory with compactification radius quantized by the tube circumference. 
This is an interesting conformal symmetry because the only other systems in nature with 
a quantized compactification radius, that we are aware of, are the chiral edge states in the 
FQHE m. 



17 



Chapter 2 

Lowest Landau Level I: Composite 
Fermions 



We begin this chapter by developing the Haldane-Pasquier formalism in Section 2.1. 

In the rest of the chapter we apply the composite fermion formalism to bosons in the 
lowest Landau level. In particular, we obtain an effective theory for an incompressible 
quantum Hall liquid of bosons with one attached vortex at general filling. As discussed in 



the introduction, theories based on the flux attachment, or Chern-Simons approach [10], 
have not been entirely satisfactory. More recently there have been several attempts to 
obtain an effective theory microscopically |l^, [H], ^0|, which too have had difficulties. In 
this work we avoid the Chern-Simons approach and follow an alternative that was developed 



for the Fermi liquid-like state of bosons at v = 1 [18|. 

The effective filling factor in our composite fermion model can be obtained as follows. 
As we will see below, the composite feels an effective magnetic field B = Bl + Br, where 
Bl is the physical field felt by the underlying particle and Br is an arbitrary field felt by 
the vortex. The total charge of the composite is q* = 1 + Br/B^. In the introduction we 
showed that the filling factor can be written as v = p<&q/Bl, where p is the density of the 
underlying particles and <3?o is the magnetic flux quantum. Since the composite particles 
must have the same density as the underlying particles, the effective filling is v e s = p&o/B, 
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which can be rewritten as 



The original Jain construction [13] begins with particles of charge +1 at v = 1/p and 
attaches <p flux tubes. On the average, the flux tubes renormalize the real magnetic field by 
Bl — ► q*Bi with q* = 1 — (pu. In our framework, q* can be any real number, so there is a 
family of theories, i.e. anyons fl50| , for any given v parameterized by q* and an integer v c q . 

Strictly speaking, Jain's model does not require fermionic statistics for the underlying 
particles, so the bosonic case is a valid quantum Hall state. In fact, there have been recent 
theoretical proposals that quantum Hall liquids of bosons are realizable in rotating Bose- 
Einstein condensates ]|51|j. In this case, the magnetic field is due to the angular velocity 
u; roughly, the velocity is modified in the rotating frame by v — > v + to x r, which is like 
minimal coupling of a magnetic field (V x (lu x r) is a constant in the same direction as to). 

In section |2.2| we obtain an effective mass by extracting a kinetic term from the interac- 



tion Hamiltonian. This problem has been central to the recent work in references [15, 19]. 
The simplest starting point is a Hartree-Fock approximation. Next, we use a conserving ap- 
proximation, which restores the constraints, in order to calculate the correct density-density 
response function. The self-consistent diagramatics consist of summing ring and ladder di- 
agrams; it is essentially the same framework that we use for paired states of fermions in 
Chapter |B| of this thesis. 

2.1 Formalism 

2.1. a Single Particle in a Magnetic Field 

A convenient framework for a quantum particle in two dimensions is an operator description. 
We begin with the simplest case of a single charged particle in two dimensions and a 
perpendicular magnetic field. 
The Hamiltonian is 
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where p is the canonical momentum, A is the vector potential for a magnetic field B, in 
the z direction, and q and m are the charge and mass. The units are set to fi = c = 1. 
The kinetic momentum is defined by 

7T = p - qA (2.3) 

with the corresponding commutator 

[TTf 1 ,TT u ] = ie^qB (2.4) 

where \x and v are space indices, x and y, and e^ v is the Levi-Civita symbol. The dynamics 
of 7r follow simply, 

TT^= i[H,7T,j] = UJce^TTy , (2.5) 

where uj c = qB/m is the cyclotron frequency. Therefore ir precesses. 

As the particle executes cyclotron motion it is located by the guiding center operator 

R = r + zxvr— (2.6) 
qB 

which obeys 

[R fl ,R v ] = -ie fll/ -^ . (2.7) 

The operators tt and R commute. Projection to the lowest Landau level is accomplished 
by replacing the particle's coordinates by the guiding center. Note that the coordinates no 
longer commute, a consequence of frozen degrees of freedom. The appropriate generator of 
translations, or "pseudomomentum" , is defined by 

K = qBz x R , (2.8) 

and obeys [K^, K v \ = ie^ v -^. The planar coordinates of a quantum particle in a magnetic 



field are a well-known example of non-commutative space [91]. 

Noting the commutation relations of tt and R, we can define two independent harmonic 
oscillator operators: 



— K , a) = J—K (2.9) 
2qB V 2qB V ' 



1 , 1 



7T . 



2qB V 2 Q B 



20 



tf = J——w (2.10) 



where tt = ir x + iir y and tt = ir x — iir y , and similarly for K. The Hamiltonian then becomes 
the familiar harmonic oscillator Hamiltonian: 

H = u c (b*b+^) . (2.11) 

Its eigenstates are known as Landau levels (LL). Each level is macroscopically degenerate 
since a commutes with b and drops out of the Hamiltonian. The complete set of eigenstates 
is labeled by two integers, m and n: 

\n,m) = _=_|0,0) (2.12) 
yrrv. yn\ 

In this convention, tt (or b) is a purely inter-Landau level operator, and K (or a) is intra- 
Landau level. 

To obtain the wavefunctions, we will use complex coordinates, r = z = x + iy, and the 
symmetric gauge, A = — x B = \B(— y,x). The single particle operators become 

tt = -2idz-^-iz (2.13) 
K = -2idz + ^-iz (2.14) 

If we restrict ourselves to the lowest Landau level (LLL), then the wavefunctions, Y>o,m(r), 
are annihilated by b. The general solution to this first order partial differential equation 
(up to gauge transformations) is 

^o,m(r) = fm(z)e- lzl2/ < , (2.15) 

where i 2 B = l/\qB\ is the magnetic length. The great simplification is that f m (z) must be 
an analytic function in z. Further requiring that ai/'o.o = yields /o and the rest of the 
fmS are generated by applying a}. The result is the set of states spanning the LLL: 

Um(z) = . 1 z m e -\z\'/4fi B ^ lg) 

v /27r2 m m!^™ +2 

It should be remarked that the intra-LL ladder operator has a very simple action on the 



u m : a) =- : 



fy/2t B and a = (£ B /V2)d z . 
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2.1.b Two Particles in a Magnetic Field 

In this section, we introduce a bound state of two oppositely charged particles in a perpen- 
dicular magnetic field. This formalism will be useful in interpreting the Haldane-Pasquier 
approach. 

Each component is characterized by its charge % and the magnetic field that it feels, 
B{. For convenience, we define 

B % = Ql B (2.17) 

We will assume that B\ > and B2 < 0, guaranteeing the existence of a bound state. In 
these units the charge is dimensionless. If we fix q\ = 1 then the total charge is 

1* = (2-18) 
B\ 

There are two sets of guiding centers, R^, and pseudomomenta, K^, which are defined in 



the same way as in section 2.1. a. The algebra consists of two copies of the single particle, 
for example 

[RuijRi/j] — ^iiv^ij~^r 1 (2.19) 
Bi 

and so on. 

It turns out that this algebra can be mapped exactly into a single particle in an effective 



magnetic field B p2|, 



B = B 1 + B 2 . 

If we construct the effective translation and momentum operators by 

K = K l + K 2 (2.20) 
1 „ 1 



vr = ^B^B 2 (—K 2 -—K 1 ) (2.21) 



then the commutators of this algebra are exactly that of a single particle, equations (£ 
2.7 , 2.8). For example [vr^,^] = ie^B. The physical picture becomes clearer once we 



define an effective position r by eqn. (2J3). Solving for r, we find 



B -(Ri-R 2 ) V B • (2.22) 
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In the limit B — > 0, r becomes (Ri + R2)/2. This is not surprising from a classical point 
of view; two opposite but equal charges travel in a straight line due to E x B drift, the 
"guiding center" moves off to infinity, and their position is given by a point exactly midway 
between them. 

To proceed with the wavefunctions of the composite, label the real-space coordinates of 
each particle by z and rj. The differential operators are given by 



Ki 

K 2 



-2id~ + 



(2.23) 
(2.24) 



B-2 



Note the relative minus sign, which comes from assuming B\ > and B 2 < 0. The ladder 
operators are given in terms of ir and K by equation ( |2.10| ). The lowest eigenfuction is 
determined from atpofi = 6^o o = 0: 



/2B 



-28 z 



2K 



z + 



-2d v + ^r, 



/2B 



i 

B 2 



-2^-^T? 



■00,0 = 

V>o,o = 



(2.25) 



By analogy with a single particle, we expect that the solution is an analytic function in z 
and 7] times two Gaussian factors. Indeed, if we choose 

then equations ( |2.25| ) are solved by 

<K*, V) = 9 /, ^ /2t% > ■ (2-26) 

Note the asymmetry between B\ <-> B 2 , stemming from the sign of B. We have implicitly 
assumed that B > 0, but if B < then the Gaussian factor in eqn. ( p. 26 ) would be £ B . 

The complete set of states is generated by a) and &t just as it was for one particle in 
eqn. 



A! V/d 

The ipXfj, are linear combinations of the independent particle basis 



(z,r)\mn) = u m (z)v n (rj) 



(2.27) 



(2.28) 
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where u and v are LLL single particle states corresponding to the two magnetic lengths Ib 1 
and Ib 2 i as in eqn. ( 2.1 €| ) . 



2.1.c Fock Space, Operators, and Constraints 

Having constructed the basis functions for a particle-vortex pair in the previous section, we 
move on to the Fock space for a many-particle system. 

We begin with canonical fermionic or bosonic operators which are matrices with two 
indices, c mn , with (anti-)commutation relations 

[ c mni c n>m']± = fimm'finri (2.29) 

The left index, m, runs from 1 to N^, the number of available states in the LLL. The right 
index, re, runs 1 through N, which we will interpret later as the number of vortices. In the 
thermodynamic limit, the filling factor is v = N/N^. Strictly speaking, this construction 
must be carried out on a finite geometry (e.g. m runs from 1 to Nj, + 1 on the sphere), 
but we will ignore this subtlety here since we will only be interested in the thermodynamic 
limit. 

The anticommutation relations are invariant under independent transformations on the 
left and right indices: 

c^U L cU R , (2.30) 

where Ul and Ur are x and N x N unitary matrices. These transformations are 
generated by the left and right "density" operators 

Pnn' = c nm c mn' 
m=l 
N 

Pmm' = c nm c m'n (2-31) 
n=l 

The left density pr will represent the physical density, as we will see below. The right 
density p R specifies a set of iV 2 constraints, which we use to define a set of physical states, 

(pL> ~ <W) l^phys) = . (2.32) 
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It will be shown shortly that the set of l^phys) do indeed give the correct Fock space. 

Because the p R generate the unitary group XJ(N)r and there is a phase factor, U(l), 
common to both p R and p , the physical states must be singlets under SU(./V)r. The 
physical states solving this constraint are linear combinations of 

l^phys ' e c n 1 m 1 c n 2 m 2 c n N m N \ (J / ; ^Z.OOj 

ni---njv 

where |0) is the vacuum with no fermions. The Levi-Civita symbol e ensures that these 
states are singlets under SU(A")r. If the c's are fermions their anticommutation relations 
ensure that | ^^h ys mjv ) is symmetric under the interchange of any pair rrii «-> raj. Therefore, 
the physical space is equivalent to N bosons each of which can occupy any one of Nm states, 
i.e. the Fock space of bosons at filling v = N/N^. Had the c's been bosonic operators rather 
than fermionic, the result would have been a Fock space of fermions at the same filling. 

We now construct the many body operators. By mapping into a single particle, we 
showed that the effective magnetic length is 

111 

JT ~ ~p2 1T~ (2.34) 

As a reminder of the physical picture, we have made the notation change from 1,2 to L, R. 
Our convention guarantees that £b is positive since £b l < £b r (equivalently, \Bl\ > \Br\ 
and B = B L + B R > 0). 

In real space the matter field is defined by 

C (Z,V) =J2 u m( z )u R (v) Cmn (2.35) 
mn 

where u m (z) = z m exp{ — \z\ 2 /A£ Bl } and u^rj) = z m ex.p{ — \r/\ 2 /4£ Br } (apart from normal- 
izations). A unitary transformation connects the \mn) independent particle basis to the 
| p\) bound state basis in section 2.1.b| . Accordingly, the fermions c mn transform into 



Ca\ = ^2c mn (mn\p\) , (2.36) 

mn 

where 

| M A) = -=-=|0,0>. (2.37) 
y/p ] - VA! 
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The overlap {mn\p\) is obtained from the definitions in eqns. fl2~27l ) and ( ggf ). To write the 



density operators, we observe that they take a plane wave form in the operator language: 
Pq = J2i e" 1 ' 11 *'^ and p^ = J2i e iq ' Ri ' L , where Ri,_R and Hi,L are the two guiding center 
coordinates of the i'th particle. [] In second quantization, the left density becomes 

Pq= E ^VA^AIe^l/ZA') (2.38) 

and similarly for p^ . The matrix element can be calculated by solving for Rj^l in terms of 
K, 7T, giving 

R L = A (K+^7r)4 

R K = A(K + ^vr)4. (2.39) 

We have introduced the shorthand notation, Aa = —z x a (for the vector a). Because 
[K,7r] = 0, the plane wave factors into an intra- and an inter-Landau level piece. We write 

Pq = E P<lW)plW)A^'\' 

Pn = E P*W)f$W)<i*W > ( 2 - 4 °) 
where the p-coefficients are defined by 

Pq(/V) = Ou|exp{^|qAK}|//) 

pL(\\\') = (A|ex P |^|^qA7r}|A') (2.41) 



pJ(A|A') = (A|exp|rf|^qA^}|A' 



The operation a A b stands for a • Ab. Note that /5 q is identical to f q for a particle in field 
B and therefore follow the same orthonormality properties and commutation relations. Pq 
and Pq are similar but with an additional factor of £b l /^b r or £b r /£b l in the phase of the 
commutator. 

Since K and tt are nothing other than harmonic oscillator operators, the p's can be 

calculated explicitly. The one that we need later is p L , so we use it as an example. First, 
*I was reminded by R. Shankar that he had guessed the same density expressions based on a small q 
limit of the Chern-Simons formulation 
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rewrite A ir in complex coordinates as ^(^7r — qW). Secondly, recall that b ~ ir and 
| A) ~ ^ A |0), which reduces the calculation to a harmonic oscillator matrix element. The 
rest is straightforward and we quote the final result: 

(A|Y) = —^=—j== e^ 2 {-&i)\d q f e"« . (2.42) 

A bunch of factors were absorbed into q on the left-hand side to avoid repetitiously writing 
them on the right. 

The final step in the construction is the Hamiltonian. In coordinate space it involves 



only the diagonal components of p 18 



H = -j d 2 Zl d 2 z 2 V(n - r 2 ) : p L (z 1 ,z 1 )p L (z 2 ,z 2 ) : (2.43) 
Or in Fourier space (Appendix), 

where V(q) is the ordinary Fourier transform of V(r). As required, this Hamiltonian is 
both translationally and rotationally invariant. By construction, the right density p R is a 
constant of the motion because the Hamiltonian acts only on the left indices, i.e. 

[H, p%] = . (2.45) 

The system of Hamiltonian plus constraints is the starting point for Read's analysis of 



composite bosons at v = 1 [18]. 



2.2 Hartree-Fock Approximation 

For convenience, we restate here the Hamiltonian that we derived in the last section. 
1 f (i^q 

^=2^ J 72^ F (q) F q(^i A i'/ i 3 A 3|^2A2,/x4A 4 )4 iMi 4 3M3 c At4A4 c M2A2 , (2.46) 
where the matrix element is given by 

^q(/"lAl,/X3A 3 | / U2A2,/i4A4)=Pq(/Ul|/i2)/0-q(^3|^4)/5q(Al|A2)p^ q (A 3 |A4) (2.47) 
~ — Iql 2 /2£ 2 

and V(q) = V(q)e b l is the apodized potential. The vertex is shown in Fig. |2.1| 
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Figure 2.1: The vertex p q (/zi| /i2)p q (Ai|A2). The dotted line represents the interaction V"(q) 

The Hartree-Fock (HF) approximation consists of replacing pairs of fermion operators 
by their expectation value at zero temperature: 

( c h\ c ij-'x)o = <W^AA' 0(A m ax - A) , (2.48) 

which fills X m ax + 1 Landau levels of composite bosons. In Jain's mapping p^], this corre- 
sponds to an effective filling fraction v e s = X ma x + 1. At finite temperature, the 0-function 
is replaced by the Fermi distribution f(e\ — /u c ), where fj, c is the chemical potential. Since 
the chemical potential is restricted to e\ max < \i c < £\ ma x+i but is otherwise arbitrary, we 
will drop it in the following. 

Expanding the Hamiltonian around this ground state allows us to sum over the intra- 
level indices /ij. In the following, we ignore the direct term, which only shifts the chemical 
potential. In the exchange term, completeness of the basis within a Landau level, 
= 1, gives or S^ 2/X3 . Completeness of the intra-LL basis is tantamount to 

translation invariance. Rotation invariance, on the other hand, shows up in V(q), which is 
required to be a function of only |q| for an isotropic system. This gives terms diagonal in 
A, <5aiA 4 or 5\ 2 \ 3 , since the others vanish when we consider the explicit expression for f L in 
equation ( |2.42| ). The HF Hamiltonian is now 

#o = E £ a4^a (2.49) 

with the exchange energy 




(2.50) 
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For concreteness, we choose the simplest non-trivial case: we fill only the LLL, X ma x = 0, 
and take a hard-core repulsive interaction, V(q) = V(0). Using the explicit form of from 
equation ( |2.42|) the energy becomes 



2A 



g^j y(0)p, (2.51) 

where p = 1/2tt£ b is the density of the composite fermions. Note how the apodized potential 
now has £ B , not t\ , in the Gaussian, which is a consequence of the relation 1/1% = 
1/£ Br . Because £ Bh < £% R , the energy vanishes asymptotically as A — » oo. Furthermore, the 
exponential form implies a linear dependence on A at small A, exactly the kind of behavior 
that one would expect for non-interacting particles in a magnetic field. The cyclotron 
frequency, from equation (2.11), is given by uj c = l/m£ B so we can identify an effective 
mass with the gap A by \/tn* = (si — EQ^j£ B . More generally when Xrnax 

> 0, the low 

energy physics is dominated by transitions between the highest occupied LL and the lowest 



unoccupied one, which gives 

1 



rrr 



(ex 

m ax 

+i - £a 

rn ax ) 



A . (2.52) 



Further justification for this interpretation will emerge as we consider fluctuations around 
the ground state. At any rate, our calculation provides a framework to calculate m* in the 



LLL. Note that equation (2.50) shows that m* has contributions from the lowest LL's (of 
composite particles) and identifies it with a particular integral over the interaction. 

As a stand-alone approximation, HF does not preserve the constraints because the 
commutator [Hq, p^] ^ for all q / 0. In the next section we augment HF in a fully 
self-consistent manner to restore this symmetry. 



2.3 Conserving Approximation and the Ward Identity 

Our goal is to find a perturbative series such that all correlation function that involve 
the constraint p^ — /5<5q,o vanish. In other words, the constraints would vanish to any 

29 



order of approximation. This is guaranteed by an exact (non-perturbative) Ward identity, 
which is derived below. We illustrate the method by calculating the response functions, or 
generalized susceptibilities, of p R — p R , p R — p L , and p L — p L . The latter is related to the 
physical quantity of interest, the compressibility. 



It is well known in the theory of metals [21, p3j that if a Hartree-Fock approximation is 



used for the two-particle Green's function, then a fully self-consistent approximation that 
conserves charge includes ladder and bubble diagrams in the response functions. In fact 



this method preserves the constraints in the composite boson problem at v = 1 as well [18|, 
and we will show that it works here, too. 

In imaginary time, the response functions take the form 

X AB (q, i0J n ) = (^(iiO^p^i-iiOn)) , (2.53) 

where A, B stand for R or L and uo n = 2nir/(3 are Matsubara frequencies. We implicitly 
keep only the connected part, thus dropping a term containing (p )'s. The fundamental 
diagrams are those that are irreducible, i.e. those that cannot be separated by cutting an 
interaction line. For a short range interaction, these are the qualitatively relevant pieces 
pl| , |53| ) so we will not perform the bubble sums explicitly here. In any case, they are easily 
obtained as geometric series of the irreducible parts (l^, 53]. 



The form of the conserving approximation in our case states that the irreducible response 
functions, x^f, are to be calculated by including the ladder series with the HF Green's 
function lines. We sum the series by solving Dyson's equation; the next few equations will 
describe the structure of the theory. 



First, let us recall the HF Green's function 53], 



Go(x,wv) = -. 1 V ( 2 - 54 ) 



^ = -^EE/(^^(q)Pq( A l A, )p L q (A'|A)g(A',^) 

= " E / (S^pfoO l^( A l A ')l 2 fey - Mc) , (2.55) 



where lo v = (2u + l)vr//3 is a fermionic Matsubara frequency, and f(£y — Pc) is the Fermi 
distribution with respect to the LL index (p c and v should not be confused with the intra-LL 
index and the filling fraction). 
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Second, the renormalization of the vertices, A A , by the ladder series can be written as 
a matrix Dyson equation: 

Aj^yfa^n) =p£(/iA|//A')- (2-56) 
dp k 

~ 53 / ^^^( k ) F k(^A,M2A 2 |/iiAi,/i / A / )pA 1 A 2 (^n)A^ lAliM2A2 (q,iw ri ) 

where has been defined in equation (|2.47| ) and T> is the frequency sum over the internal 
Green's functions, 

T> Xl \ 2 (iuj n ) = -^Y^QQi^ioJy + iw n )Go(Xi, (2.57) 

/(£A 2 ~ He) - /(gAi ~ Pc) ^ 

£a 2 — £ Ai — iu n 

Because the Green's function is independent of the intra-LL index fj,, it is convenient to 
define a purely inter-LL vertex, A by 

A M,^'A' (q, = Pq(/V) M\' (q> W n ) ( 2 -59) 

so that Dyson's equation becomes 

A A 4 A ,(q, i ^) = ^(A|A')- (2.60) 



d 2 k 



F(k)p^(A|A 1 )p z : k (A 2 |AOPA 1 A 2 (^n)e iqAk£2 fl A< 1 lA2 ( q ,^ ri ). 



(2vr)2 

The phase factor is due to the magnetic translation commutator algebra — see the discussion 
immediately following equation fl2.42| ). 

The response functions are given in terms of the renormalized vertices by 

x4 B (q>^n) = -/9o^A A 1 lA2 (q,iu; n )pA 1 A 2 (^n)p i? q (A 2 |Ai) , (2.61) 

A, 

where po = l/2iii 2 B is the density of particles per LL, coming from the trace over fx's. The 



diagramatic structure is shown schematically in fig. 12^2 

The key to showing that all response function containing A^ vanish, i.e xS^ = Xh^ = 
= 0, is the Ward identity for A R , which we now derive. The first principles derivation 
follows standard field theoretic techniques pT|. Consider the exact vertex in real time: 



AjU ll/l2 A a (q,Mi,t2) = (^{Pq(*)c m Ax^i)4 2M2 (*2)}) , (2-62) 
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Figure 2.2: Diagramatics of the conserving approximation, (a) The exchange self energy; (b) the 
ladder series renormalization of the vertices, A; and (c) the response functions in terms of A. 



where T is the time ordering symbol. By taking the time derivative dt of both sides, using 
dtfjq = 0, and then Fourier transforming back to frequency space, we find the exact Ward 
identity 

icj n Af lA2 (q,iw n ) = g~ 1 (\ 1 ,iu; n +iLJ l/ )pR(\i\\2) - p% (Ai|A 2 )^ _1 (A 2 , iu v ) . (2.63) 

Here, Q is the exact Green's function, and iuj v on the right-hand side cancels identically, but 
is introduced for convenience. The two terms are due to differentiating the time ordering; 
physically, they are due to c\ti) and c^) acting as sources in equation (2.62). For the 



particular HF and ladder series that we use here, we can verify the Ward identity by 
substituting Q$ for the exact Q and plugging the whole expression into the right-hand side 
of the Dyson equation ( 2.60| ). Upon using the definition of Qq from equation ( 2.55| ), we find 



that the ladder series satisfies the Ward identity. 

Therefore our diagramatic scheme preserves the constraints, and we can be sure that 
the physical quantities that we calculate in this approximation will be consistent. 

2.4 Response Functions 



It is now straightforward to use the Ward identity in the response functions, equation ( 2. 61] ), 
to verify that 



X £f (q, iu n ) = xSr (q, »wn) = Xtr (q. = (2.64) 
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and in fact for any correlator containing A R . There is one proviso in this procedure, that 
is discussed in detail in reference [p!8[ , having to do with the constraints at iui n = 0. Our 
procedure is only valid at non-zero frequencies because we divided the Ward identity by 
iuj n to isolate A^. A complete proof requires more care, but we will not pursue this here. 
In any case, there is no problem with taking the limit to — > 0, which requires only small but 
non-zero frequencies. 

It remains to calculate the physical density-density response, x LL ■ To this end, we first 
rewrite the vertex in a more symmetric fashion by introducing the scattering matrix, T, for 
the ladder series. Although the internal Green's functions do not carry momentum, which 
comes in only through the vertices, we can nonetheless absorb some of the momentum 
dependence into T by taking advantage of translation invariance again. We define 

fAiA' 1 ;A 2 A' 2 (q>^r l )=Xl PlO^l I Ai,/4 A^; M2A2, M 2 A 2 (^V-qC/^M) ( 2 -65) 

The Dyson equation for scattering, also known as the Bethe-Salpeter equation, takes the 
form 

fA 1 A' 1 ;A 2 A^(q,^n)= / ^V^A ( Ai I A 2 )*£ k ( A' 2 1 \[ ) e^l - (2.66) 
"E / 7 ^^(k)^(A|A 2 )p^ k (A' 2 |AOe^ k ^P A v(^ n )f AlVi;AV (q,i W „). 

AA' ^> 

It is convenient to view r AlA ' . \x' as a vector with components labeled by AA', while AiA^ 
and (q, ioo n ) are parameters. Then the problem reduces to inverting a matrix in the indices 
(A 2 A' 2 ; AA'). 

The formal structure of this matrix equation is elucidated by reducing it to the eigenvalue 
equation, 

E M A 2 A' 2 ; AA' (q, i"n)A Xl ^ . AA' fa, = U Xl \[ (q, i^n ) A\i X[ ; A2 A 2 fa> iuJ n) , (2.67) 

AA' 

where the kernel is 

M A2A , ;AA Kq,^ = 5 AA2 A^ (2.68) 

and u, A are the eigenvalues, eigenvectors. We have obtained an exact zero eigenvalue 
solution of this equation at iuj n = 0. Using the properties of the commutatator [p k , p^] (c.f. 
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equation ( 2.42| ) and the immediately following discussion) and the definitions of e x (2.55) 
and of V\\i (|2.58|) , we find that 



A, 



AiAi; AA 



(q,0) = e x p*{\\>!)-pK{\\>!)e x , , 



^AiA'^Q) ) = °- 



(2.69) 
(2.70) 



The similarity of this solution to the Ward identity ( 2.63j ) suggests that the existence of the 
scattering zero mode is related to the vanishing of correlators containing p . 

Another advantage of the Bethe-Salpeter equation is that we can rewrite the response 
functions symmetrically, 



+Po Pi ( A i I A i K ( iuJ ™ ) r Ax a; ; a 2 y 2 (q, iv n )T>\ a x > 2 {iu n ) p B q ( A 2 1 A' 2 ) 

A^ 

where \n B is the bare bubble 



(2.71) 



X^(q,^n) = -poE'°q( A l A/ )^AA'(^n)p B q (A , |A) 



(2.72) 



AA' 



Fig. 2.2 illustrates the summation. 



X = 



Figure 2.3: (a) Ladder series for the scattering matrix r and (b) response functions in terms of Y. 

To obtain the momentum expansion of % '■> consider the expansion of p from its 
definition in eqn. ( |2.42| ), 



Pq(A|A') = S xx > + <A| iq Att |) + 0(q 2 ) . 



(2.73) 



The first term, diagonal in AA', cannot contribute to the response because the transition 
amplitude V xx > is purely an inter-LL operator. The expanded response becomes 
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X"(q,w„) = -a^(A , |qA7r|A)P A A'(^n)(A|qA^|A / )+ (2.74) 
AA' 

+a^(Ai|q Avr|Ai)P Al v i (itJ n )f Al v i;A2 ^(q,itJ n )P A2A / 2 (ia; n )(A2|q AvrlA's) + . . . , 

where a is an overall constant. 

To obtain xtr through C(q 2 ), we need only T(0,iu> n ). At q = 0, the Bethe-Salpeter 
equation is 

f AiA' i; A 2 A 2 (0,iw„) = Vx 2 A 2 ;AiAi - X! ^A 2 A 2 ; AA'^AA' (^n)f AlA ' i; AA / (0, iw n ) , (2.75) 

AA' 

where the interaction matrix element is 

^a-a.a; = / 7^2^( k )/ 5 k(Ai|A 2 )p z : k (A / 2 |A / 1 ) . (2.76) 
Rotation invariance at q = requires that the matrix elements vanish unless 

Ai + A 2 = A'} + A2 . 

Now let us make two further restrictions that afford an exact solution for T. First, we 
choose Xmax = as we did above to illustrate the exchange energy. Secondly, we work 
at zero temperature where the Fermi function is f(e\ — /u) = Q(\ ma x — A), restricting the 
P-amplitude to 

1 



A A - lLO n 



V xo (iu n ) = - * (2.77) 

A A + lLO n 

A A = ex- e 



Along with rotation invariance, these restrictions allow us to solve equation ( 2.75| ) for the 
scattering matrix (at q = 0) 

VOX; OA 



OA; OA 



AO; AO 



1 + Vq\; 0aA)A 
^AO; AO 



(2.78) 



1 + V X o- AO^AO 

The two channels above represent a particle in the O'th LL propagating on one leg of the 
ladder diagram and a particle in the A'th LL on the other, and vice versa. These channels 

35 



do not mix in our example. In fact, we will only need A = 1 for the lowest order term in 
X LL because of the vertices (A|q A 7r| A'). At this point we find a crucial identity for A = 1: 

Vbl;01 =VlO;10 = A x , (2.79) 



which is easily proven by comparing equations ( p. 76 ) and ( |2.50| ). Plugging equations ( p. 78 ) 



( 2.79j ) into equation (2.75) for the response function, we find 



Xirr(q 5 ^n) = — I <0| vr| 1> I \V m +V m --—— \ \q\ z + 0(|q| 4 ) 

= + O(|q| 4 ) (2.80) 

(an overall factor has been left out). Thus, the lowest order term in the density-density 
response is of order |q| 4 . 

This is the main physical result that we wanted to reproduce within our composite 
fermion framework. Its main content is that the system is incompressible, i.e. the com- 
pressibility, k, vanishes. The connection of compressibility to the density-density response 



is contained in the definition [54] 



K = lim y iL (q,0) . (2.81) 

q^O 

For the irreducible diagrams that we have considered so far, \ LL = Xbr: so that k = 0. We 
expect that the order of limits is consistent with our calculation of the ladder series, which 
is a Taylor expansion around q = at finite frequency. 

The full response function \ LL can be obtained from the irreducible part by a bubble 



summation [21, 53 



LL _ Xi] 



LL = _ ( 2 g 2 ) 

i + v(cDxt L 



tirr 



For a short range interaction, such as ours, this geometric sum has no qualitative effect, so 



that k remains at zero. A well-known early work by Girvin, MacDonald, and Platzman [55] 
contains general arguments for the momentum dependence of the density response function 
for incompressible liquids in the LLL, and is consistent with our result. 

Before we close this section and move on to the effective field theory, the interaction-gap 
identity of equation ( 2.7S| ) is worth a couple more words. A slightly more general case is 
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when \ m ax > 0, 

Vx.A+l; A,A+1 = ^A+1,A;A+1,A = £A+1 — £\ , (2.83) 

where A = \ m ax- We expect that the compressibility will vanish again, although we have 
not performed this calculation explicitly. This identity seems to have an analog in the Fermi 
liquid- like state of bosons at v = 1 There, the system has a divergent compressibility 
due to a fixed Landau parameter, F\ = —1, and m* is also coming from an integral of 
the interaction below the Fermi surface. In both cases, the identities lead to the correct 
compressibility because there is no bare kinetic term due to the LLL projection. 

The form of the left and right density response functions, suggest a physical interpreta- 
tion of the vertices. To lowest order in q, all response functions vanish, and we can take any 
linear combination p^ — xp^ for the physical response without changing the result. Suppose 
we choose the weighted combination 

PS - Pi + fjrP? • (2-84) 

In operator language the densities are p^ = e* q ' R,A . Using the operator mapping in equations 
fljTgg ) and (|]22|), the momentum expansion of the new vertex leaves 



B 1 
q B L B L 

(1 + q • r) - -Liq A vr + 0(q 2 ) . (2.85) 



B ^ ,1 



The first term is the first order term of a plane wave for the composite particle with charge 
q* = B/Bl at position r, which is consistent with the two-particle mapping of section 2.1.b| . 



The second term is interpreted as the dipole moment. The charge of the composite does not 
show up in the response functions, but presumably would come out if backflow corrections 
are included as in the work of Lopez and Fradkin [^7]] . 
If the expansion is exponentiated, we obtain 

L _^ iq-r ( B 2 



Pq-e^ r ( — -^qAvr) . (2i 



This form agrees with the work of Shankar [15|, which starts from a different approach 



using the Chern-Simons theory at the outset. As the first line of eqn. ( ggp shows, this 
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particular choice of x makes the physical density a purely intra-LL operator. As such, it is 
obvious that the density-density response vanishes to 0(q 2 ) in the ladder apporximation, 
since the vertices contain only inter-LL transitions. Further, our weighted combination 
agrees with Read's |l8| density for v = 1 where Br/Bl = —1. Nonetheless, we stress that 
our derivation does not specify x by itself. 



2.5 Effective Theory 

In this section we will show that the ladder series can be replaced with a dynamic gauge 
field, yielding an effective theory much like that for bosons at v = 1. It is not a Chern- 
Simons field theory, but the familiar relation of density to the curl of a gauge field will 
appear. 

In accordance with the previous section, the low energy physics of bosons at v ^ 1 is 
that of composite fermions filling an integral number, z/ e fj = u/q* , of Landau levels. Let us 
take an ordinary fermion field c(x, t) in a static magnetic field B = V A A such that exactly 
v e g levels are filled, and couple it to a nondynamic gauge field a: 

H cS = f dtd 2 x-^\(-iV - A-a)c| 2 -/ic t c. (2.87) 

J ZTYl 

The fermion density is p = v e g/27r£ B , or equivalently p = v/2it1 2 Bl . The chemical potential 
p is tuned to lie between the uppermost filled LL and the lowest empty one. m* is the 
only parameter in the theory and is obtained from the smallest energy gap of the original 
problem as defined in equation ( |2.52| ). The action of the theory is an action for a as well 
as for c, S , but in contrast to Chern-Simons theory, there is no kinetic term for a; it is a 
"strongly coupled" gauge field in the language of field theory. There is also the term clqc^c, 
however we will choose the temporal gauge in which ao = at finite frequency, so that its 
fluctuations do not affect the response. 

The gauge symmetry of H e $ is ordinary U(l), which can be viewed as the long distance 
limit of the global U(A^)^ symmetry of the right coordinates. The gauge invariant density 
of this model, c^c, is identified with the constraint p R , which fixes 

p R = = p . (2.88) 



38 



This condition is the long distance limit of the full constraint that was constructed in 



Section 2.1.c . To obtain an expression for the physical density, consider the gauge invariant 
momentum density before a is included, 



g(r) = — {c f (V - iA)c- (V-HA)c f cj . 



(2.89) 



The single-particle version is the tt operator of Section 2.1.a, This allows us to rewrite the 



density suggestively. The first term in equation (2 



, g*e iq ' r , is a plane wave for a charge 



of magnitude q*; at tree level its expectation value can be replaced by p, which leaves 

p = p-<Z*4VAg, (2.90) 

where q* = B/B^. Similarly, at tree level the gauge potential is related to the momentum 
density by a = g/p. Therefore the physical density becomes 



P 



2^ VAa < 



(2.91) 



where we used u e g = v /q* . Equation (2.91) is precisely the fluctuation piece of the Chern- 
Simons equation ( |1.2| ), despite the absence of a kinetic term for a! It should be borne 
in mind, however, that we have imposed the particular linear combination p — xp R with 
x = B/Bl, which is responsible for this appealing result; in principle, any coefficient of 
V A a is obtainable in this way. 

Let us now consider the correlation functions. The basic conductivities are <x 
(gigj) /m* 2 . Because of the Onsager relation a xy = —a yx and isotropy a xx 



a 



fit/ ■ 



it is 



convenient to use complex coordinates g = g x + ig y and g = g x — ig y so that the expectation 
values (gg) and (gg) vanish. Our aim is to compare the susceptibility 

X LL (q) = (5p(q)5p(-q)) = £% L (q A g(q) ■ q A g(-g)) 



(2.92) 



to the ladder series in the previous section. We will show that within a random phase 
approximation (RPA), the responses are identical (at least to 0(q 2 )). 

The RPA has been applied in the context of the quantum Hall effect by several authors 
[18, |57| ]. It is a bubble sum for the gauge field fluctuations. The basic terms, shown in fig. 



2.4, consist of a diamagnetic and a bubble piece. The gauge field correlator is 
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2m* 



g 

2m* 



9 
2m* 



AAA/V 



(a) 



(b) 



Figure 2.4: Bubble summation for the gauge field propagator, (a) The shaded circle includes the 
diamagnetic coupling and the (gg) bubble; (b) The thick wavy line represents (aa) 



(a(q)a(-q)) 



(2.93) 



where {g(q)g(—q)) is the bare bubble, which we can evaluate in the single particle basis: 



(9(q)9(-q)) = - E (^\7r\^X')V X y(iu; n )(fi'X'\lf\^X) 



(2.94) 



The calculation of (a(q)a(—q)) is analogous, but with {g(q)g(—q)) ■ The matrix ele- 
ments of 7r produce a factor 2v c q1 2 b since it = \/2£b b is the inter-LL ladder operator and 
T>\\/ connects only states on opposite sides of the Fermi surface, which restricts A, A' to 
Xmax-, X m ax + 1- The sum over fi's gives the density per LL, l/2ir£ 2 B , with the end result 

1 



_ ~p A 



m* A + iuj n 

where one factor of m* has been replaced by 1/£ B A. 



(2.95) 



Now, the RPA g — g response consists of the bubble sum shown in fig. 2.5 



X = 



Figure 2.5: x LL as a bubble sum. The wavy line is the gauge field propagator from equation (2.93) 



See also fig. \IA 



(g(q)g(-q)) = {g(q)g(-q)) 



(2.96) 



2 (g(q)g(-q))o^ Hq)a(-q))^-; {g(q)g{-q)) 
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Using (gg) = (gg) and the identities ( |2.93| ), ( |2.95| ) in equation ( p.96| ) and plugging the result 
into the expression for x LL > equation (|2.92[) , we find that 

X LL (q)=0 + O(q 4 ) . (2.97) 

This is the same result as in the conserving approximation of the last section. The structure 
of the RPA is such that, the gauge field propagator replaces the scattering matrix T, and 
the two channels Toijoi, riO;io correspond to (aa), (aa). Diagramatically, the second bubble 



term in Fig. 2.5 is exactly the ladder sum in Fig. \2.S\ (b). 

Another way to test the effective theory is by integrating out the fermions. Since there 
are ordinary fermions filling f e g- Landau levels, we expect a Chern-Simons term in the 
effective action for a. 



The RPA prescription in Fig. 2.4 implies that 



£[a] = ~ 2^2 ^ a ^9i9j) aj + ^ Yl aiai ' ( 2 - 98 ) 
i,j i 

The diamagnetic piece combines with the bare bubble amplitude in equation ( |2.95| ) to give 

2 

m* az+w 2 a i a i- ^ n ^ ne nm h uj — > 0, this diagonal term vanishes. On the other hand, the 
cross term is proportional to tu n in the same limit, leaving 



£[a] = -^-e lj aiUJ n aj , 
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where e 13 is the Levi-Civita symbol. Since the original problem was gauge invariant, the 
complete Lagrangian must contain the scalar potential a$: 

C[a] = -Mie^A a ^ OA ; (2 99) 

where /i, u, A = t, x, y and dt = iuj n . This form is correct to leading order in q, u and shows 
the correct Hall conductivity of the composite fermions. 
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Chapter 3 

Lowest Landau Level II: 
Composite Bosons 



In the previous chapter, we considered composite fermions. In this one, we will consider 
composite bosons. In the former case, we derived an effective theory from first principles. 
However, this is not possible for more than one attached vortex and we follow a phenomeno- 
logical approach instead. The main feature that we find is the magneto-roton excitation, 



which was predicted and analyzed by several authors [55, 58]. Our analysis provides a 
physical picture of this excitation. 



3.1 Formalism 

In this section we generalize the Haldane-Pasquier formalism by considering p objects at- 
tached to the underlying particle. The physical picture that will emerge is similar, with 
the composite particle being a bound state of p vortices and one particle. We specialize 
to v = 1/p so that each vortex carries charge 1/p to maintain neutrality and the magnetic 
lengths are related by 

t?B R =Pl\ h - (3-1) 
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3.1. a Fock Space 

For fractional fillings v the number of vortices in the composite particle is p. The matter 
operators are now p + 1 rank tensors c m;ni ... n , with 

[Cm; m ■■■n p i c n ' 1 --.n^;m']± = ^mm'^nin' 1 ' ' ' ^n p n' p (3-2) 

m runs from 1 to and the n's from 1 to N. The left (physical) density is obtained by 
tracing over all the right indices 

Pmm' = c ni---n p ;m c m;n 1 ---n p (3-3) 
ri\ ■■■n p 

However, there are now p right densities: 

Pnn 1 = c ni---ni_inn i+ i---np;m c m;ni---ni_in'ni + i---np , (3-4) 

m; rij , h a 

where a = 1, . . . ,p and h a implies that there is no sum over n a . Similarly, there are p copies 
of the constraints in eqn. ( |2.32 ), which require | *^^ s mjv ) to be a singlet in each of the p 



right indices. A straightforward extension of the p = 1 case (eqn. ( |2.33| )) shows that the 
basis is 

\mfm-m N \ . . . ,<VV"7jvt J ...J |n\ hfl 

l^phys / _ 6 6 C «i--a P ;mi C /3i--/3 p ;m 2 C 7i-"7 P ! ™iv l U / ' v 0,t -V 

where there are N Greek indices of the type af3 • • • 7. The Levi-Civita symbols and an- 
ticommutation relations of the c's ensure symmetry under the interchange of any pair of 
physical indices m; «-> nij, and we are left with a bosonic Fock space of N particles in 
orbitals. Again, had we started with bosonic c's, we would have ended up with a Fock space 
of composite Fermions. 

3.1.b Physical Operators and Constraints 

The field operator has p right coordinates rji, 



c(z,r] 1 ,...,ri p )= u m(. z ) u n 1 (m)---u^ p {v P ) c m . ni ... np . (3.6) 

m, rii 

It is convenient to change the vortex coordinates to the so-called "center-of-mass" (or Ja- 
cobi) coordinates, which have been used in few-body problems in the context of atomic 
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physics [p9[| . This reduces the p-complex to one center-of-mass coordinate t\ cm = + 
• • • rj p ) and p — 1 relative coordinates £ a (a = 1, • ■ • ,p — 1), which are linear combinations 
of the r/j. We will denote the linear transformation by 



(3.7) 



where a = 0, . . . ,p — 1 and a = stands for cm. Generally, Greek indices will be used for 
the £'s and Latin indices for the 77's. 

One of the nice properties of Jacobi coordinates is that ^/pR is orthogonal: 



|£ cm | 2 + l6l 2 + ---|£ P -il 



^(\vi\ 2 + --- + \v P \ 2 ) 



(3.8) 



For the special where each vortex carries charge 1/p (y = 1/p), we can use this property to 
set the magnetic length of each £ to l\ = £% /p- The utility of this transformation is that 
all coordinates now have only one magnetic length, £b l - 

The particular way in which Jacobi coordinates are constructed is well-illustrated by 
two special cases, p = 2 and p = 3, both of which we will utilize below. For two vortices 
there is only one relative coordinate, so the Jacobi system is 



£1 



2 fa _T ?2) 



(3.9) 



The normalizations are chosen so as to preserve the normalization in eqn. Q3.8[) . Specifically, 
the 1/2 factor in £1 is the reduced "mass", m (with 1/m = 1/mi + 1/771-2) of the two 
vortices — each vortex is taken to have unit "mass" (rrii = 1). 

For p = 3, the coordinates are arranged so that £1 is a vector from rjx to 772 , and £2 
connects 773 to the center of mass of 771 and 772- Each £ is normalized by the square root of 



the reduced mass of the two objects which it connects. Fig. 3.1 illustrates this construction. 
The particular linear combinations are 



Scm 

V ^ 2 J 



( 



1 

71 



1 


1 


1 


\ 


v/3 


v/3 






1 


1 







72 


V2 




1 


1 


pi 


/ 


v/6 


V6 


V 3 



m 
V % J 



(3.10) 
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Figure 3.1: The composite with three vortices attached. The underlying particle is at z with charge 
+1, and the vortices are at ?y; with charge —1/3. The (unnormalized) Jacobi vectors are dotted lines. 
z and the center of mass of the vortices, £ cm , are connected by Ak, which we will later interpret as 
a dipole moment. Since the composite is neutral, it drifts with momentum k. 

There is a certain degree of freedom inherent in assigning the £. For instance, we could 
have chosen £i to connect rji to 773, and £2 to connect 7]% to the center of mass of 771,2 • 
However, each choice obeys the normalization condition in eqn. (j3.§|). Viewed classically, 
|£cm| 2 is a constant in the absence of external forces, so the normalization implies that 
|£i| 2 + • • • + |£p-i| 2 is fixed and that all choices of Jacobi sets can be transformed into each 
other by a member of SO(p — 1). Although this is a fundamental and generally useful 
property of the coordinate system, we will not use it in this thesis, but point it out for 
completeness. The most important property for us is that eqn. ( |3.8| ) scales the £ in such a 
way that all coordinates have the same magnetic length. 

For higher values of p, the Jacobi coordinates can be constructed recursively by grouping 
the vortices into pairs and connecting the centers of mass, and then repeating the process 
with the centers of mass. The details of the general procedure are described elsewhere 
however, here we only need the p = 2, 3 cases. 

Now, the combination of the particle at z and the center of mass of the vortices at £ cm 
is just like a particle and vortex of equal but opposite charge, which can be treated by the 
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noncommutative Fourier transform (Appendix). The transformed field operators become 
c k ({?}) = / ^zfld 2 ^ c(z,r}[---rf p ) M&n>*)Il ■ ( 3 - n ) 

J i=l a=l 

The (5-functions implement the Jacobi coordinate transformation of eqn. ( |3.7|) through 
£a = Rairfi- Since the £ a themselves are coordinates with a magnetic length £g L , the field 
can be cast into the complementary "spin" basis 



<*({£}) =E c W-,<(£i) • ( 3 - 12 ) 

The quantum numbers cr a are nonnegative integers. They are angular momenta of the 
relative coordinates, but we will refer to them simply as "spin". In this chapter we consider 
spinless fermions and bosons so there should be no chance for confusion. 

The left density is integrated over all right coordinates and has the same momentum 
structure as in the p = 1 case: 

^ = / (2<j2 ^ e2lkAq ^ C L|q, CTl ... CTp _ a Ck+iq,^-^! • ( 3 - 13 ) 

The right densities are more complicated, 

p* = (3.14) 

/ 7S2 /n^^ T q(^-^?cm-^)e i<kAqft ^ <_l q ({C})c k+ i q ({e'}) 

Notice the magnetic translation in £ cm — 77^, which is the vector from the center-of-mass to 
the i'th vortex. 

Finally, the Hamiltonian looks just like the unequally charged case, eqn. ( |2~44| ), in the 
previous chapter, 

Again, only "left" operators appear since H is a physical quantity and p^ are constants of 
the motion. 
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3.1.c Remarks on Many-Particle Wavefunctions 

In this section, we augment the physical picture of the composite fermions and bosons by 
outlining the many-body wavefunctions that they can describe. In particular, we recover 
Laughlin's function at v = 1/p 

Consider first N bosons at v = 1/p. There are ./V particle coordinates Z{ and pN vortex 
coordinates rj s i with s = 1, . . . ,p and i = 1, . . . ,N. The constraints (Section 3.1. b| ) require 



that the rj s i dependence of the wavefunction be that of a full Landau level for each s [18], 

*phys(^l,Wl,l'-- ->Vp,U-- -,ZN,Vl,N, ■ ■ ■ ,V P ,n) = /(*Li-- • > *n) \\ {~l s ,i ~ Vs,j) 

s; i<j 

The last factor is the product of p Laughlin-Jastrow factors (the Gaussian factors have 
been left off). Another way of writing it is the product of p Slater determinants of the 
matrices mfj = {uf{v s ,i)}- Each determinant is a Vandermonde determinant and is the 
unique totally antisymmetric wavefuction annihilated by the corresponding constraint. 

The simplest ground state of composite bosons at zero temperature is the single-particle 
condensate 

(ckov^-i) = VP 5 kj0 n^ a ,o • (3-16) 

a 

In coordinate space this can be rewritten suggestively as 

(^o(2,?7i,...,7?p)) = nK z ,Vs) , (3-17) 

s 

where 5 is the delta function with the magnetic length of the vortices, Ib r - This puts the 
vortices on top of the particle, as one would expect from the lowest state due to electrostatic 
attraction. Projection onto the physical basis gives 

N 

■■-,Z N ) = /LK%* e- lv ^ 2/ <n 11 ( %i - rj SJ ) f[(M^l,i, ■ ■ (3-18) 

s,i s;i<j i 



Using the condensate wavefunction ( 3.17] ), / is the product of p identical factors, 



f(zi, ...,Zn) 



i i<j i 



(3.19) 
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Each factor is a Slater determinant of the matrix m%j = {u^(zj)}, which is the same as a 
Vandermonde determinant of the mat] 
exactly the Laughlin state at v = 1/p, 



Vandermonde determinant of the matrix Vij = Zj times an overall factor. The end result is 



f P (zi,. ..,z N ) = H( Zi - ^.) p IK |Zi|2/<i - ( 3 - 2 °) 

i<j i 

Thus, we reinterpret the Laughlin state as a composite boson condensate. 

However, there is a slight surprise if f p is the unique characteristic of a Laughlin state. 
It turns out that the condensate is not required to carry the spin a a = 0; any linear 
combination of the internal states {u^ a (£, a )} will project onto the same f p . In other words, 
the simple condensate in eqn. ( 3.16| ) is generalized to 

<c k (£i, • • • = v^k,o ufo, • • • , II e- |?Q|2/4 ^ , (3.21) 

a 

where u is a suitably normalized anti-analytic function of the £ a . In terms of composite 
bosons, the Laughlin state is infinitely degenerate (for all p / 1). We will come back to this 



question in Section 3.2 when we impose the constraints and model the fluctuations around 
the ground state. 

Despite the infinitely degenerate ground state, the fundamental quasiholes can be rep- 
resented unambiguously. We modify one delta function in equation ( 3.1 7| ) by 5(z,rj s ) — > 



J2m u m+i( z ) u m(Vs), moving the particle at z radially from the origin by one angular mo- 
mentum unit. The effect on f p is an overall factor of z ii which is the form of the quasihole 
wavefunction given by Laughlin j|. This can be seen by writing one of the factors in eqn. 
( 3.19| ) as a Slater determinant; the shifts u^Zi) — > u^ +1 (zi) are equivalent to multiplying 
every column i by Zj and rescaling every row m by an overall numerical factor (due to the 
normalization of u m ), which only changes the determinant by an overall constant. Quasi- 
holes can be moved around by applying magnetic translations f , which are described in the 
Appendix. 



3.2 Energy Functional 

Section p.l.cj explained how the Laughlin wavefunction is a condensate of composite bosons. 
The salient feature was the degeneracy of the ground state. Let us rephrase this in terms 
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of the Hamiltonian ( |3.15| ), which is a function of only the left density, p is 

pL f^y e ¥^l L J . (3 . 22) 



(27r) 2 ^ k-iq,<Ti-a- p _a k+ 5 q,(7i-<Tp_i 

In eqn. ( Eml) of the previous section, we pointed out that the condensate is infinitely 
degenerate because the Laughlin state does not depend on the wavefunction of the internal 
vortex coordinates. In other words, the expectation value of the Hamiltonian is unaffected 
by the internal "spin" state as long as the condensate is {ck >ai ... a - p _ 1 ) cx <5k,o> which is enough 
to give (p^) = ~p. Unlike the composite fermion case in Chapter ||, there is no unique mean 
field ground state to expand around. This degeneracy should disappear once the constraints 
are imposed. 

Here we will follow an alternate route which allows us to include the constraints at the 
outset. We will construct a phenomenological Lagrangian that respects all the symmetries 
and include the constraints as Lagrange multipliers. The resulting Landau-Ginzburg theory 
can then be systematically expanded about a unique mean field solution. 

The simplest rotationally invariant contribution to the Lagrangian consists of a momen- 
tum and a spin piece 

Cm + Cj = -J2 f^|k| 2 + -a) i,...^^...^ , (3.23) 

k,o- D 



\2M 2 J k<J i- CT P-i 



where M, J are constants and a = J2a=l a a- We may guess that the full Lagrangian should 
be 

£ = J! c Lv--a>_i < ^-Ck<7i---<Tj,_i + £m + £j + £ constr , (3.24) 

k,cr a 



where the first term is the usual time derivative [53] and £ C onstr is the Lagrange multiplier 



term that imposes the constraints. However, this expression is not gauge invariant. C must 
be invariant under the symmetries of c, which preserve the commutators of eqn. (|3.2|) , 

v 

c » U L c [J Ur, , (3.25) 

i=l 

where Uj^ is a unitary matrix acting on the i'th right index and Ul acts on the left index of 
the matrix c m;ni ... np . We will find that in order to preserve these symmetries, it is necessary 
to introduce p gauge potentials that act on the right indices. The bare term Cm + £j will 
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acquire the vector components at, \x = x,y, and the Lagrange multipliers, A*, will play the 
role of the scalar potential a l Q . Since the mass and spin terms are decoupled in Cm + £j, 
we will consider their gauge invariant forms separately in the following two subsections. 



3. 2. a Mass Term 

To write the mass term in coordinate space we need the transform of ik, which is the analog 
of V in ordinary space. To this end, define the LLL coordinate operators, Z and Z 

Z = z5(z,£ cm ) and Z = S(z,^ cm )^ cm , (3.26) 

which are adjoints of each other, Z* = Z. In Cartesian coordinates Z = Z x + iZ y and 
Zj 1 = Z^. It is straightforward to verify that the fk are eigenoperators of Z, Z, that is 

[Z * fk] = — ifcfk and Z * f k = ikf\ i , (3.27) 



where the *-commutator is defined in eqn. Q7.40 . As usual, k = k x + ik y and k = k x — ik y , 
so that — e^ u [Z u * fk] = kn fk- Therefore, —it tlv Z v is the analog of the ordinary derivative 
operator when acting on f . To see how Z acts on the matter field, we need to consider 
the Fock space operators in more detail. 

In general, c m;ni) ... jnp can be transformed on the left by some matrix M mm / or on the 
right by some matrix A n / n . that acts on the i'th right index. In coordinate space, M = 
M(z,z') and A = Aj(r^, 7^). We will continue to use the *-operator, but its meaning must 
be clarified when acting on the right coordinates since there are p of them. We define 

(cA.K.^,...,*)-/^^....,* 0.28) 



The index on the operator will always indicate the coordinate on which it acts, rendering this 
notation unambiguous. An example is c* Z, where Z acts on the center of mass coordinate: 

(c* Z)(z,f} 1 ,...,r} p ) = - f d 2 ■q' 1 c{z,rf 1 ,..., rj v )rf^S{rf x , %) + ■■■ 



p J 

+ - f d 2 n' p c(z,r] 1 ,...,rf p )n' p 8(ri' p ,rj p ) (3.29) 
p J 

Let us define the coordinate operator from the right, v l , by 

v % = rjiS (r]i,rfi) and v l = 5{r} i ,rf i )rf i , (3.30) 
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such that 

c*Z=~y^c*v i . (3.31) 
p £ — 1 

The Cartesian components of the u's are Hermitian, vV = vt. Acting on the left, Z * c is 
the same as before, 

= zc(z,rj 1: . . . ,f} p ) . (3.32) 

Now, the Fourier transform, eqn. ( |3.11| ), together with the property in eqn. (|3.27| ) afford a 
definition of the "derivative" operator 

Oc = --[Z* c] and dc=-[Z*c] . (3.33) 

We translate this definition into Cartesian coordinates by using the conventional relations 
d = (d x — id y )/2 and d = (d x + id y )/2, 

d^c = -ie^[Z u *c] , (3.34) 

where e^ v is the Levi-Civita symbol. The non-commutative Fourier transform of d^c is 
—ikuC. Further like ordinary derivatives, obeys the Leibnitz property d{a *b) = da*b + 
a * db, which follows easily from the Jacobi identity [A, [B, C]] + [C, [A, B]\ + [B, [C, A}] = 0. 
The unitary symmetry in eqn. ( |3.25 ) requires that the derivative operators are covariant. 



The inifinitesimal version of the transformation on the right is obtained from the product 
of expanding each = 1 + i A» , 

c i — ^ c + ic * A 

c f ^ c f - iK * c f , (3.35) 

where A is a Hermitian operator, 

A = J2H4,%) (3-36) 

i 

and its conjugate is A = J2i ^-iiVi'Vi)- A standard procedure from field theory can be used 
to make d^c covariant |60[] . Since Z fl *c is already covariant under right transformations, let 
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us consider each term in c* separately, which are of the form c* . First, we introduce 
a gauge potential a^(?7i, %) which transforms according to 

+ . (3.37) 

The derivative of Aj is constructed like d^c, 

d fM A i = -ie» v [vi*A i ]. (3.38) 



Then the combination c* {ef^v^ + aV) transforms according to the rules in eqns. (3.35) and 



c*(e^ + <) ~ (c + ic*A i )*(^vl + al + i^[vl*A i ] + i[a t li *A i ]) 

*2\ 



c * (e"X + a\) + ic * (e^vl + oj.) * A, + 0(Af ) (3.39) 



showing that it is covariant. 



The total contribution of the gauge potentials to e^ v Z v = ^ v v % v is just the sum of 
the individual gauge potentials, which we term a cm , 

< = ;E°i- ( 3 - 4 °) 

" i 

The covariant derivative becomes 

D^c = d^c-ic*a c ™ 
(L> M c)t = ^ c t + ia7*ct. (3.41) 

These derivatives obey [Z)^,^] = iG^, where G^j, = d^a^™ 1 — d v a c ™ — i[a c ™ * a£ m ] is the 



field strength [ SO | . By analogy to standard non-Abelian gauge theory, we will choose the 



transverse gauge for each a 1 

= . (3.42) 



The derivative is again defined by eqn. fl3.38|) . Using the non-commutative Fourier trans- 
form, the gauge condition is exactly q • a q = 0. 

The fully covariant mass term can now be written in coordinate space as 

1 



(iV)t * D„c 

-±-Tr (V + ia™ * ct) * (d,c - ic * af) . (3.43) 



2M 
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At mean field, the saddle point of the full Lagrangian 

<«n = ° 

(Ck^^ap-x) = ^4,011^,0- ( 3 - 44 ) 

a 

After quadratically expanding Cm about the saddle point, we find that the terms linear in 
o^ m are proportional to c^a^ m , which vanishes by the gauge choice. The quadratic term in 
a^ m decouples from the matter fluctuations, and we are left with 

5C M = ~Tr(d^cU d»c + pa™*a™) . (3.45) 

In the next section, we will find a similar decoupling of the matter and gauge fluctuations 
in the spin term. Even if the ordinary kinetic energy k 2 /2M is generalized to an arbitrary 
polynomial in k, the same decoupling would hold at the Gaussian level. Gauge fields will 
be felt, however, as Lagrange multipliers in the constraint terms. 

3.2.b Spin Term 

Turning now to the spin term, we will suppress the left index in this subsection. 

The preceeding subsection defined the operators v i = rjiS^i,!]'^, whose sum was related 
to the center of mass coordinate operator. For the relative coordinates, the natural basis is 

w a = — i RajV^ and vf = R a ,jV^ , (3.46) 
3 j 



where R is the Jacobi transformation £ a = J2j Ra,jVj that was constructed in Section 3.1.b| . 



The i prefactor ensures that, in Cartesian coordinates, is related to e^v^. This also has 
the effect of making the anti-Hermitian, iv?* = —w^, because is Hermitian. 

In the single particle basis u^(^ a ), w a and w a act like i a /V^(-B L and ((.B L /V^)d^ a , 
respectively — see eqn. ( |2.16| ) and the immediatly following discussion. Therefore, the 
quadratic form w a * w a is exactly cr/2, and the spin term can be written as 

Cj = -J^Tr(V* *c f *c*w a ) (3.47) 



■ J H Tr ( 



* d * c* 



53 



Cyclicity of the trace was used to separate w * w in the first line. By changing to Cartesian 
coordinates in the second line, we ignore a constant term in the energy. This contribution 
is from the commutator [w a * ur] = —25 a p to * = (w a * uf" + w a * w a )/2. 

Just as the <9's were not covariant in the last subsection, the w a are not, either. Recall 
from that discussion that each c * e^u* is made covariant by e^ u vl \— * e^ u v l u + a^. We thus 
introduce p — 1 gauge potentials a^, 

a« = Y,Ra,i4 (3.48) 

i 

and the covariant operators Wff, 

c * Wf? = c* w" 1 . + c* af. 
W£i*J = -w°*J + a^*J. (3.49) 

As in the previous section, we work in the Cartesian basis \i = x, y with W x = (W + W)/2 
and W y = (W — W)/2i. The correct spin term is now 

a 

= J^Tr(y%*J -a«*c ] )*(c*w*+c*afj . (3.50) 

a 

The part linear in is equal to 

- JTrc f *c* [w% * a%] . (3.51) 

The commutator is a sum of terms of the form e^ v [v\, * a^J = id^a 1 ^ because [v\, * ajj] = 
unless % = j. According to the transverse gauge choice, all terms in eqn. ( 3.5 H ) vanish. In 
fact we can prove a more general statement for a polynomial spin dispersion. Any term in 
the energy of the form a™ can be written as 

<ct 1 ... ap _ 1 c ai ... ap _ 1 =TrcUc*(W«Uw2r • (3-52) 



The linear piece in always looks like ( 3.51 ) and vanishes in the transverse gauge. For 
our purposes, we will stick with the simplest non-trivial example, Jo. 
The saddle point is 

«> = o 

<Pui-<t,_i) = 11^,0 (3-53) 



a 
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where the left coordinates are suppressed. Quadratic expansion about the saddle point gives 
the perturbation of Cj as 

5Cj = -j^2^(K* cUc *< + p a X) • ( 3 - 54 ) 

a 

The second term is also proportional to J2i * because y/pR is orthogonal. It 

remains to treat the constraints gauge invariantly, the subject of the next section. 



3.3 Constraints 

We will treat the constraints by including p Lagrange multipliers, X l (7]i,r]' i ). The Lagrangian 
then includes the term 

v 

£constr = i^Tr \* * p* ■ (3-55) 

i=l 

The p Ri were constructed in eqn. ( 3.15 ). Gauge invariance under the transformation in 
eqn. ( |3,35| ) requires that 

A* i — > A* + i[Ai * A*] . 

The saddle point for A* is 

=~Kmrt) (3-56) 

or (Aq) = 5qfi in momentum space. Expansion of the fluctuations is analytically feasible 
if the p Ri are expanded in powers of q. It can be done exactly for p = 2, but we will also 
consider p = 3 to finite order. 

At p = 2 there is only one relative coordinate, £ = (r/i — f?2)/2. The saddle point is 



<Cq(0) = V^q,0^(0 

(A q ) = 5q, (3.57) 

The displacement of the vortices from the center of mass was £ C m — ??i = ~C an d £ C m — = £ 
according to the Jacobi mapping in eqn. (|3.9| ). Plugging this into the right densities, eqn. 
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T5| ), and expanding in powers of c about the saddle point we get 

+ ^ | d 2 ecl q (C) e ^ l?|2+ " 9? -^ 1912 , (3.58) 



where c is the fluctuation about (c) and 8p R = p R — p. l\ L is set to one above, and in what 
follows. If we make use of the identity 

e -w-m-\\tf = s(z, -iq) = £ , (3.59) 

CT=0 

then 

cr 

fyq 2 = V^Xl'V"' 7 ^) + C -qa^<r(^) • (3-60) 
a 

This is an explicit expansion in powers of q since u a (q) ~ g CT (aside from an overall Gaussian 
factor). The superscript L has been dropped from u as we are setting £b l = 1. 

Moving on to p = 3, there are two relative coordinates, 6,2 that are given by eqn. ( 3.1C| ). 
The matter field is c q (6>6) with expectation value 



(cq^i,^)) = \/p<5q,oWo(6)^o(6) • 



As in the p = 2 case, we obtain 6m — Vi i n terms of the Jacobi coordinates from eqn. ( 3.1Q ), 



put the result into the definition of p Ri in eqn. ( 3.15 ) and expand about the saddle point 



in powers of c. For example, 6m ~ Vi = ~ \ §6 ~~ Jofei which gives 



^ = (3.61) 

/ d 2 6d 2 6cq(6,6)«o(6)«o(6)e-3l v4€i+v46l 2 -^(v / f« 1 +v / I^)-|kl 2 + ^. c . 

= 2vr| cf^ 2 c^ 1 ^ 2 )S^ 1 ,-iq^)S^ 2 ,-iq^) + h.c. , 
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where we used the delta function identity, eqn. ( 3.59 ). The procedure is analogous for p R2 ' 3 , 
and we obtain in spin space 

^Pq 1 = VpYI 'Vio-^o-i (-iq \/s/2)u a2 (-iq \jl/2) + 

C1CT2 

VpJ2 c -w l a 2 u cA-n\f^)ua 2 (-iq\fy?) ; 

(71(72 

qo-io-2 u o'i (iq \j3/2) u a2 ( iq ) 

(71(72 

J 3 = c q° CT 2 Uct 2 ^) + c -qo<7 2 n ^ (iqV2) . 

(72 (72 

3.4 Density Response and the Magnetoroton 

The results of the last three sections add into the fluctuations of the Lagrangian density, 

p 

5C = Tr (c f * d T c) - 5C M - 6Cj - iJ2 <W_ q <* , (3.62) 

i=i 

where the first term is the usual Berry phase, r being complex time, and the other pieces 
are given by eqns. ( |3.45| ), Q3.54 ), and (|3.55| ). 



We can drop the fluctuations in the vector potentials, a cm and a 1 since they are decoupled 
and will not affect the physical quantities. Then in momentum space, 5Cj + 8 Cm become 

5C M + 5Cj = Y, Ek,a 4(7i-a p _i c k(7i-,7 P _i , (3.63) 

(7of 

where 

^k,(7 = £k + Ja 
k 2 

= 

2M 

J rj = Ja . (3.64) 

As observed previously, both and J a can be arbitrary polynomials in k and u, respectively, 
but the present form suffices for our purposes. 
Let us rewrite the Lagrangian as 

5C = 5<j)\ G^ 1 54>k (3.65) 
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/ 



(3.66) 



kc 

\ c ka J 

The propagator matrix G carries the physical information, including the correlation func- 
tions and the spectrum of excitations via the zeros of its determinant. We consider the 
p = 2, 3 cases in the following sections. 



3. 4. a The case p = 2 

It is a lengthy, but straightforward, calculation to obtain the determinant by Gauss-Jordan 
elimination. It proves convenient to use the basis A = A 1 ± A 2 , and we find, 



Det G" 1 = Peven^odd 



(3.67) 



where 



oo 

n 



U) r . 



odd 



<a=0 
' oo 

n[ 

,<7=0 



E 



\u2a\ Ev2(T 



k(T =0 ^,2(7 



(3.68) 



E 



k,2cr+l 



OO I 1 2 yji 

E P2ct+1| ^k,2o-+l 
T?2 



.(7=0 ^,2(7+1 



and u CT = u a (k). 

The general features of the spectrum emerge already if we keep only the terms up to 
a = 3. We find that the lowest mode is 

1/2 



lUJ r , 



( k 4 £ k ,Q + 8E k ■ 



(3.69) 



which is plotted in fig. 3.2. In keeping with an incompressible fluid, there is a finite gap at 
k = 0. In general, the modes have energy gaps equal to a J with a > 2. The lowest mode is 
a mixture of only cr=even channels, the next one up is a mixture of <r=odd, and so forth, 
alternating in even/odd mixtures. Because the gap would vanish if the constraints were not 
included, or if J — > 0, the low energy physics is dominated by the internal vortex excitations 



of a. The main feature in fig. |3.2| is the dip. According to the approximation in eqn. (3.69) 
it appears whenever J > 1/M. We have plotted the spectrum by including terms up to 
cr = 4, i.e. to 0(q 8 ), and the shape is virtually identical to fig. EO. In proper units, the dip 
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Figure 3.2: Lowest mode at v = 1/2 obtained from the action expanded to 0(k 3 ). k is in units of 
inverse magnetic length. The parameters are m — 3 and J = 1. 

is always located at |k*| ~ 1/£b l - It would seem to correspond to the magnetoroton that has 
been proposed by several authors (5^, |58fl . However, ours is the first analytical observation 
of this phenomena; earlier predictions relied on numerics. Based on the composite fermion 
picture, k is proportional to the dipole moment, or separation between the particle and the 
center of mass of the vortices. It is tempting to conclude that the dip is an electrostatic-like 
feature, occurring when the relative orbit radius of the vortices is equal to the distance from 
the particle to their center of mass, i.e. £/2 ~ £ cm . 

The density-density response function, x LL {l) can be- calculated within the above frame- 
work as well. The connected part of this correlation function is 

x LL { q ) = {5p L q 5p L _ q ) , (3.70) 

where q = (q,itJ n ) and 5p L is the fluctuating part of p L (eqn. Q3.22D ). Expanding about 
the condensate in eqn. ( 3.5 7| ) gives, 



S Pq = vM c q,0«0 + cL qi o«o) • (3-71) 

One way to get the lowest order term in x is to substitute the expansions for p Rl < 2 . A simple 
linear combination from eqn. ( 3.60| ) yields 



Uspf+dp^) = VP E (c^u a (-iq) + cl n(7 u a (-iq)) (3.72) 

cr=0,2,4,... 

= ^Pq+VP E ( C ^U a {-iq) +C t _ q(J U (7 (-^)) . 



<r=2,4,. 
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On the average, the left-hand side is zero because the constraints vanish at mean field, 



(Spf) = 0, so that bp\ ~ C(|q| 2 ) and 



x LL (q) = o(|q| 4 ) 



(3.73) 



just as we had before in the composite fermion case. 

The same result can be obtained by an explicit calculation that extracts the field corre- 
lators from eqn. fl3.66| ). An expansion in minors (Cramer's rule) of the propagator matrix 
G leads to 

(iuj n + e k ) Q + \u \ 2 



where 



( c k,o c k,o) 

( c k,0 C -k,o) 

Q = 



(ujl + el)Q + 2\u \ 2 e k 
KP 

{ujI + 4)<2 + 2 KI 2 £k 

2\u a \ 2 E Ka 

Plugging these correlators into x LL {k) = p((ck,o + c[ k0 ) 2 ), we find 

Ikl 4 



E 



(3.74) 



(3.75) 



o-=2,4,6, 



lim x (k,iw„ = 0) 

k— >0 



16J 



(3.76) 



This result is consistent with an incompressible quantum Hall fluid [55]. Note that the 
limit is independent of M and is valid so long as J ^ 0, which is consistent with the vortex 
excitations dominating the long distance physics. 



3.4.b The case p = 3 



The expansions of the constraints, eqn. ( [3.62 ), is prohibitive in general. To get an idea of 
the spectrum, we cut the expansion off after tri 2 = 2. The spectrum, as obtained from the 
zeros of G is 



ILUr, 



^k,0^k,2 U 4 S k , 2 +4£ k ,o 

771 rp ^ k 4 E kil +4£ k , 2 
^k,l^k,2 ^ 



1/2 



1/2 



(3.77) 



,k 4 E ki2 +4E ki i, 

The second branch in eqn. (3.77) is doubly degenerate. Presumably this degeneracy would 



be lifted in a better approximation. The two modes are shown in fig. 3.3. With the help 
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of Mathematica, we have analyzed the spectrum in more detail and find that the gaps 
are Ja, a > 2. A non-zero gap, of course, is consistent with an incompressible fluid. Once 




Figure 3.3: Lowest two modes at v — 1/3 obtained by expanding the action to 0(k 2 ). The higher 
mode is doubly degenerate, k is in units of inverse magnetic length. The parameters are m = 3 and 
J = 1. 

again, the main feature is the magnetoroton gap. It seems to be a general feature for all p, 
provided that J is large enough relative to 1/M. 
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Chapter 4 

Paired States of Bosons in Two 
Dimensions 



In this chapter we analyze p- and d-wave pairing of bosons from two points of view. First 
as an exact gound state of a particular Hamiltonian in the FQHE, where the bosons are 
really composite bosons. And second, in the framework of BCS theory of paired states of 
ordinary bosons. 

We find that the permanent is an anti-Skyrmion texture sitting on the transition between 
ferromagnetic and helical ordering, both of which are single-particle condensates. The 
transition is of second order so that we can write a continuum quantum ferromagnet action 
pl| for the local magnetization in its vicinity. An analogous description on a lattice is 
straightforward and will be discussed briefly. The results are compared to numerics with 
excellent agreement. 



In Section 4.2 we move on to the Haffnian. As we did for the permanent, we first 
construct a three-body Hamiltonian for which the trial wavefunction is exact, and from 
which numerical diagonalization techniques can extract the spectrum. Unfortunately in 
this case we cannot calculate the spectrum analytically and we resort to an effective BCS- 
type hamiltonian at the outset. A mean field analysis suggests three possible phases in the 
presence of an attractive channel at I = —2: (i) a single particle condensate, or Laughlin 
state, (ii) a pure pair state and (iii) a charge-density-wave phase. We conjecture that 



G2 



the Haffnian lies on the special point separating (i) from (ii). In support, our numerical 
evidence suggests that the Haffnian is compressible and that it may contain incipient long 
lange correlations in the density. 



4.1 The Permanent: p-wave Pairing 
4.1. a Analytic Structure of the Ground State 

Before discussing the actual calculations, we briefly summarize the relevant properties of the 



permanent state. Detailed analysis and related states may be found elsewhere [28, 3C, 62]. 
For the moment we begin by choosing to put the system on a sphere with a magnetic 
monopole in the center |14|| . 

For v = 1/q, the permanent state is a spin-singlet ground state of spin-1/2 fermions for 
q=odd and of spin-1/2 bosons for q=even. The former case is the relevant one here, the 
simplest being q = 1. At this filling factor, the microscopic three-body Hamiltonian is the 
projection onto the manifold of states spanned by the closest approach of three fermions. It 
takes an especially compact form on the sphere, where each particle in the lowest Landau 
level (LLL) has orbital angular momentum N^/2: 

H = V Pijk{\^-1,\) . (4.1) 

i<j<k 

The arguments (3N^/2 — 1, 1/2) are the total orbital angular momentum and z-component 
of spin, respectively, of three fermions at closest approach. is the projection operator 
onto these states for each triplet of particles. Note that the angular momentum is allowed 



to be half integral because the relevant single particles states are monopole harmonics [63]. 
Due to symmetry, three fermions can never have total orbital angular momentum 3N ( p/2, 
so H is properly regarded as a projection onto triplet states of momentum greater than 
or equal to 3N ( f ) /2 — 1. To specify a unique ground state, the total flux, N$, through the 
sphere must be fixed. At = (A — 1) — 1, the permanent is the densest zero-energy 
eigenstate of H. The densest state being the one with the smallest N^. For our purposes, it 
is most convenient to project stereographically onto the plane, where the i'th particle has 
the complex coordinate Zi = X\ + iyi. Denoting the spin component of the i'th particle by 
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°i = T«i iii casts the permanent into the form: 



t i \ - i r Tp(2i-l)lp(2i) ~~ lp(2i-l) Tp(2/ i yr > > 

^permOlOl, • • • , Z N/2 a N/2 ) = V TT — ^ ^ ^ ^ ^ ^ JJ>* - (4-2) 

P i=1 «p(a»-i) - z P(2i) ^ 



The prefactor is the permanent of an x N matrix, which is a determinant with the sign 
of the permutation omitted, and P stands for all permutations of ./V objects. The second 
factor is the usual Laughlin-Jastrow ansatz with the Gaussian factors omitted. The BCS 
pairing structure of the prefactor is manifest in the above form; the spatial part of the pair 
wavefunction is in an orbital angular momentum / = — 1 eigenstate and the spin part is a 
singlet. Therefore the prefactor is totally symmetric, representing a paired wavefunction of 
bosonic coordinates, which are non other than the composite bosons. 

The factors in the denominator ensure that the relative angular momentum of each 
disjoint pair (zi — Zj) is reduced such that each projection by P^ gives zero. In other 
words, H penalizes those states which do not appear in the wavefunction. Addition of one 
flux quantum preserves this property and there is a space of zero energy states including 
the Laughlin state. Equivalently, reducing the flux by one quantum through the Laughlin 
state creates an anti- Skyrmion, which is a uniform spin configuration on the sphere (a 
"hedgehog") costing zero energy. Stated in yet another way, extra flux in the permanent 
creates quasiholes (or edge states), which belong to a degenerate manifold of zero-energy 
states |3(], [34|] . The advantage of the q = 1 Laughlin state is that it is a Slater determinant 
of single particle states, which in this case is a Vandermonde determinant: 

(A 

i y o j . i<j 

where (1 0)j is the spin state of the i'th particle, so that the total z - component of spin 
is S z = N/2. The spin-wave states correspond to superpositions of the degenerate states 
defined by S z = N/2—1. Notice that the prefactor of the Laughlin-Jastrow factor is trivially 
constant and symmetric under particle interchange. By analogy to the permanent (4.2), we 
interpret the prefactor as a wavefunction of composite bosons. 

In the following subsection we switch to the plane and construct the spin wave excitations 
in the LLL using magnetic translation operators, with which we obtain the exact spin wave 
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spectrum at all wavevectors. We do not utilize the composite boson picture yet and all 
results in the next section are completely exact. 



4.1.b Magnetic Translations on the Plane 

A single magnon mode with momentum k is represented in terms of spin density operators 
by = J2jLi e~ lk ' rj a~ , which periodically flips one spin of a completely polarized state. 
However, all operators must not have any components in the higher Landau levels. To 
impose this constraint, it proves necessary to project all operators into the LLL using the 



Bargmann-Fock representation |69], 65]. As a consequence, a magnon is composed of flipped 
spins along the direction perpendicular to k. 

The fundamental operator that we will need is the projected one-body density operator 



Pk( r i) = e _ik ' ri . Here and throughout this chapter the overbar denotes projection into the 
LLL. This operator is discussed in detail in the Appendix. Because translations no longer 
commute in the LLL due to broken time translation symmetry in the presence of a magnetic 
field, a charged particle picks up an Aharonov-Bohm phase as it traverses a closed circuit. 
However, there is a one-to-one correspondence between the ordinary translations in zero field 
and the magnetic translations. In zero field, /3k( r i) translates the i'th particle by a distance 
r, and the corresponding magnetic translation, Tk(i), is a translation by l B z x k, where Ib is 
the magnetic length and z is the direction of the magnetic field. The cross-product appears 
as a consequence of the Lorentz force. The explicit relation is 

TkW = e^-pkfa) , (4.4) 

where the magnetic length has been set to unity — as it will be throughout the chapter. The 
ubiquitous Gaussian factor can be thought of as the momentum space version of the most 
localized wavepacket in the LLL, or the projected delta function. The total spin and charge 



density operators are easily obtained from (4.4): 



_ N -4* 

i=l 

N & 

?k = E e ""^K, (4.6) 
i=l 



65 



where a a is the a'th Pauli matrix. The projected spin density flips one spin and translates 
the resulting state by z x k, which is a quasihole-quasiparticle pair. The limit of large 
k corresponds to a large pair so we expect that the spin wave spectrum approaches an 
asymptotic value in this limit. 

The Laughlin state ( |4.3| ) is a Slater determinant, which allows the density-density cor- 
relation to be determined analytically Q at v = 1. In terms of the magnetic translations 
in ( [4.5D , the pair correlation function (projected onto the LLL) is: 

2 2 2 

(*L| - 7W^ l*L> =5>-V e -V(T q (i)T_ k (j)) = -(P> 2 W^> 

or 

£<r q (i)r- k (i)> = -<P> 2 5 kq , (4.7) 

where the colons remove self-correlations by normal ordering. Equation (fO|) is the building 
block for much of the subsequent analytical results. 

The non-commutativity of magnetic translations is irrelevant in ( |4.7| ) because the prod- 
ucts involve different particles. However, the translations do not commute for the same 
particle. Rather, they form a representation of the magnetic translation group with the 
algebra: 

q A k 

[Tq(i), Tk(j)] = 2i(%r q+k 00 sin — ^— (4.8) 

r q (i)r k (i) = r q+k (i)e^ Ak , (4.9) 

where the wedge product stands for (q x k) • z. 

Using the projected spin density operator we can easily construct a single magnon by 
applying the spin density operator ( |4.6| ) to the Laughlin state: ^ k = 5 k ^L. 



4.1.c Spin- Waves and Instability; Exact Results 



The translationally invariant version of the three-body Hamiltonian (4.1) on the plane can 
be written as: 

H = v E - r Mn ~ r o) ( 4 - 10 ) 

Of course, all observables are to be calculated after projection into the LLL. Since the 
Hamiltonian conserves the total spin, as well as being translationally invariant, the state 
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with one magnon is an exact eigenstate of H with energy given by the usual expression p6| : 



{* L \SZ H, 5 k 



(4.11) 



{* L \St S k |* L > 

Within the expectation value of H, the two gradients reduce the orbital angular momentum 
of a triplet in and in ^ k by one, hence the identity with the spherical Hamiltonian Q4.ip . 

H can be rewritten in terms of the real-space density operator p(r) = 5{r — rj) in a 
similar manner: 



H = V J d 2 xd 2 xW : V 2 p(x) p(x') p(x") : <J(x - x'^(x' - x") 



(4.12) 



qp 



The last line is obtained by a Fourier transform and by the relation between the density 



and magnetic operators (4.5). Substitution of the last line into the energy expression ( 4.11 ) 
and use of the magnetic operator algebra Q4.8|) and ([Op reduces the dispersion to: 



u; k = —V 2_, Q~ r 1 ' 
q>p 



2 e 2 x 



(4.13) 



E (-r-^K- P U)r p (k)) L ( e J P Ak + e^-P) Ak + e^ Ak - 3) . 

Rather than calculating the triple-density correlator directly, observe that everything inside 
the expectation value is already in the LLL. Thus, the magnetic translations can be replaced 
by density operators according to the rule in ( |4.5| ). After Fourier transforming back into real 
space, we are left with terms like (:V 2 ./9(x) / o(x / )/9(x"):)(5(x — x / )<5(x / -x" + kx z), which can 
be evaluated by writing the density operator in second quantization, p(x) = (^(x)?/;(x)), 
and using Wick's theorem. For the spin-polarized v=\ state, the Green's function is known 
exactly [|[, being, in the symmetric gauge, 

($(xl)$(x!)) l = (p) e -il x - x 'l 2 e f xAx ' . (4.14) 

The end-result is the exact spin wave spectrum: 



w k = C 



1 - e" fe2 /2 (* + i 



(4.15) 



where Q is the area of the system, (p) = 1/2tt1 2 b , and the asymptotic value is C = lQV(p)Q. 
Equation ( 4.1 5| ) is the central result of this section. It should be borne in mind that, 
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although it was obtained for the Laughlin state ( [4.3[ ), it is also valid for the permanent 
(|4.2| ), since the two are degenerate in the presence of our three-body interaction (4.1). The 
comparison to exact numerics is shown in Fig. 4.1 , 



0.2 

0.15 
=k 0.1 
.05 






Figure 4.1: Spin- Wave Dispersion: Points are numerical data for N = 10, iV^ = 9 (computed by 



E. Rezayi). Solid line is the spectrum in (1.15), with C determined by matching to the data at the 



highest k. For comparison, the dotted line is the spectrum for a two body interaction ( 4.1 6| ) with 
V 2 = 0. 

At large k, uj^ approaches an asymptotic value as expected. This is the energy of a 
widely separated quasihole-quasiparticle pair, as in earlier work on two-body interactions 
floTfl , However, the novel feature of ( 4.15| ) emerges at small k, where uj^ ~ k 4 , in contrast to 
the usual quadratic dependence. In other words, the spin stiffness is exactly zero when the 
interaction is the three-body Hamiltonian fl4.10| ). This is precisely what was conjectured 
earlier based on numerical evidence |3(|. Here, however, we have an exact calculation 
verifying that claim, and the comparison with the original data is quite good, as shown 



in Fig. 4.1. Therefore, the permament is poised on the brink of an instability because lu^ 



becomes negative as soon as the spin stiffness dips below zero. 

The degeneracy of the polarized Laughlin state and the unpolarized permanent can also 
be understood in light of the long wavelength behavior of the spectrum. At zero spin 
stiffness, it should cost no energy to create a slowly varying spin texture of unbounded 
size. In fact numerical calculations of the spin-spin correlator {S z (x)S z (y)) on the sphere 
confirm that the spins are approximately anti-aligned at antipodal points, indicating an 
anti-Skyrmion texture with long range order (Fig. W^). In the next subsection we will see 
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that pairing of composite bosons predicts precisely this ordering. Skyrmions are not new in 
the FQHE |38], and are typically associated with excitations when extra flux is added. 
Here, however, they appear at one fewer flux quantum. 




-0.3-1 — . — , — . — , — . — , — . — , — . — , — . — , — . — , 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



Figure 4.2: (S Z S Z ) of the Permanent and, for comparison, of the Laughlin state with one fewer 
flux. A is the magnetic length. (Computed by E. Rezayi) 

The three-body spectrum obtained above cannot exhibit an instability because it is 
always stable — i.e. projection operators generally have no negative eigenvalues. However, a 
negative quadratic term can be restored by including a short-range, two-body interaction. 
For concreteness, we model such a potential by the simplest, non-trivial expansion, 

H 2 = E ^(xi - xj) + ^ 2 Vf(5(xi - x,-) , (4.16) 

in place of H. There is no restriction on the sign of V2, but Vo should be positive. Except 
for the special case V2 = 0, for which H2 is thepseudopotential Q at v = 1, the Laughlin 
state is no longer an eigenstate. However, due to translational invariance, the state with 
one magnon is an eigenstate, so that its spectrum may be calculated exactly once again. 
Projecting H2 onto the LLL and following the same procedure as before yields: 

w k = -^EV^ 2/2 (e iqAk + e^ qAk -2) , (4.17) 

2 q 



where V q = V + \V 2 q 2 . Eqn. ( pTp agrees with earlier calculations using many-body 
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techniques p7| . Unlike the three-body case, the coefficient of the quadratic term in the 
dispersion does not vanish in general. However, by tuning the interaction parameters, the 
stiffness can be forced to zero, and the two- and three-body interactions will behave similarly 
at large distance. For example, the ratio V^/Vq = —2/3 mimics the three-body spectrum 
(|4.15| ) at small k. 



When V2/V0 < —2/3, uj^ looks like the familiar "Mexican Hat" potential, with a 
minimum at the wavevector Q 2 = IV0/V2I. At zero temperature, all the magnons (or 
quasiparticle-quasihole pairs) condense into the momentum Q. Nonetheless, the total ve- 
locity of the superfluid is zero because the up-spins condense into Q/2 and the down-spins 
into —Q/2. We illustrate this explicitly by a Hartree-Fock function in the LLL with the 
variational parameter Q. At this point this is more of an ansatz than anything else, but we 
will justify this in the next subsection. For the moment, consider the single particle orbitals 

Mr) = ± { efQ n r ) ^' 2|2/2 - ( l ^ Z? ' 2 > ( 4 - 18 ) 

^ { e-^ r ) { p_ Q/2 (r) ) 

where PQ/2(r) is the one-body density operator that is to be projected into the LLL accord- 
ing to the rule in ( |4.4| ). The Slater determinant, ^q, of these orbitals exhibits helical spin 
order and reduces to at Q = 0. In terms of the magnetic operators, \I>q is: 



^ 1 ,..,^)=nV° 2/4 TQ/2W n^-^)e-^^ 4 . (4.19) 




Let us rewrite as the spinor operator times the purely spatial part, V^q = Tq^q. The 
energy Eq of can be calculated in exactly the same way as the spin-wave energy (by 
using ( fLlip ). We find that Eq is identical to loq, i.e. 

^7+ 



(Sq 



H2, Sq 



(Sq S 



Q 



q J q/i 
H2, Tq 



Eq = " _ ' ■ (4.20) 

<7q Tq) q 

In fact, this identity holds generally for any translationally invariant interaction that can 
be expressed as a product of charge density operators. In particular, it is true for both of 
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the three- and two-body interactions, illustrating that our picture of a magnon condensate 
is internally consistent. 

The wavevector Q plays the role of an order parameter which increases continuously 
from zero as one crosses over from the ferromagnetic into the helical phase. On the other 
hand, the magnetization undergoes a first order transition. The boundary between the two 
phases is defined by p s = 0. In the following subsection, we prove that helical order is 
incipient in the ground states of the permanent by adopting the composite boson point of 
view. A similar procedure on the sphere will show that the analog of the helical ordering is 



exactly the anti-Skyrmion, as was suggested by the numerical data in Fig. 4.2 



4.1.d Spin Order of the Permanent 



Let us begin by recalling the general theory of paired bosons in zero magnetic field [23]. An 
effective Hamiltonian which captures this physics is of the BCS type: 



Keff = E 



£kcL c k<x + ^ (Ak c -k|C kT + A k 4 T cl ki 



(4.21) 

where £ k = £k — M and £k is the single-particle kinetic energy and A k is the gap function. 
In the fractional quantum Hall effect, the quasiparticles entering K e ff are the composite 
bosons, c kcr , which see no magnetic field. We assume that e k ~ k 2 /2m* at small k, where 
m* is an effective mass of the composite (see for example Chapter |2[). For p-wave pairing, 
we take A k to be an eigenfunction of rotations in k of eigenvalue I = — 1. At small k the 
generic form of the gap function is thus 

A k ~ A{k x - ik y ) (4.22) 

where A is a constant. Although there is no explicit single particle condensate in K e ft, it 
will be emerge naturally below as being equivalent to pure p-wave order. 

More rigorously, one must solve the self-consistent gap equation when the interaction 
contains an attractive I = — 1 channel. Consider a non-singular interaction which has a 
power series expansion at short distance: V^(k) = ao + «2^ 2 + • • • • The gap equation is 

A k = -^E^( k - k ')^- (4-23) 

k' k ' 
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where the quasiparticle energy is 



£ k = ,/^-|A k | 2 (4.24) 



Note the minus sign in contrast with the familiar fermion case. This is a complicated non- 
local integral equation, but it separates if we assume that the gap function contains only 
one angular momentum channel: A k = |Ak|e 4 ^, where (j) is the polar angle of k. The 
interaction expands similarly into angular momenta: 

oo 

V{k-k') = V m {k,k')e- im ^-^ (4.25) 

m=— oo 

where the coefficients Vi(k, k 1 ) depend only on the magnitudes of k and k'. It is straight- 
forward to show that the leading order behavior of V; in k is 

Vi(k,k')~k l (4.26) 



Substituting the expansion ( 4.25j ) into the gap equation ( 4.23 ) yields precisely the p-wave 



gap ( |4.22| ) at leading order in k whenever V—i(k, k') is negative and Vi = for I ^ —1. For 
general /-wave pairing the gap is proportional to {k x — ik y ) at long distance. 

The quasiparticle energy E^ contains important physical information. When \i < and 
A is small enough, the gap is Eq = fj, and 2Eq is the energy needed to break a condensed 
pair. On the other hand, as the gap closes the paired state ceases to exist and a single 
particle condensate appears |23fl . Precisely at n = the spectrum is unstable for any finite 
A. However, as we will show shortly, the pure pair state is really a single particle condensate 
and the solution of K e ff is fully consistent when we expand about this new minimum. 

In the absence of any single particle condensates, a pure pair state of spin-1/2 bosons is 



lQ) = ^ exp {^ 5kC kT c -k| | 10) (4-27) 



where M is the normalization given by 

To ensure that |0) is normalizable it is necessary that |(/k| 2 < 1 f° r an k and p-wave order 
requires that is antisymmetric in momentum space: = —g-k- Therefore, due to 
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bosonic statistics, the lowest allowed spin state of a pair is the spin singlet — so 
\Q) has total spin 5 = 0. 

In real space, the (unnormalized) component of the pair wavefunction with N particles 
(N even) is 

N/2 

*(riai, . . . , r N a N ) = J2Ii # ( r P(2i-i) - r P(2i)) (tp(2i-i)lp(2i) - ip(2i-i)tp(2i)) ( 4 - 29 ) 

P i=l 

where g(r) is the inverse Fourier transform of and P runs over all permuations of N 
objects. This is the form of a permanent of an N x N matrix, which looks like a determinant, 
but with the sign of P omitted. In fact, it is the analog of the Pfafnan (which is the 
determinant of an antisymmetric matrix) for paired fermions. 

If we are dealing with the FQHE, then the bosonic operators Ck<r really originated as 
composite bosons. During projection to the LLL, one typically picks up a cutoff factor on 
(7k of exp(— /^|k| 2 /2), where Ib is the magnetic length. We will neglect this factor in all 
that follows with the understanding that K e ff and <?k are valid at long distance. With this 
caveat, it is easy to see that if 

with A a constant, then the asymptotic behavior at long distance of the inverse Fourier 
transform is 

g(r) ~ - . (4.31) 

z 

Comparing this with fl4.29| ) , we find that at long distance, the pair wavefunction is precisely 
the permanent prefactor of the LLL state defined by ^perm (|4.2| ). 

While the Laughlin-Jastrow factor is taken care of by the transformation to composite 
bosons, it should be borne in mind that the price of this projection is an extra constraint or 
a fluctuating (Chern-Simons) gauge field. We are neglecting these effects at the mean field 



level. Analogous questions for pairing of composite fermions have been raised recently [72]. 

Now let us consider the occupation number at wavevector k: (nk) = (c^Ckf) + (c k j,Ckj,). 
From the form of \Q) ( [4.15 ) this can be written as a function of g^ only, 

<n k > = T^T|2- ( 432 ) 

1 - m\ 
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Substituting the asymptotic behavior (|4.3C| ), we can write the total number of particles as 

The occupation numbers fall off algebraically at small k and even more quickly (exponen- 
tially) on the scale of k > 1/£b- The condition of normalisability of the ground state 
guarantees that |A/fe| 2 < 1 for all k — in other words, |A| must be less than the minimum 
wavevector, |k m j n |. To make sense of this expression we now impose boundary conditions 
that compactify the plane into a torus. For simplicity, consider an L x L torus in the 
xy-plane; generalization to L x x L y and modular parameter r is straightfoward. 

On the torus there are four degenerate ground states for the permanent Hamiltonian 



(4.2) corresponding to periodic (++) or antiperiodic, (H — ), (— +), ( ), boundary condi- 
tions in each of the two directions. When both directions are periodic (++) the minimum 
reciprocal lattice vector allowed is k m j n = and iV becomes sharply peaked (in fact it is a 
delta function) at hq. One of the spin directions (say, j) is singled out by the correlations 
and we are left with nothing other than a Bose condensate with 

M = y/NSkflS^ (4.34) 

which is the spin-polarized Laughlin state. On the other hand, in the antiperiodic sector, (-1 — ) 
or ( — h), |k m j n | = tt/L. Inverting (|4.33| ), we find that |A| = |k m j„| — 0(1/N) and in the 
thermodynamic limit one of these two sectors is macroscopically occupied. For example, in 
the (-1 — ) direction, k m j n = (tt/L, 0). Since the total momentum of the condensate must be 
zero, we occupy (+7r/L,0) and (— tv/L,0) with equal probability: 

[n [n 

(ci<t> = y y^k,k mm and (c k ;) = W y^,-k mm (4.35) 

In real space this is precisely the helical winding with Q = 2k m j n which corresponds to 
the spins winding exactly once over the length of the system L in the x direction. For 

the remaining antiperiodic sector, ( ), |k m j n | = \/2it/L and the winding is along the 

diagonal. 

In any case, we have shown that three of the ground states of the permanent Hamiltonian 
are really single particle condensates with helical spin order and the remaining one is the 
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Laughlin spin polarized state. As we tune through the transition, one or the other long 
range order takes over. This justifies our Hartree-Fock ansatz ( 4.20| ). In the following 



subsection (4.1.e| we shall expand about the single particle minimum by using an effective 
Landau-Ginzburg theory instead of K e ff- 

Finally, we repeat the analog of the above on a sphere, which provides a nice intuitive 
picture of the anti-Skyrmion texture. Recall that the permanent state on the sphere has one 
fewer flux quantum than the Laughlin state. Thus, composite bosons live on the surface 
of a sphere with a magnetic monopole of strength one. The appropriate single particle 
states are monopole harmonics, Yl,a/(#> 0)) with angular momentum L > 1/2 and M is the 
magnetic quantum number in the range — L < M < L [63]. A pair wavefunction must be 



rotationally invariant just as it is translationally invariant on the torus or the plane. The 
unique bilinear scalar is 

oo L 

E E Y LM {e,(j>)Y L ^ M (d'A')f L {W\L,M-,L,-M), (4.36) 

L=-l/2M=-L 

where (0, 0|L, M; L, —M) is the Clebsch-Gordan coefficient for coupling \L, M) and \L, —M) 
into the orbital singlet |0, 0) and /l is a function of L only. This particular Clebsch-Gordan 
coefficient is equal to a function of L times (— 1) , i.e. the only M-dependence is in 
the phase factor (— 1) [74]. If we combine all of the L dependence of /l and the angular 
momentum coupling into a single pair amplitude, g L (—l) M , then we can write the many- 
body paired state as 

-1) M 4,mt c I-M|||0)- (4-37) 

Note that M is half-integral so the pairing amplitude gL(—l) M is antisymmetric under 
M — > — M which is consistent with bosonic statistics and with the planar symmetry k — > 
— k. The relationship of angular momentum to linear momentum is L = |k|i?, where R is 
the radius of the sphere (Haldane in ref. |l[). Therefore, by analogy to ( 4.30| ), we expect 




that g ~ 1/L at small L, although we have no explicit proof of this statement. Fortunately, 
the exact form of g L is not important for the following. 

The lowest Landau level of the composite bosons has only two available states: [1/2, ±1/2) 
By analogy to the torus, the spins ought to condense into these two lowest states with equal 
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probability. In real space this says that 

< Ci (M)) = ^1/2,-1/3^,0) = -e*Vl-cos0 

, (4.38) 

(c T (0, 0)) = Y 1/2tl/2 (9, 4>) = VI + cos 

The expression for the first few monopole harmonics may be found in ref. p3| ] , and a factor 
of y/N/2 has been omitted from the right hand sides of ( 4.38| ), With the above condensate 



it is easy to calculate the expectation of the spin density (Si) = (c^af^cp), where is 
the i'th Pauli matrix with spin indices a and (5: 

(S z ) = — cos 9 

(S x ) = -sm6 cos 4> (4.39) 
(S y ) = — sin 9 sin (j> 



This is precisely an anti-Skyrmion spin ordering as the numerical data in Fig. |4.2| shows. 

In summary, the mean field theory of composite bosons is completely consistent with the 
analytical and numerical results of the previous subsection. Furthermore, it seems that pure 
p-wave pairing of bosons in two dimensions should be viewed as a single particle condensate. 

4.1.e Effective Field Theory Near the Transition 

An effective continuum quantum ferromagnetic action (CQFM) for the polarized FQHE 
has been proposed recently by Read and Sachdev |6l]]. The idea is to write a sigma model 
for the local magnetization |6q1 , n. Their model includes terms up to momentum squared, 
the higher order terms being irrelevant in the renormalization group sense. We modify this 
CQFM by stabilizing it in the helical region (p s < 0) by adding terms that are quartic in 
momentum. At the end of this subsection, we will briefly discuss this theory on a lattice. 
The CQFM Lagrangian density is 

Co [n] = U(n) • d T n + ^ (Vn) 2 , (4.40) 

where r is complex time and the time derivative term is the Berry phase. A is the monopole 
vector potential such that V„ x A = n. Although the model defined by Co is unstable 
when the stiffness is negative, higher order derivative terms can stabilize this region. The 
symmetry broken helical phase contains an SO(2) x SU(2) residual symmetry; £0(2) for 
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rotations of Q in the plane and SU(2) for the spins. There are three terms at leading 
non-trivial order that obey this requirement: 

where a, b = 1, 2 and i, j = 1, 2, 3 and repeated indices are summed over. A renor- 
malization group analysis shows that mode elimination in Cq generates a combination of 
only the first two terms: C,j = J(2d a nidbnid a rijdbnj — d a riid a nidbnjdbnj) , which is associ- 
ated with spin wave scattering fl6l| . The third term, Ck = K (V 2 n) 2 , is associated with a 
second-nearest neighbor interaction on a lattice. If only these two terms are retained, then 
the total CQFHM Lagrangian is given by 

C[n] = C [h}+Cj[h} + C K [h] . 

Although Cj cannot introduce any k 4 terms into the dispersion since it arises from mode 
elimination, it will emerge that Ck is sufficient to reproduce all of the long wavelength 
features found in the previous subsection. We choose to keep Cj at this point for added 
generality. 

Using spherical angles, n can be parameterized by cj), its direction in the plane and 9, the 
fluctuation out of the plane. Small deviations ($, ip) from helical ordering with wavevector 
Q are given by 9 = cp = Q • x + ip. The fluctuations in the spherical angles should obey 
|Vy?|, |Vi?| <C |Q|. In these coordinates, n = (cos ~d cos(Q-x+( / 9), cos $sin(Q-x+9j), — sin ■#), 
and the ferromagnetic phase is recovered when Q = 0. The Berry phase reduces to the 
simple expression, i{)d T ip. 

The mean field energy density of the helical state is given by 

MQ 2 ) = \ [psQ 2 + {J + K)Q A ] , (4.41) 

which has the desired shape when the stiffness is negative. The spectrum of C can be found 
by including fluctuations to second order in ip and i?. The Green's functions (ip{— k, — u>)ip(k, u)) 
and {'&(— k, — w)-#(k, u)) both have poles at 

iw k = (4.42) 



< 



\ { [4J(Q • k) 2 + K(k 2 - Q 2 ) 2 ] [4(J + K)(Q- k) 2 + Kk 4 ] } 1/2 p s < 

(i/a 2 + \Kk 4 ) Ps > o, 
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where Q is the momentum which minimizes Eq. For certain ranges of the parameters (J, K), 
the spectrum is always positive, and C describes a stable system. For simplicity, consider 
the helical state with J = and K positive. Then the spin wave energy of the CQFHM 
simplifies to 

iu k = y | k 2 - Q 2 1 ^4(Q-k) 2 + fc 4 (4.43) 

As required, ~ k A when p s zero. The three-body spectrum obtained in the previous 
subsection can be reproduced by identifying K/2 with the coefficient of the k A term in 
( 4. 43] ). On the other hand, when p s is negative, we obtain the reasonable behavior, at small 



k: 

[ k 2 if k_LQ 

iuj^ ~ < (4.44) 

I k if k || Q 

This may be expected because the spins are alligned ferromagnetically perpendicular to Q 
but anti- ferromagnetically parallel to Q. 

4.1.f Lattice Model 

Before leaving the spin waves and moving on to the charged excitations, we briefly sum- 
marize a lattice model of spin-1/2 bosons that exhibits the helical transition. We represent 
the spin sector by two types of hard-core bosons hopping on a lattice with a short-range, 
spin-dependent interaction. 

Defining bi a to be the bosonic destruction operator on site i, consider the following 
Hamiltonian: 



Hi a t = ~t b L b j° ~ vYl Hi + U Y1 ni ( ni ~ !) + J Y1 ( Si ' S i ~ \ n i n i ) ' ( 4 - 45 ) 

<ii>,CT i i (ij) 



where is the total number of bosons per site, and {ij) denotes the sum over nearest 
neighbors. The [/-term is required since we are dealing with hard-core bosons; the singular 
gauge transformation which mapped the fermions to bosons maintains the repulsion at 
the same site. The final, antiferromagnetic, term is very much like the fermion t — J 
interaction |66| ]. In the continuum limit, it reduces to the rotationally invariant interaction, 
|e°" r '0 o -VV' T | 2 . We will not discuss Hi at further, save to point out its mean- field features. 
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This two-component model is similar to the one-component lattice boson model con- 



by large t/U or special values of [i where H[ a t has particle-hole symmetry, one is free to 
consider bose condensation of the particle fields. If one also makes the self-consistent re- 
striction t/U < Jn/2U, n being the average number of bosons per site, then the free energy 
of the helical phase is lower than that of the ferromagnetic phase. As in the continuum, 
this winding is described by the condensation: 



where iVj = = N/2 and Nl is the number of lattice sites. The optimal condensate 
wavevector lies along the diagonal of the lattice: Q = Qx + Qy with Q determined by the 
solution of cosQ = 2t/Jn (in units of the inverse lattice constant). 

As the hopping decreases, provided that one is in a given region of fJ,/U, the ground 
state crosses over into a Mott insulator and it is no longer valid to argue based on Bose 
condensation. However, at points of particle-hole symmetry, the superfluid persists down 
to infinitesimal hopping. For the Hall liquid, the insulating phase is most relevant since the 
charge excitations must be gapped. It would be interesting to map out in detail the phase 
diagram of this magnetic superfluid to insulator transition. 

4.2 The Haffnian: d-wave Pairing 

The previous section has analyzed the spin sector of the permanent in some detail. Its spin- 
wave dispersion was shown to be soft, allowing a magnetic transition. We now want to ask 
the question whether there is a ground state wavefunction whose density sector exhibits an 
analogous behavior. To this end, we introduce a d-wave paired wavefunction, or "Haffnian", 
describing hard-core, spinless bosons at filling factor v = 1/2. In the FQHE, it is a d-wave 
paired state of composite bosons. We will then argue that the Haffnian is compressible 
and sits on the phase boundary between incompressibility and non-uniform charge density 
order. At present we have no effective field theory that captures this behavior. 



sidered earlier by M. Fisher et al. [f75ll . In the superfluid regime, which is characterized 




(4.46) 
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Although d-wave paired states have been proposed in a different contex by Wen and 
Wu their Haffnian wavefunctions describe incompressible states of fermions. As our 
proposal is somewhat different, in the following subsection we present the relevant con- 
structions in some detail, mainly along the lines used to investigate non-abelian statistics 



[30]. 



4. 2. a Analytic Structure of the Haffnian 

The Haffnian builds in pairing into the Laughlin state, much like the Pfafnan of Moore and 



Read f2§: 



= V2 1 IIC* - *i) 2 ■ (4-47) 

p (*P(1) -ZP(2)) 2 ■■■(ZP(N-I) -Z P{N) ) 2 f^ 

This describes an even number, A, of spinless bosons at half-filling and flux A^ = 2(N — 1) — 2 
(on the sphere). The prefactor is known as a Haffnian in linear algebra and is also the per- 
manent of the A x A (A > 2) matrix M {j = l/fe - zj) 2 (i / j, Ma = 0). The Haffnian, 
Pfaffian and determinant are related by several identities, which may be found in, for in- 
stance, Greiter et al. |3l|] 

To construct the parent Hamiltonian for ^?Hf, it is convenient to work on the sphere, 
where each particle has orbital angular momentum N^/2. Using the same notation as in 



the permanent (4.1), Unf is a sum of three-body projection operators: 



H Hf= £ V Pi jk (3N^/2) + V 2 P ijk (3N lf> /2 - 2) + ^^(3^/2 - 3) . (4.48) 

^ Hf is the unique zero-energy eigenstate of Huj at A^ flux. The proof of this statement 
proceeds by showing that the maximum angular momentum of any triplet in ^ Hf is 3N l / ) /2— 



4; the details are in Section 4.2.b. Note that projection onto angular momentum 3A^/2 — 1 is 
absent. This is a consequence of the symmetries of Clebsch-Gordan coupling; three spinless 
bosons of angular momentum L cannot be in a total angular momentum state of 3L — 1. 

Alternatively, ^ Hf can be rewritten explicitly as a paired state of composite bosons. By 
analogy with the permanent, the order parameter (on the plane or torus) is an eigenstate 
of angular momentum with eigenvalue I = —2, i.e. A ~ A(A: X — ik y ) 2 to leading order in k. 
Likewise the many-body state is a BCS wavefunction of d-wave bosons in two dimensions 
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and the asymptotic behavior of the pair state reproduces the 1/z 2 nature of the Haffnian 
prefactor ( [4.471) . 



4.2.b The Haffnian Hamiltonian on the Sphere 

To see that ^Hf does not contain any triplets of total angular momentum greater that 
3N ( p/2 — 4, we generalize Haldane's original argument for two-body interactions p4|1 . This 
argument is easily applicable to n-body interactions. The factors (zj — Zj) on the plane 
correspond to (uiVj — UjVi) on the sphere, with (ui,Vi) being the spinor coordinates of Zi. 
The total angular momentum of a triplet (ijk) on the sphere is one-half of the maximum of 
the coherent state operator Sjjfc = f2 • (Lj + Lj + L^), where f2(</>, 9) is any direction on the 
sphere and Lj is the angular momentum of the i'th particle. (Lj + Lj + L^) commutes with 
any factor involving only particles i,j,k; i.e. all factors (u a v b - u b v a ) with a, b £ {i,j,k}. 
Translating this to the plane, we can find the maximum of Sijk acting on ^Hf almost 
by inspection. When all factors involving exactly one of Zi,zj or z^ are multiplied out, 
the result will be a polynomial with terms of the form zfz^z^. The maximum value of 
1/2{A + -B + C) is exactly the maximum value of the total angular momentum of the triplet 
(ijk). In particular, this maximum is 3iV^/2 — 4 for ^Hf, so it is certainly a zero energy 
eigenstate of Hjjf. 

It must still be shown that ^Hf is the unique ground state. To this end, we will use 
the method in Appendix A of Milovanovic and Read J64]]. Consider the behavior of the 



Haffnian prefactor in (4.47) as three particles (ijk) approach each other, the other particles 
remaining far away from the three. The Laurent series must contain terms of the form 
(zi — Zj) q%i (zj — Zk) qjk (zk — Zi) qki with q a b positive or negative integers. By continuity, the 
total function must contain this Laurent factor for any position of the particles. In order 
for ^hj t° be analytic, as it must be in the LLL, each q must not be smaller than —2, or 
Q = Qij+Qjk+Qki > —6. In particular, ^> Hf is annihilated by Pjjfc(3A r 0/2) if the inequality is 
strict, Q > —6. Another way to see this is on the plane: Pjjfc(3iV^/2) is the projection onto 
the closest approach of a triplet — all three particles are clumped at the north pole, taking 
the maximum L z value — which takes the form 5^ 2 \zi — Zj)S^ 2 ' (zj — Zk) on the plane. Thus, if 
Q = —6, the delta function interaction does not annihilate ^ Hf- Reducing the total angular 
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momentum (on the sphere) by one corresponds to restricting further the particles' closest 
approach, or increasing the lower limit of Q by one. In this way, we obtain the requirement 
Q > —3 in order for Pij k (L) to annihilate the ground state whenever L > 3iV^/2 — 4. The 
extreme case Q = —2 corresponds to the densest eigenstate. There are four possible such 
factors: 

1 1 (Zj - Zj) 2 (Zj - Zj) 

(Zi - Zj) 2 ' (Zi - Zj)(z k - Zj) ' (Zj - Z k ) 2 (z k - Zi) 2 ' [Zj - Z k ) 2 (z k - Zi) 

Symmetrizing these factors with respect to (ijk) leaves only terms like the first one. There- 
fore Hh f automatically requires a pairing structure of its ground state that is given by the 
Haffnian. 

Other zero-energy eigenstates are obtained by multplying in factors symmetric in the 
particle coordinates, being formally allowed since they can only increase Q. These states 
describe quasiholes and edge excitations and are less dense than the ground state since they 



contain added flux. They are enumerated explicity in the Section 4.2.c; here we verify that 
they are indeed zero-energy eigenstates of Hjjf. Without loss of generality, pick a definite 
triplet, (ijk) = (123), for convenience. There are several cases to check, corresponding to 
the possible terms appearing in some quasihole state ( 4.50[ ): (i) the pair (21,22) is broken, 
(ii) the pair involving z% and another particle, say Z4, is broken, (iii) both pairs (21,22) 
and (23,24) are broken, and (iv) neither pair is broken. For illustration, we check case (ii), 
the others being done similarly. Applying the generalization of Haldane's argument, the 
maximum degree of the triplet (123) is 2(N - 3) + 2(N - 3) + 2(N - 3) + 2n + (n - 2). The 
first three contributions come from the terms (2, — z a ), where i = 1,2,3 and a 7^ 1,2,3, 
the fourth term is due to the quasihole operator $(21,22,25,...), and the last takes into 
account the maximum orbital quantum number of the unpaired boson, z^~ 2 . Using the 
flux condition = 2(N — 1) — 2 + n, leads to one-half the maximum degree (or maximum 
total orbital angular momentum of a triplet) being 3^/2 — 4, which is consistent with 
requirement that it be less than 2>N^/2 — 3. Note that it has tacitly been assumed that 
n > 2; if n = 1 then the unpaired boson is in the zeroth orbital and half of the maximum 
degree is 3N,p/2 — 3 — n/2. The other three cases can be checked straightforwardly by such 
counting, proving that the quasihole states ( 4.5Q ) are annihilated by Hnf- 
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4.2.c Zero Energy Eigenstates 

In addition to the Haffhian, there are other zero energy eigenstates of Hjjf generated by 
adding n = 1,2,... flux quanta. The structure of these excitations is more complex than 
that of the familiar Laughlin quasiholes, having an infinite degeneracy in the thermodynamic 



limit. In previous work [Q, this degeneracy has been suggested to provide the necessary 
manifold of states for nonabelian statistics, which we discuss below. Proper construction 
of the zero energy states is useful for understanding the numerical spectrum of H B j, so we 
will go through them in some detail. 

The quasihole in the paired state, like in the Laughlin state, is generated by one flux, but 
it is built of two vortices at w\ and W2 each carrying one-half flux quantum. Let n be the 
number of flux added to the Haffnian, i.e. = 2(N — 1) — 2 + n, with w±, . . . , W2 n being the 
positions of the vortices. Denoting B/2 as the number of broken pairs in the Haffnian and 
{mi,77t2, • • • ,tib} as the quantum numbers of the orbitals into which the upaired bosons 
are placed, the explicit form for the manifold of zero-energy eigenstates is 

^ mi ,m 2 ,...,m B (z 1 ,z 2 , ■ ■ . , z N ; Wi, w 2 , ■ ■ ■ , w 2n ) = (4.49) 

B (N-B)/2 ^ f s 

Y T\ Z mk TT ^MB+2l-l),Z a{B+ 2l);W l ,...,W 2 n) yj, g ^ 2 ^ 

aeS N k=l ^ 1=1 {Z„{B+21-1) ~ Z u{B+2 l)) 2 i< - 

$ is the quasihole operator given by: 

n 

${z 1 ,Z2;W 1 ,...,W2n) = E YI(Z1 ~ W T (2r-l))(z2 ~ W r{ 2r)) ■ (4-50) 

t£S 2 „ r=l 

It can be verified that the $ mi ,m 2 ,...,m B {z\, Z2, ■ ■ ■ , zn', ti>i, w 2 , ■ ■ ■ , W2 n ) are in fact zero- 
energy eigenstates of H^f (Section [4.2.a| ). If no pairs are broken (B = 0), $ builds in 
two vortices — each within half of the bosons — for each of the n flux. Hence the interpreta- 
tion that the vortices carry charge 1/4, or l/2q for more general filling factors. The unpaired 
bosons are labeled by the orbital quantum numbers {mi, . . . ,m#}, which must satisfy the 
condition < mi < ... <m B < n — 2. The upper limit follows from the constraint on 
and the ordering is simply to avoid overcounting upon symmetrization. This is equivalent 
to putting B bosons into n — 1 orbitals, for which the multiplicity is 

B + n-2 
B 



8:] 



(4.51) 



Note that for n = 2 all of the bosons coming from broken pairs are in the lowest orbital, 
z°, and there are no unbroken pairs at n = 1. 

There is an additional degeneracy coming from the positions of the quasiholes them- 
selves, which is calculated by expanding in the io's using the elementary symmetric poly- 
nomials e m (w) = J2i 1 <...<i m w h ■ ■ ■ w i m - The e m have the property that, for m = 0, . . . , j, 
they form a basis for the algebra of all symmetric polynomials in j variables. Thus, zero- 
energy states may also be obtained as linear combinations of the e m . When <3? is expanded 
in this way we obtain all the symmetric polynomials in w\ , . . . , u>2n in which the degree 
of any w is no greater than (JV — B)/2. The total number of linearly independent states, 
for a fixed B and a fixed set of mj's, is at most the total number of linearly independent 
symmetric functions of w in the expansion of which establishes the upper bound on the 
positional degeneracy of the quasiholes. That this is also the correct degeneracy, without 
overcounting, is proven elsewhere [30]. We can now write down this number by regard- 



ing the vortices as some kind of bosonic particle, interpreting the e m (w) as the states for 
2n bosons occupying the (N — B)/2 + 1 orbitals {1, . . . , w^ N ~ B ^' 2 } (recall that $ appears 
(N - B)/2 times due to the product Ilifi~ B)/2 in (|450|)): 

/ 



(4.52) 



V 



(N - B)/2 + 2n 
2n 

Throughout this construction it has been tacitly assumed that N is even, so B is necessarily 
even as well. The construction for N odd proceeds with only slight modification, but the 
counting in ( [4.5l| ) and in ( f4.52 ) does not change, since B has the same parity as N. In either 



case, the total number of linearly independent quasihole states is obtained by multiplying 
the two combinatorial factors and summing over the allowed values of B: 

( 



E 



B + n-2 
B 



(N - B)/2 + 2n 
2n 



(4.53) 



To complete this description, one should check that all states in (4.50) exhaust all zero- 



energy eigenstates at fixed iV and n and that they are linearly independent. Since this is 
somewhat involved and is discussed at length elsewhere [pOj] , we will omit it here. 

As a special case, notice that for two added flux {n = 2) the Laughlin state, for bosons at 
v = 1/2, is recovered when all pairs are broken [B = N). Therefore, up to two flux quanta, 
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the Haffnian and the Laughlin states are degenerate eigenstates of Huf- It is convenient 
to adopt the composite boson interpretation of the quasihole positional degeneracy ( 4.52| ), 
Rewriting this combinatorial factor as 



(N - B)/2 + 2n 
(N -B)/2 



(4.54) 



affords an interpretation as (N — B)/2 composite bosons in 2n + 1 orbitals. On the sphere 
at n = 2 this is the correct degeneracy for an L = 2 angular momentum multiplet, indepen- 
dently of both iV and B. The unpaired bosons are all forced into the lowest orbital m = 0, 
which is manifested by the degeneracy factor (|4 .51 ) reducing to unity. Thus, the manifold 



of degenerate states is composed of N/2 states (one for each of B = 0, 2, . . . , N — 2) carrying 
angular momentum L = 2 and one Laughlin state (for B = N) carrying L = 0, which is 
a rotationally invariant. The main point of the expression ( 4.54j ) is that the zero energy 



eigenstates of Hhj contain both d-wave pairs and unpaired composite bosons in the m = 
orbital 

The degeneracy at L = 2 is a feature that is reminescent of the permanent. In that 
case, an analogous argument at n = 1 leads to a Laughlin state at v = 1 and a set of 
N/2 states carrying L = 1. The (N — B)/2 paired bosons were interpreted as spin waves 
|30|| , which we have established to condense into helical order (|4.35| ). Spin wave excitations 
are, of course, gapless, but one may expect that the L = 2 modes in the Haffnian are 
massive, being density excitations. However, we argue that this is not the case and H^f is 
compressible, indicating that it sits right on the transition from an incompressible state to 
one with non- uniform charge density order. 

4.2.d d-wave Pairing of Spinless Bosons 

In the previous subsection we showed that the lowest single particle orbital may be macro- 
scopically occupied by breaking pairs ( 4.54| ). A proper description of this system must 



therefore include a single particle condensate at the outset. As soon as there is such a 
condensate there is also isotropic scattering out of the superfluid, which gives rise to an 
s-wave pair amplitude. 
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On the torus or plane, the full Hamiltonian is given by 

K = E ( ^ - H c U + \ E ^(q)^ +|q i_i q c k _ i qCk , +|q (4.55) 

k \ / k,k',q 2 2 

We assume that the lowest state k = is macroscopically occupied with an amplitude 
$ = (co) and expand around this condensate: Ck — > $ + Ck- In addition there is the 
possibility of s- and d-wave pairing of the anomalous correlators (2kC_k). The resulting 
effective Hamiltonian is 



K eff = E' ^k4 C k + 1 ( A k4 c -k + A k: c -kCk 
k L 1 



(4.56) 



where the prime on the summation indicates that the k = term is to be omitted when 
> 0; otherwise the summation is over all wavevectors. The chemical potential is shifted 
by the (cqCq) contribution from expanding the interaction so that the single particle energy 
is now £k — k 2 /2m* — fj, + 2~l/(0)|<3?| 2 . For consistency, it is necessary that V^(0) is positive, 
which is a standard assumption in interacting Bose systems. It is also common to assume 
that the exchange contributions do not modify m* significantly. 

The pair order parameter must satisfy the self-consistent gap equation in the presence 
of a condensate: 

Ak = $Mk)-^EV(k-k')| k f (4.57) 

k' k ' 

To make sense of this complicated integral equation, we shall assume that Ak is dominated 
by the (s+d)-wave symmetry and, furthermore, that its leading order behavior is 

i A, ~ A s 

A k = A£ + A£ with { _ (4.58) 



where A s,d are constants. Recall that the interaction can be expanded into angular momen- 
tum eigenstates (4.26) and that each channel has the asymptotic behavior V[(k, k 1 ) ~ k l . In 



other words, it is implicit in the above form ( ^4.58|) that the long distance property of each 
order parameter is dictated by its respective angular momentum channel. 

The standard procedure is to diagonalize K e jf. This is the familiar Bogoliubov trans- 
formation that replaces the original particles and holes by quasiparticles ak: 

a k = u k c k + «kcL k (4.59) 
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Commutation relations are preserved if |uk| 2 — |^k| 2 = 1 and symmetry requires that Uk 
ti_k and ^k = v_k- The following solutions for and put K e ff into the form K e ff 
Ek-^ka^k : 



4 = US-VlA 

-k = U4r + l) ( 4 - 6 °) 



1 A A 



2 \E k 

El = e~ |A k | 2 . (4.61) 

There is a kind of gauge freedom in that both and can be multiplied by a k— dependent 
phase factor, without changing the physics. We adopt the convention that «k is real and 
positive, while carries the possible d-wave phase factor. 

If both s- and d-wave amplitudes are nonvanishing, an interesting feature emerges im- 
mediately. The quasiparticle spectrum contains an anisotropic piece coming from the cross 
terms in |Ak| 2 , 

El = £k - |A S + k d k 2 e- 2i ^\ 2 , (4.62) 

where k is the magnitude of and <fi is the polar angle of k. Choosing the constants A s,d to 
be real for the moment, the anisotropy is 

- 2A s A d k 2 cos 2cp (4.63) 

which has two minima at <p = 0, n. There is a relative U(l) phase degree of freedom 
between the pairing amplitudes, which is equivalent to an SO (2) rotation of the minima. 
In other words, the transformation A s — > e 2ia A s is the same as 4> — > <fi + a. It is interesting 
that this anisotropy arises not from an attraction in I^(k), which is isotropic, but rather 
from the coexistence of two distinct pairing symmetries in momentum space. 

A unique ground state wavefunction exists for each solution of the gap equations. Its 
generic form is like that of the permanent ( [4.15 ) but with an allowance for the single particle 



condensate 

1 / 

—7 exp < 



exp j$ C t + £'f) k 4<i k j |0) . (4.64) 



The normalization is the same as in Eq. ( 4.28; ) but with an extra factor of exp(|$| 2 ). 



The reciprocal space pair wavefunction, g^, can be expressed in terms of the Bogoliubov 
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parameters Uk and as 



5k = — (4.65) 



The generic form of the iV-body paired wavefunction has already been discussed in Eq. 
( 4.29| ) and the Haffnian exists only when $ = since it is a pure paired state. Furthermore, 
the asymptotic behavior of the real space pair wavefunction g(r) ~ 1/z 2 is obtained when 

The iV-body state of composite bosons constructed in this way is exactly the Haffnian 
prefactor of Eq. ( 4.47] ). And, the same remarks following Eq. (4.31), concerning the 
validity of this mean field approximation, apply here. 

Unfortunately, unlike the permanent case, the occupation number (rik) = (c k Ck) does 
not contain any new information since it is constant: 

M - ■ 

However, the required form of does constrain the parameters at which the Haffnian can 
exist. At small k, the asymptotic behavior of c/k ( [4. 66 ) is consistent with A s = only if Uk 
and i>k have the same asymptotic behavior. 

To understand where this point lies we must understand the phase diagram as derived 
from the self-consistency condition ( [4.571 ), It is simplified by assuming the asymptotic 
expansion for Ak in Eq. fl4.58| ) and by expanding V(k — k') into angular momentum 
eigenstates as before ( |4.25| ). This decomposition separates the BCS gap equation ( 4.57 ) 
into two coupled integral equations, 

A* = *V( k )4y>„(^) A - + Ad r^' (4.07) 

AV - -l E V_ 2 (M') ^f±^ ■ (4.68, 

The assumption that the s-wave amplitude is constant then implies that k = in the first 
equation to leading order. The coefficient Vo(k,k') reduces to V{k!) at k = 0, and the 
leading order term of V-2(k,k') is proportional to k 2 so the second equation is consistent 
(c.f. Eq. ( 4.26| )). Further simplification is possible by choosing a particular gauge for the 
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phases of the condensates. A convenient choice is to pick a real A d and to absorb the phase 
of A s into (/)'; this only rotates the inhomogeneity of the spectrum Q4.63|) . Since the phase 
of $ is locked to that of A s , all order parameters can be chosen to be real and positive. 
Finally, since the inhomogeneity is even under parity, (ft' — > —(ft', the odd terms s'm2(ft' in 
the numerator cancel and only cos 2<ft' appears. 

The crudest criterion for the solutions of the gap equations is whether or not $ vanishes. 
Consider first the situation when $ = and the fluid consists of pure pairs. In principle both 
s- and d-wave pairing can coexist since the gap equations allow this generally. However, if we 
make the reasonable assumption that V(k) is positive everywhere then only A d is allowed. 
This may be seen from a closer inspection of the ^'-dependent part of the integrand in 
( 4.67| ). Let us split it into two disjoint pieces characterized by cos 20' negative or positive. 



The corresponding E^i is always larger in the former case because the inhomogeneity ( 4.63 ) 
has a negative sign. Therefore, the total contribution of the A d term to the right hand side 
is always negative. That is to say 

deft' V(k')^^>0 . (4.69) 
Jo Ek> 

Since A s is also positive in our gauge there is no consistent solution of ( 4,67| ) when $ = 0. 



The only possibility is a pure d-wave gap that satisfies ( 4.68 ) with A s = 0. In this phase 
H < and 2Eo = 2\fi\ is the gap to breaking a pair (see also Thus, when there 

is no single particle condensate the system is specified by the following gap equation and 
quasiparticle spectrum: 

1 A d h' 2 

^ 2 = — 



2^ By 



m* 



— ) 2 -(n 

2m* J V J 



(4.70) 



A d is real in this gauge. 

On the other hand, when $ > there is a new set of considerations, which require 
Vo(k,k) = for consistency. Firstly, in order for K e ff to be a stable approximation to 
the full Hamiltonian, (cq) must be a minimum of the free energy. The resultant constraint, 
sometimes known as the Gross-Pitaevsky equation, pins the chemical potential to the single 
particle condensate by fi = ^(O)^ 2 (in the gauge where $ is real). Furthermore, the 
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spectrum is dominated by the gapless Goldstone phonon mode at low k, requiring Eq = 
e - |A | 2 , or 

- li + 2U(0)$ 2 = A s . (4.71) 

In conjunction with the pinning condition li = V(0)<I> 2 , this implies that A s = ^(O)^ 2 . 
But, using the positivity condition in Eq. ( 4.69[ ), this is consistent with the gap equation 



(|4.67| ) only if we reinstate the k-dependence of A s or if Vo(k, k') vanishes. We shall assume 
the latter as it is much more tractable. Ultimately, this condition is traced back to the 
asymptotic behavior of the gap function (4.58), which is somewhat restrictive. However, we 



now have a tractable and fully consistent system. Summarizing these results for $ > 0, 

A s = V(0)$ 2 = li 
A d k 2 = -lyy_ 2 (M') A * C ° S20/ + Adfc/2 (4.72) 



^ = 2 M^- Adcos2 ^ + a )) fc2 + 



« A 2 



k l 



2m* , 

We have used the gauge in which all order parameters are real, which shifts by a, the 
phase of A s . The single particle condensate constraint A s = li has also been used in the 
last expression. 

We can now map out the phase diagram in the space of \x and A d by looking at the 
spectra in Eqs.( 4.70 ) and ( |4.72|) . Fig. |4^ illustrates the possible phases. The dotted line 



indicates the instability of both spectra when A d > l/2m*. For clarity, in each of the 
four regions there is an inset of the typical spectrum E^. The spectrum in region (IV) is 
presumably stabilized at higher momenta (which is represented by the dashed curve) so that 
the instability is really a kind of "d-wave roton". The spectrum in region (III) is clearly 
unstable, indicating a new phase. As li increases from (IV) into (III), the pair gap and 
probably the roton gap, too, go to zero. Although we have no explicit calculations, it is 
probable that the phase in (III) is a charge-density-wave at the roton wavevector. On the 
other hand, there are no instabilities in regions (I) and (II). (I) is a pure, isotropic, d-wave 
condensate, in which the gap is half of the energy required to break a pair. As the gap 
closes, and we move into (II) at fixed A d , a single particle condensate develops so that all 
three condensates, A s , and A d exist; accordingly, there is a linear phonon mode and an 



90 



anisotropy due to the coexistence of both pairings. As A d increases from (II) to (III), the 
anisotropy is strong enough that the gap closes at finite momentum, in addition to k = 0. 
This is consistent with the transition from (IV) to (III). 

The Haffnian itself can only exist at (i = 0, which can be seen from the asymptotics of 
and Since the pair amplitude is a pure I = —2 eigenstate at k — > when the Haffnian 
is the asymptotic form of the many-body wavefunction, A s must vanish. Furthermore, 
the requirement that \g^\ = \v^/u^\ — > 1 as k — > restricts the relative behavior of Kk 
and t>k, and the explicit solution in Eq. ( [4.61 ) implies that \i = 0. The thick line in the 
figure marks the transition region, fi = with a stable spectrum (A d < l/2m*), where the 
Haffnian represents the long range behavior. E. Rezayi \76\ has analyzed numerically the 
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Figure 4.3: Phases of bosons in two dimensions. 



Haffnian Hamiltonian in eqn. |4.47> He found that for Vq negative the ground state is a 
paired state, while for Vq positive it is a Laughlin state, which is a Bose condensate. This 
is consistent with the Haffnian being on a phase boundary. 
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Chapter 5 

Fermion Pairing: Quantum Hall 
Effect for Spin 



In this Chapter, we provide a detailed derivation of the Hall conductivity for spin transport 
in the d- and p-wave pairing of fermions. In the FQHE, the fermions are really composite 
fermions, and we ignore gauge field fluctuations. This is equivalent to showing that the 
induced action for the system in an external gauge field that couples to the spin is a Chcrn- 
Simons (CS) term. In the d-wave case, the system is spin-rotation invariant, so we obtain 
an SU(2) CS term, while in the p-wave case, there is only a U(l) symmetry, so we find a 
U(l) CS term. In both cases, the Hall spin conductivity is given by a topological invariant. 
Within the BCS mean field approach, using a conserving approximation, this topological 
invariant is the winding number of the order parameter in momentum space and is an 
integer, which is the statement of quantization. We argue that the quantization in terms of 
a topological invariant is more general than the approximation used. 



5.1 BCS Hamiltonian 



Considering first the spin-singlet paired states, we use the Nambu basis where the symme- 
tries are transparent. Define 



V2 



(5.1) 
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with 

/ 



(5.2) 



so that \I> transforms as a tensor product of particle-hole and spin-space spinors. We consider 
an interacting system and approximate it as in BCS theory, then with a minimal coupling to 
the gauge field, we use a conserving approximation to obtain the spin response. In Fourier 
space, we should note that 

*k = ^( °\ I- 




,t 

In the Nambu basis, the kinetic energy becomes (with K = H — fj,N) 

K o = X^k(4t c kt +4| c k|) 

k 

= E^tfo®^, (5.4) 

k 

where = |k| 2 / (2m)— /i is the kinetic energy, containing the bare mass m, and the products 
in the spinor space are understood. Products like a z (£> I act on the Nambu spinors, with 
the first factor acting in the particle-hole factor, the second in the spin-space factor. The 
interaction term, for a spin-independent interaction V, is 

*int = \ E V n *t+q(^ ® ^)*k*L_q(^ ® ^ = • (5.5) 
kk'q 

Here the colons : . . . : denote normal ordering, that is all the c^s are brought to the left. 
In the BCS-extended Hartree-Fock approximation, the effective quasiparticle Hamiltonian 
(for later reference) is 

K cS = E*k[&(<7*®/)+ReA k (<T x ® J) 

k 

-ImAk^^/pk. (5.6) 

This is for singlet pairing, where A_ k = A k , and not just for d-wave. Here £ k is plus 
the Hartree-Fock corrections. If we define a vector 

E k = (ReA k ,-ImA k ,6c) (5.7) 

then the quasiparticle energy £" k = |E k |, and 

K oS = E *k( E k • ° ® /)*k- (5-8) 
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5.2 Spin Response in a Conserving Approximation 

In the Nambu notation, it is clear that K = Kq + i^int, and K e s, are invariant under global 
SU(2) rotations that act on the spin-space, that is the second factor in the tensor products. 
The spin density, the integral of which over all space is the total spin and generates such 
global transformations, and the spin current densities are given by 

^o(q) = ^£<- q / 2 ('®^)*k +q /2 (5.9) 
k 

J "(q) = ^E^L q / 2 (^®^)*k +q/ 2, (5.io) 

k 

where i = x, y is the spatial index, and a = x, y, z is the spin-space index. Spin conservation 
implies the continuity equation, as an operator equation, 

dJ^/dx^ = 0, (5.11) 

where fi = 0, x, y, and the summation convention is in force. 

So far we have not introduced a gauge field for spin. Since the spin is conserved locally, 
we can turn the symmetry into a local gauge symmetry by introducing an SU(2) vector 
potential, and making all derivatives covariant. The effect on K is to add the integral of 




(5.12) 



The gauge field is to be used solely as an external source, with which to probe the spin 
response of the system, and then set to zero. 

If we now consider integrating out the fermions, then we can obtain an action in the 
external gauge fields, which can be expanded in powers of A a . The zeroth-order term is 
the free energy density, times the volume of spacetime, and the first-order term vanishes by 
spin-rotation invariance. The second-order term corresponds to linear response: the second 
functional derivative with respect to A*, at A® = 0, yields the (matrix of) linear response 
functions. In particular, the spatial components yield the conductivity tensor in the usual 
way. Therefore we consider the imaginary-time time-ordered function, 

K = ^(J a Jq)Jt(-q)), (5-13) 
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where time-ordering is understood, and from here on we use a convention that p, q, etc. 
stand for three- vectors p = (po,p), and further po = iuj is imaginary for imaginary time. 
For [i = v = i = x or y, an additional "diamagnetic" term nS ab /Am is present in 11^, which 
we do not show explicitly. As consequences of the continuity equation and the related gauge 
invariance, Ii^ u must be divergenceless on both variables, q^H^u = Qu^^u = 0. To maintain 
these when using the BCS-Hartree-Fock approximation for the equilibrium properties, one 
must use a conserving approximation for the response function, which in this case means 
summing ladder diagrams (compare the charge case in Ref. |2l|| , pp. 224-237). We used the 
identical method in treating composite bosons at v ^ 1 in Chapter ||, Chapter |2[ 

One begins with the BCS-Hartree-Fock approximation, which can be written in terms of 
Green's functions as (we consider only zero temperature, and / dpo is along the imaginary 
Po axis throughout) 

G-\p) = Po -e p a z ®I-i:(p), (5.14) 
d 3 k 
(2tt) e 



S Cp) = i I 7^{<y z ®I)G{k){a z ®I)V{k-q). (5.15) 



Note that G(p) and S(p) are matrices acting on the tensor product space. The equations 
are solved by 

G~ 1 (p) = R> - E p • cr ® /, (5.16) 

(we write 1 for 1 J) as one can also see from the effective quasiparticle Hamiltonian K e g, 
and A p obeys the standard gap equation. 

In the response function, the ladder series can be summed and included by dressing one 
vertex, to obtain (again not showing the diamagnetic term) 

n£(g) = -iJ-^tT[rf(p,p + q)G(p + q) 

xr b Ap + q,p)G(p)], (5.17) 



where 7^ is the bare vertex, 



lt(p,P + q) = - {P ^ Q)i <T z ®<Ta, (5.19) 
2m 



95 



and is the dressed vertex satisfying 

f d 3 k 

K(P + q,p) = lt(p + q,p)+ij 7^31^ ® IG ( k + l) 

xT h v (k + q,k)G{k)a z ®IV{p-k). (5.20) 

At small q, we can obtain useful information about this function from the Ward identity 
that results from the continuity equation. The particular Ward identity we use here is an 
exact relation of the vertex function to the self-energy, and the conserving approximation 
(the ladder series) was constructed to ensure that it holds also for the approximated vertex 
and self energy functions. 



Following Schrieffer's treatment 21], we consider the vertex function with external legs 
included: 

A£(n, r 2 ,r 3 ) = (JJl(r 3 )* (n)& (r 2 )) , (5.21) 

for spacetime coordinates r%, r 2 , r^. Applying d/dr^ to both sides and using the operator 
continuity equation, we obtain the exact identity in Fourier space 

q^(p + q,p) = h ® a a G- l ( V ) - ^G-\p + q)I ® a a . (5.22) 

Since G _1 is trivial in the spin-space indices, it commutes with I ® a a . Hence at q — > 0, 
the right-hand side vanishes, so T(p + q,p) has no singularities as q — > 0. This differs from 
the charge case, for example, where this calculation (using the ladder series approximation) 
leads to the discovery of the collective mode p9[] . Since the spin symmetry is unbroken, no 
collective mode is necessary to restore this conservation law, and so there is no singularity 
in the vertex function for spin. 

One can verify that the Ward identity is satisfied using the BCS-Hartree-Fock G _1 and 
the ladder series for T. At q = 0, this yields the important results 

T>,P) = -\d^G- l {p)I ® a a , (5.23) 



or explicitly, 



T? = -~diG-\p)I®a a , (5.24) 
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where d\ and <9 M stand for d/dpi, d/dp^ from here on, and the extra minus in the first 
relation is consistent because implicitly q^T^ = gcTo — qiTi. 

We now calculate II at small q. To zeroth order, use of the Ward identity shows that 
the J-J function gives zero, except when ji = v = i. In that case, it reduces to a constant 
that cancels the diamagnetic term also present in just that case. Hence we require only the 
part first-order in q. In the expression for II above, we first shift p — ► p — ^q, so that q no 
longer appears in any bare vertices, but does appear in the Green's functions on both sides 
of the ladder, between the rungs which are the interaction lines. Hence to first order, we 
obtain a factor ±i<9 M G = =f^GdG~ 1 G in place of G in one position in the ladder. Since 
there may be any number of rungs (including zero) between this and either of the vertices 
at the ends, the terms can be summed up into a ladder dressing each vertex, evaluated at 
q = 0. Hence we obtain to first order 

n£(9) = ~J -^tv[r° (p, P ) qx d x GTUp,p)G(p) 

-K(p,p)G(p)T b u (p,p)q x d x G] . (5.25) 

Using the Ward identity, this becomes 

n£(ff) = ±qx f -^tr {{I ® * a ){I ® aJGdpG- 1 

x \Gd\G~ 1 , Gd v G~ l \^ (5.26) 

Since the G's are independent of the spin-space indices, the explicit cr's factor off, and 
the result is 5 ab times a spin-independent part. The latter can be simplified using the 
BCS-Hartree-Fock form of G, by writing the latter as 



The spin-independent factor contains e^x since it is antisymmetric in these labels. Keeping 
track of the signs, we find for the quadratic term in the induced action 

with M. given by the topological invariant 

M = j l^-Ep • (diE p x djE p )/El. (5.29) 
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The right hand side is exactly the Pontriagin winding number to, and is an integer as long 
as E is a continuous, differentiable function of p; it is 2 for the d-wave case. 

To ensure SU(2) gauge invariance, the CS term should include also a term cubic in A, 
with no derivatives. For this term we evaluate the triangle one- loop diagrams with three 
insertions of J, with each vertex dressed by the ladder series. Setting the external momenta 
to zero, the Ward identity can be used for all three vertices, and the result can be seen to 
be 

1 f d 3 p 



n^ c A (0,0) = --J-^tr^I^^Gd.G- 1 



(2vr) 3 

x {(/ ® a b )Gd u G~\ (/ ® a c )Gd x G- l }\ . (5.30) 

The anticommutator { , } arises since the result must be symmetric under permutations of 
the index pairs fi, a, etc. The product a a abcr c , when traced over the spin-space indices, yields 
a factor 2ie a i, c , which is antisymmetric, and so the remainder must contain e^ u \ to maintain 
symmetry; the rest of the structure is the same as before. Hence the full result is the SU(2) 
CS term, which we write in terms of the 2x2 matrix vector potentials = \a a A a ^ 

k f q . 2 



^ J (fxe^MA^Ax + —A fl A u A\). (5.31) 



Here k is the conventional notation for the coefficient of such a term, in this same normal- 
ization; if we wished to quantize the theory by functionally integrating over A, we would 
need k = an integer. In our case k = M./2 = 1 for d-wave. 

For the spin-triplet case with an unbroken U(l) symmetry, we must use the fact that 
A_k = — Ak. For example, in the two-dimensional A-phase, as occurs in the 331 state 
in the double-layer FQHE system with zero tunneling, the pairs are in the isospin S z = 
triplet state + and the U(l) symmetry generated by S z is unbroken; we recall that 
the underlying Hamiltonian is not assumed to have a full SU(2) symmetry. The effective 
quasiparticle Hamiltonian becomes, in the Nambu-style notation, 

K eS = 5Z*lKk(^®^)+ReA k (<7 x (8)^) 

k 

- Im A k (a y ® a z )\ * k . (5.32) 

The U(l) vector potential A^ couples to S z , and the vertex functions contain I®a z , which 
commutes with the BCS-Hartree-Fock Green's function G. The tensors appearing in the 
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three terms in K e g obey the same algebra as the three in that for the spin-singlet case 
(where they were trivial in the second factor), and as in that case commute with I ® <j z . 
Consequently, the derivation for the induced action to quadratic order in is similar to 
that for the SU(2) singlet case above, and the traces in the Nambu indices can be carried 
out with the same result as before, to obtain the abelian CS term 

1 f dA 

and no cubic term. In this case, M is again given by the winding number m which is or 
±1 in the p-wave strong and weak-pairing phases (respectively) discussed in this chapter. 



5.3 Discussion and Generalizations 

We note that the effect of the vertex corrections we included as ladder series is to renormalize 



the q = vertices as shown in eq. ( 5.23 ) for the spin-singlet case, and use these in one-loop 
diagrams with no further corrections. This corresponds to the minimal coupling p — > p — A 
in the action, as one would expect by gauge invariance. If we assume such a coupling, 
and treat the low-energy, long-wavelength theory near the weak-strong transition as Dirac 
fermions with relativistic dispersion and minimal coupling to the external gauge field, then 
the expression for M as an integral over p covers only half the sphere in n space, and 
we would get ±1 (d-wave), ±1/2 (p-wave). The missing part results from the ultraviolet 



regulator in the field theory version of the calculation [78], or from a second fermion with 



a fixed mass in some lattice models [77|. In our calculation, the remainder is provided by 
the ultraviolet region, where Ak — * as k — ► oo. At the transition, p = 0, the map is 
discontinuous and covers exactly half the sphere in the p-wave case, so Ai = 1/2, as in 
other problems. In the d-wave case with rotational symmetry, the value of |wk/i*k| as k — > 
is nonuniversal, and hence so is the value of a s xy at the transition. This is a consequence 
of the non-relativistic form of the dispersion relation of the low-energy fermions in this 
case. We may also note that for a paired system on a lattice, as in models of high T c 
superconductors, a similar calculation will give an integral over the Brillouin zone, which 
is a torus T 2 , instead of the k plane which can be compactified to S 2 . But maps from T 2 
to S 2 are again classified by the integers, and the integer winding number is given by the 
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same expression, so quantization is unaffected. 

We can also argue that the quantization result away from a transition is exact in a 
translationally-invariant system, at least in all orders in perturbation theory. For this we 
use the form in eq. ( |5.26| ) or ( |5.30| ), where the Ward identity for the vertex has been used. 
Diagrammatically, it is clear that the exact expression can be similarly written, using the 
exact (i.e., all orders in perturbation) Green's function and vertex function. (This is also 
true when the CS gauge field interaction is included.) The Ward identity that relates them 
is exact, and the result for a s xy is of the same form as shown. The next step, the frequency 
integrals, cannot be done explicitly in this case, because the precise form of the Green's 
function is unknown, and the analogs of £k, Ak (or of u^, v^) do not exist. The latter 
do not exist because in general the poles in the Green's function, which would represent 
the quasiparticles, are broadened by scattering processes, except for the lowest energies for 
kinematical reasons. However, the form in eq. ( |5.26 ) is itself a topological invariant, as we 
will now argue. As long as there is a gap in the support of the spectral function of G, G(p) 
is continuous and differentiable on the imaginary frequency axis, and tends to I I/po 
as po — > ±zoo. Thus G _1 exists and never vanishes. Considering the spin-singlet case for 
convenience, the spin-space structure is trivial, so we may perform the corresponding traces, 
and then G or G~ l is a 2 x 2 matrix, with the same reality properties on the imaginary 
Po axis as in the BCS-Hartree-Fock approximation. (The spin-triplet case should work out 
similarly, because of the algebraic structure already mentioned.) It thus represents a real 
non-zero 4-component vector, in R — 0, which topologically is the same as S 3 . S 3 is obtained 
by dividing G by its norm, (trG^G) 1 / 2 , and the normalized G is a 2 x 2 unitary matrix with 
determinant —1, so it lies in S 3 . The k space can be compactified to S 2 as before, and the 
frequency variable can be viewed as an element of the interval I = (— 1, 1), so the integral 
is over S 2 x I. However, since the limit of the Green's function as po — > ±ioo for fixed k is 
independent of k, we can view this as simply S 3 . Thus we are dealing with maps from S 3 
to S 3 , the equivalence classes of which are classified by the homotopy group -k^(S 3 ) = Z. 
The integral we have obtained simply calculates the integer winding number or Pontryagin 
index of the map, when properly normalized (G can be normalized to lie in SU(2) without 
affecting the integral) . This establishes the quantization of translationally-invariant 
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system with a gap, at least to all orders in perturbation theory, and probably can be made 
fully non-perturbative (as the Ward identity is already). 
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Chapter 6 



Adsorption on Carbon Nanotubes 



In this chapter we switch gears into one dimension and consider adsorption on nanotubes 



as discussed in the introductory chapter, Section |L3 

6.1 Nanotube geometry 



Fig. 3.1 illustrates the way in which a nanotube is obtained by wrapping a graphite sheet. 
The hexagons are at positions R n ,m = na^ + mo_ , where a± are primitive lattice vectors 
of the honeycomb lattice. The standard convention is to identify i?o,o = Rn,m, and to 
simply label the tube (N, M). The case (N, N) is known as the armchair tube, (N, 0) is the 
zig-zag, and all others are chiral. There is a geometric frustration whenever the wrapping 
destroys the tripartite nature of the infinite sheet, which occurs when (N — M)mod 3 is 
non-zero (this criterion is familiar in the context of electronic conductivity [fH|] ) . Thus, both 
zig-zag and chiral tubes can be frustrated geometrically, whereas armchair tubes cannot. 
The adsorption sites form a triangular lattice wrapped on the cylinder, which is shown in 
Fig. |6.2| for the (7, 0) zig-zag. In the following section we specify the Hamiltonian and 
consider the simplest (classical) limit in which intersite tunneling of adatoms is prohibited. 
We will then turn on the hopping perturbatively (the quantum case). 
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Figure 6.1: An example of wrapping of the graphite sheet to make a (2, 1) tube. a± are the primitive 
lattice vectors of the honecomb lattice. The solid rectangle is the primitive cell of the tube. The 
tube can also be built up by stacking the dotted region along the axis with the solid dotted lines 
identified. Also shown is the tripartite lattice labeling A, B, C. 
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Figure 6.2: Adsorption sites on a (7,0) zig-zag nanotube 
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Figure 6.3: Phase diagram of the (N, M) tube 

6.2 The Hamiltonian and Classical Limit 

When the adsorbed gas is a hard-core boson, the lattice gas is defined by the Bose-Hubbard 



Hamiltonian [ 42 



H = -t ^2 Kbj + tybi + V ^2 n i n i ~ M X! Ui ' 

(ij) (ij) » 



(6.1) 



where is the boson density at site i, V is the nearest neighbor repulsion and t is the 
hopping amplitude. The occupation numbers rij are restricted to 0, 1 by the hard-core 



condition. There is a familiar Heisenberg spin representation [81], which identifies Sf 
1/2, Sf 



n 



b\ , and S i 



hi. The Hamiltonian is thus 



7i = -2tJ2 s-sj + sfs] + v Sf S z j- h sf , 

(ij) (ij) i 



(6.2) 



where Sf = rij — 1/2 and H = fj,— 3V is an effective external magnetic field. Throughout the 
paper we will use the spin and density representations interchangeably. The spin models 
obtained in this way are similar to recent examples of "spin tubes" [48]. 

The Ising limit of the spin models, t = corresponds to the case when hopping is 
forbidden, and already contains many interesting features. We start the analysis in this 
regime, obtaining the phase diagram as a function of the magnetic field, and then consider 
quantum fluctuations perturbatively in t/V. We summarize our results first. 

The phase diagram in the temperature-magnetic field plane of a typical tube is shown 
in Fig. |6.3| . When the index q = (N — M)mod 3 is 1 or 2, we find four lobes (solid 
lines), corresponding to two plateaus with magnetizations m_ < 1/3 and m + > 1/3. Here, 
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we use the standard Ising notation in which spin is ±1. Note that the plateaus are real 
phases only at zero temperature because the tube is one-dimensional. At finite temperature, 
the boundaries should be interpreted as crossovers. Nonetheless, deep within a lobe, at 
&bT <C V, the magnetizations are well-defined. Specifically, for q = 1, we obtain the exact 
expressions 

1 ( 2 \ 1 , 

m + = - I 1 H — — m_ = - 1 



2M + NJ 3 V 2N + M 

2M 



H c = 4 - — — V (6.3) 

V n + mJ y J 

The complementary case of q = 2 is obtained by interchanging N «-> M. On the other 
hand, those tubes without geometric frustration (q = 0) behave similarly to the flat sheet 
(dotted lines) which has only two lobes with magnetizations ±1/3 [^TJ. In the flat sheet, 
the dotted lines are second order phase transitions in the universality class of the Potts-3 
models p2| . In our wrapped case, as the tube perimeter approaches the flat sheet limit, 
one expects that the geometric frustration becomes irrelevant. Indeed, as N or M — ► oo, 
m + and m_ squeeze 1/3 as the inverse of the tube diameter and become indistinguishable. 
Beyond the lobes, where the field is strong enough to overcome all nearest neighbor bonds 
(l-ffl/V > 6 at ksT = 0), the tube is fully polarized. The filling fractions are obtained 
from the magnetizations by m = — 2(n — 1/2). The phase diagram, however, is more 
easily visualized in terms of spin since spin reversal, m <-> — m, corresponds to particle-hole 
symmetry, n <-> 1 — n. 

We have verified this prediction numerically by transfer matrix methods [pi]] for zig-zag 
tubes up to N = 11 and for the chiral tubes up to N + M = 7.F] The transfer matrix rows 



for a sample tube are delineated by dotted lines in Fig. 3.1, It should be noted that similar 
transfer matrix calculations have been carried out for the special case of unfrustrated zig-zag 
tubes (q = 0) [p^] . The motivation in these earlier works was a finite size scaling analysis 
of the solid phases on flat graphite. 

Although 7 is probably too small to be physical, we believe that the arguments in this 



paper generalize to any tube. In Fig. |6.4j we display sample data for two zig-zag tubes with 
*The programming used the Gnu implementation of Fortran 77 on a PC. The main stumbling block is 
very large contributions to the partition function at low T; the LAPACK routine library was used to handle 
numbers outside ordinary machine range. 
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Figure 6.4: Magnetization and entropy per site at ksT — 0.05V 

different q: (7,0) and (8,0). The magnetization curves show clear plateaus whose values 
and transition fields match those predicted by Eq. fl6.3|) . By increasing the temperature 
and following the evolution of the plateaus, we generate the phase diagram above. 

We find that a rather interesting feature of the zig-zag (iV, 0) tubes emerges, making 



them exceptional. The insets in Fig. 3.4 indicate an extensive entropy at zero temperature, 
which has plateaus, too. Upon enumerating the degenerate space explicitly (Section |Q| 
below), we shall show that the entropy is exactly s = (ln2)/iV and that it occurs in m+ 
for q = 1 and in m_ for q = 2. In the presence of hopping, the non-degenerate plateaus 
retain their gaps, whereas the degenerate ones become correlated states with a unique 
ground state and gapless excitations. More precisely, conformal invariance develops and the 
effective theory has central charge c = 1 with a compactification radius, R, quantized by 
the tube circumference, R = N . 

In order to understand the magnetizations and nature of the geometric frustration, it 
is more intuitive to use the original bosonic picture. As a result of hard-core repulsion 
on the infinite graphite sheet, the m = 1/3 plateau corresponds to filling one of the three 
sublattices, A, B or C, of the triangular lattice. This configuration minimizes the repulsion, 
VriiJij, while maximizing the filling, [in. It is natural to try the same for nanotubes, as we 



illustrate in Fig. 6J3 for (5, 0). Upon wrapping, however, the thick vertical lines are identified 
and the lattice is no longer tripartite. In fact, the number of sublattice sites is no longer 
equal, and there is a mismatch along the thick line, which we term the "zipper" . On the left 
we fill the A sublattice, obtaining the filling fraction n+ = 2/5, and on the right either B or 
C may be filled with the result that ra_ = 3/10. For general (N, 0) there are 2N hexagons in 
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Figure 6.5: Fillings and zipper of the (5,0) zig-zag tube. n± corresponds to m T 

the unit cell, and the filling fractions are n + = \2NJS\ /2N and n_ = [2./V/3J /2N, where \x] 
and [x\ denote the larger and smaller of the two bounding integers of x, respectively. The 
magnetizations in Eq. ( |6.3j ) follow directly by using the correspondence m = — 2(n — 1/2). 
Furthermore, due to the sublattice mismatch, the number density of adjacent particles, rib, 
may be non-zero. In the case of (5, 0), there are two broken bonds per unit cell in n + , and 
none in n_. This result generalizes to any q = 2 zig-zag tube: rib+ = 2/2N and rib- = 0. 
For q = 1, the argument goes through as before, except that rib+ = 1/2N. We summarize 
this compactly by rib + = q/2N. 

Substituting these fillings into the Hamiltonian Q6.1| ) yields two energies per site, e± (p) = 
Vrib±— fin±. The transition occurs when these levels cross: e + = e_, or 

q \2N/3] [2N/3\ 

2N ~ = '"-IN- (6 - 4) 

Solving for [i and using the correspondence H = ji — 3V gives precisely the critical field 
in Eq. ( |6.3| ). In particular, this explains why there are exactly two independent plateaus. 
Note that, for the special case of the zig-zags, the critical field depends only on q and not 
on N per se. 

In the above analysis, we have made only one assumption, namely that the zipper runs 
parallel to the tube axis. In general, the zipper may wind helically around the tube or 
wiggle sideways. However, in all the cases that we considered, the straight zipper has the 
lowest energy, and moreover, our transfer matrix computations, which are blind to this 
assumption, are consistent with our analysis. 
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The chiral tubes are different. Due to their geometry the zipper is forced to wind, but, 
again, we find that the choice of the straightest possible zipper reproduces our numerics for 
N + M up to 7. The determination of the fillings and level crossings is much more involved 
than that of the zig-zag, and we leave it for a more detailed paper. In any case, our analysis 
reveals that the plateaus in a chiral tube are not macroscopically degenerate, so that the 
zig-zags are at a special degenerate point. 



6.3 Macroscopic Degeneracy and Quantum Fluctuations 

Having understood in detail the Ising limit, we now turn on a small hopping, t <C V, that 
introduces quantum fluctuations. Deep within a plateau, the substrate is maximally filled 
since adding a particle increases n&. Consequently, all plateaus begin with a classical gap 
of order V, and we work in the Hilbert space of the classical ground states. Those plateaus 
which have only a discrete symmetry must retain their gaps, but the macroscopically de- 
generate plateaus are more complicated. 

Let us reconsider the n + filling of the (5, 0) tube in Fig. |6.5|, Notice that a particle may 



hop laterally by one site without changing n&, as we illustrate in Fig. 6.6, left. Imagine 
building a typical n + state layer-by-layer from top to bottom, with a total of L layers. Each 
new layer must add exactly two filled sites and one nearest-neighbor bond (n& = 1/5). This 
constraint implies that no two adjacent sites may be occupied within a layer; if they were, 
then, to conserve n&, two adjacent sites must be occupied in the next, and so on up the tube. 
However, this state is not connected to any other by a single hop. Similarly, the particles 
cannot hop from layer to layer because this adds another intra-layer bond. An allowed 
state can be represented as a string of occupied sites, a = {<7j}, i = 1, . . . , L, which in our 
example is a = {• • • (5, 3) (5, 2) (5, 3) • • •}. At each layer, there are exactly two possibilities 
for the following one. For example, (1,4) can be followed by (1,4) or by (2,4). However, 
the total number of possibilities at any given level is five. Fig. [0] (right) summarizes this 
structure succinctly as a square lattice wrapped on the cylinder. A typical state, then, is 
a lattice path along the tube. There is a recent Hubbard model considered by Henley and 



Zhang [83] of spinless fermions on a square lattice in which the bookkeeping of states is 
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(4,2) (2,5) (5,3) (3,1) (1,4) 

\/\/\/\/\ 
(1,3) (3,5) (5,2) (2,4) (4,1) 

/\/\/\/\/ 
(4,2) (2,5) (5,3) (3,1) (1,4) 

\/\/\/\/\ 
(1,3) (3,5) (5,2) (2,4) (4,1) 

/\/\/\/\/ 
(4,2) (2,5) (5,3) (3,1) (1,4) 

\/ \/\/\/ \ 
(1,3) (3,5) (5,2) (2,4) (4,1) 

/ \/\/\ /\/ 
(4,2) (2,5) (5,3) (3,1) (1,4) 



Figure 6.6: LEFT: Typical configuration in n+ (or mJ) of the (5,0) tube. Alternating numbering 
within layers allows a symmetric description from bottom-to-top or top-to-bottom. RIGHT: Allowed 
states as paths on a wrapped square lattice. The vertex labels may be dropped. 

very similar. 

Generalizing to (N, 0), we find N possible states in each layer and two in the succeeding 
one, and the structure of states is again that of a wrapped square lattice with N squares 
along the circumference. The dimension of the Hilbert space is the number of lattice paths, 
N2 L , so that in an infinitely long tube, the entropy per site is exactly (In2)/A r , as claimed 
earlier. Notice that constrained paths introduce correlations along the length of the tube, 
despite the absence of inter-layer hopping. 

The matrix elements of the projected Hamiltonian connect only those states that differ 
by a single hop: 

[ -It if £,-<W =L-1 
(r\7i\a) = I f"* 11 (6.5) 

I otherwise 

It turns out that this Hamiltonian is exactly solvable, being closely related to a class of 



solid-on-solid models that were introduced by Pasquier [S9|. In the following section we 
derive the continuum limit of Tl, and we confirm the result numerically in the succeeding 
section. 



6. 3. a Continuum Limit 

In the previous section, the paths a were labeled, for clarity, by a string of occupied sites on 
the nanotube. A simpler representation is to work with the wrapped square lattice directly, 
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where the path is uniquely specified by an initial point and its direction in each layer. We 
will label the topmost layer by % = 1 with i increasing by 1 with each downward move. 
There are L layers of hexagons and we impose periodic boundary conditions, L + 1 = 1. In 
order for the layers to match, L has to be even. 

To specify the initial point on a, we chose an "anchor" a on one of the iV sites in the 
i = 1 layer (a = 1, . . . ,N). Now, represent a step to the right in layer i by a fermion, cj, 
and a step to the left by a hole, Cj. A state |cr) in the Hilbert space, S, is represented by 

k> = |a)®4 1 4 , -- c ll > > ( 6 ' 6 ) 

where a is the anchor site in the first layer and i\ ■ ■ ■ i p are the layers where the path steps 
to the right. For instance, the portion of the path in fig. 6J3 is \a) = |3) <8> c\c\c\ • • • |0). 



The fermionic representation is convenient since there is exactly one step in each layer, but 
hard-core bosons can also be used. In any case, in one dimension they are equivalent. The 
number of particles (steps to the right) and holes (steps to the left) must add up to L in 
order for the path to close on itself along the length of the tube. Each path also has a 
topological character for the number of times that it winds around the tube, which must 
be a multiple of N for the path to close. These two conditions may be written as 

N p + N h = L (6.7) 
N p -N h = bN , (6.8) 

where N p h is the number of particles or holes, and b is an integer. If the particles are 
assigned a charge, then bN is the total charge. Note that 6 = corresponds to half-filling, 
N p = N h = L/2. 

Whenever it is allowed within a layer, a single hop changes the step sequence right-left 
to left-right and vice versa, which corresponds to cj +1 Cj or c|cj + i. In a layer without a 
kink no hops are possible, and the hopping terms vanish by fermionic statistics. Since we 
are working in periodic boundary conditions, the boundary terms, c\cl and c^ci, must be 
treated more carefully. A hop at this point is necessarily accompanied by a translation of 
the anchor point by \a) i— ► \a rb 1). Let us represent this operation by 

R±\a) = \a±l) (6.9) 
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with R_ = R + . Cylindrical wrapping requires a Zjv symmetry because \a db N) = \at), i.e. 
R± = R±. Putting the bulk and boundary hopping terms together, the Hamiltonian of 
eqn. fl6,5|) becomes 

rL-i 



n = -2t 



4+i c. + a_ ® 



+ /i.e. . (6.10) 



We can think of Ti as describing free fermions on a periodic one dimensional chain with a 
Ztv impurity on one of the bonds. 

Ti can be diagonalized exactly in momentum space. Going around the tube lengthwise 
contributes a phase e tkL while going around the perimeter contributes e 1 ^, with = 2-rra/N 
(a = 1, . . . , N — 1). Therefore toroidal boundary conditions require 

e ikL e i<t> = 1 (g n ) 

Or, 

where u = 2t is the velocity and n is an integer. The Hamiltonian contains the usual free 
particle dispersion, but with the allowed k given by eqn. ( 6.12] ), 



H = -4t Y cos k c l°k ■ (6-13) 

k 



If the spectrum is linearized around the Fermi momentum, (at half-filling), then at 
small k, nonzero a states cost an additional energy of {2irv/ L)(a/N) 2 . 

The a/N offset in k can be thought of as a minimally coupled vector potential such that 
the magnetic field is a <5-flux tube through the torus containing (a/N) flux quanta. In other 
words, one of the bonds along the chain (the "anchor") had Z^v symmetry, whose effect is 



equivalent to a flux tube. Fig. 6.7 illustrates this equivalence. The offset in k is like a total 
current in the fermion system. 

At this point, one can see two topological effects of the torus. First is the Z^r flux tube, or 
total current. As we have seen, its contribution to the energy near \kp\ was (2irv / 'L)(a / 'N) 2 . 
Second is the path winding along the length of the tube, or total charge, eqn. (|6.7|). Its 
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Figure 6.7: The left ring shows the Zjy impurity on the anchor bond (wavy line). The right ring 
shows the equivalent alternative, where the impurity is replaced by a flux tube through the torus. 

contribution to the energy near \k F \ is similar, (2irv/L)(bN/2) 2 . The total energy due to 
these topological sectors is 

2ttv ( a 2 b 2 N 2 \ 



AE, 



a,b 



+ 



(6.14) 



L \N 2 ' 4 

This expression is familiar from the Luttinger liquid model of one-dimensional spinless 



Fermions [88]. 

We can now obtain the continuum limit of our model. It is well known that free fermions 
in one dimension are equivalent to free bosons. The corresponding Lagrangian is 

1 



C = -[v-\d^) 2 -v(d^) 2 ] 



(6.15) 



where (p is the bosonic field. £ is a conformally invariant theory with central charge c = 1. 
We conjecture that the topological effects that we described above come from compactifying 
ip on a circle of radius R, 



+ 2itR 



(6.16) 



By compactifying the boson, topological modes (or zero modes) appear. In field theory, 
they are conventionally obtained from electric and magnetic monopoles. The energy of the 
zero modes is 



L \R 2 



(6.17) 



where a and b are integers labeling the fundamental cycles on the torus. Comparing E® b 
to A.E a> i, ( |6.14j ), we find that R = N. The ordinary phonon, (oscillator) modes exist on 
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Figure 6.8: The gap scales as 1/N 2 L. 

top of each topological sector and simply contribute the usual phonon energy, so that the 
complete dispersion is 



E = El b+ 2 ^-\n 



(6.18) 



The overall picture of a compactified boson with central charge c = 1 is consistent with 
the solid-on-solid models of Pasquier [|S9|]. 

6.3.b Numerics 



We have diagonalized the original Hamiltonian, eqn. (6.5) numerically with periodic bound- 
ary conditions for system sizes up to N = 11 and L = 10. Due to the sparseness of Tim we 
were also able to obtain the ground state energy up to L = 16. We will fix It = 1 in what 
follows. 

We find that the degeneracy is lifted and the ground state becomes unique and uniform. 
The ground state energy, Eq(L), follows Eq ~ — 0.61L — 0.31ttc/L. The lowest N— 1 excited 
states are given by A a = a 2 A/(N 2 L), with A = 12.9 ± 0.5, which is shown in Fig. 6.8 
for a = 1,2,3. All of these levels are doubly degenerate. This ground state energy and 
spectrum are in perfect agreement with free bosons compactified on a radius R = ( I N, as we 
described in the previous section. 

Right- and left-moving oscillator modes of energy oj n = vk n , where k n = 2irn/L, appear 
in the spectrum, but for a < N the zero modes are the lowest. Our spectrum in Fig. 6.8 
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corresponds to E® b with 6 = 0. The modes with non-zero b are very high in energy and 
are washed out by our small system size. To fix £, we look at higher low-lying levels (which 
also scale like 1/L). We find that the iV'th excitation energy is independent of TV" and 
quadruply degenerate. This can happen only if the iV'th zero mode, E± N0 = 2irv/( 2 L, is 
degenerate with the lowest oscillator mode, uj±\ = 2irv/L, which fixes £ = 1. Thus, the 



compactification radius is R = N. The velocity can be read off from the slopes in Fig. |6^8 
as v = A/2-/T. Within our accuracy, v = 2. The rest of our spectrum is consistent with 
these parameters. For instance, we find a unique, zero-momentum state with a = b = 0, 
which consists of one right- and one left-moving oscillator mode with n = 2n/L at energy 
E = 2E\. Note that £ is fixed only by counting degneracies, not by fitting any parameters. 

6.3.c Higher Order Corrections 

The preceeding discussion is valid to first order in t/V. The next terms are of order t 2 /V 
and involve virtual transitions to adatom configurations that are not in the degenerate 
subspace S. The generic form is 



v 



{ij)(ki) 



V s , (6.19) 



where Vs is a projection operator into S. Another way of writing the second order pertur- 
bation is the familiar form, 

- E\ - -Go- 
where |o"), \a') G S while |A) S. E\ is the energy of the virtual state. E a = E a / are, of 
course, constant, and all energy differences are due to the nearest neighbor repulsion Vriirij. 

There are three types of virtual processes: (i) single particle hopping from a to a ^ a', 
(ii) two particle correlated hopping from a to a ^ a 1 and (iii) single particle diagonal hopping 
from a back into a. For concreteness, consider process (iii) in the (5, 0) state in fig. |6.6| . 
Whenever there is a kink in a, such as in the third layer from the top, the contribution to eqn. 
( 6.20j ) from all virtual hops is — (35/6)4t 2 /y. For example, the adatom on site 3 can hop into 



any one of its six neighbors with the energy denominators 1/2 + 1/2 + 1/2 + 1/2 + 1/2 + 1/3 
(in units of t 2 /V). Similarly, the adatom on site 5 contributes 1 + 1 + 1, for a total of 
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35/6 (it is forbidden to hop one site over to the right because the resulting state is in S). 
On the other hand, if there is no kink, the contribution is —8 • At 2 /V. The criterion for 
a kink in layer i is 2 [1/4 — {hi — l/2)(rtj + i — 1/2)] = 1, where hi = cjcj; otherwise this 
quantity vanishes. Similarly, the absence of a kink is synonymous with the nonvanishing of 
2 [1/4 + {hi — l/2)(nj_|_i — 1/2)]. Thus, the total diagonal contribution to TC can be written 



rC- 
St 2 

~v 



EE 



35 
~6~ 



II; 



At 2 . 



13 ^ _ _ 
— 2^ rnn i+ i 



rii+i- 



+1 



11; 



^ + i-2) + 4 



\*)(*\ 



11; 



(6.21) 



3 ^ ' 3 

1 

For general N, the correction scales like N. The essential term in the last line of eqn. ( 6.21 ) 
is the first one. This four-fermion interaction renormalizes the radius R by corrections of 
order t/V. 

Let us return to processes (i) and (ii). An example of (i) is the adatom in the third 
layer from the top, site 5, hopping to the second layer, site 1, and then back into the third 
layer, site 1. The intermediate state is not in S. This process serves only to renormalize 
t because its amplitude is the same for all kinks. An example of (ii) is the particle in 
the fourth layer, site 5, hopping to site 4, followed by the particle in the third layer, site 
5, hopping to site 1 in the same layer. This correlated hopping occurs in a configuration 
containing the sequence particle-hole-hole or hole-particle-particle, which corresponds to a 
next-nearest neighbor interaction c| +2 c J+i c i+i c i + / l -c We have not analyzed all such terms 
in detail, and it is possible that there is a delicate cancellation of the terms (ii) and (iii) 
when the fields are linearized around \kp\, but we consider it more likely that they do not 
cancel so that R is renormalized at order t/V. 

One can also consider the extreme limit in which t > V. In this case, the XXZ Hamil- 
tonian in eqn. (|6.2| ) is simply the XY model on a cylinder. Let us denote the spin angle 
relative to the cylindrical surface by ip{x, 9), where x is the coordinate along the tube and 
9 is the coordinate around the perimeter. Uniqueness of the wavefunction requires that ip 
has the periodicity (p{x, 9 + 2ir) = tp{x, 9) + 2irm, where m is an integer. The low energy 
excitations are purely along the length of the tube; excitations around the perimeter will 
cost an energy on the order of 1/N, which is large compared to 1/L. Thus we can freeze 
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(a) (b) (c) 




Figure 6.9: The fillings for N = 2. n + = 1/2 (a) and n_ = 1/4 (b). Each adatom can live at 
either site in its layer because each site is connected to every site in the neighboring layers, (c) 
shows the triangular lattice in the plane; the horizontal double bond is due to the periodicity around 
a cylinder, (c) is exactly the geometry of the spin ladder studied by other authors (albeit with 
different coupling). 

the 9 coordinate, and the energy density is proportional to |0 X (/9| 2 . Since the periodicity 
is still tp = (p + 2irm, we end up with a free boson compactified on radius Rxy = 1- The 
question is how the adsorption regime t <C V, which is also a compactified boson but on 
radius R = N, is reached. 



6.4 Special Case: N = 2 

Before concluding with the effective theory, we should point out that the geometry of the 
(2, 0) tube is special; all sites in adjacent layers are interconnected. As a result, all of its 
plateaus have an extensive entropy, and we find that hopping opens a gap in both plateaus. 
Fig. |6.9| illustrates this exception. At either filling, the adatom in each layer is free to hop 
to either site — both configurations are iso-energetic because each site is contiguous to all 
sites in the neighboring layers. Hence both plateaus are macroscopically degenerate. In 
the presence of hopping, each adatom lives in a double well potential, which has a gap of 
order t. In fact, this tube can be written as a spin chain that has been studied at isotropic 



coupling[87], — 2t = V . Two plateaus were found in this case, and it is tempting to speculate 
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whether the two regimes are connected adiabatically. 



6.5 Discussion and Conclusion 

One observable consequence of conformal symmetry is that the low temperature heat ca- 
pacity is fixed by cplfl: 

„ -Kvk 2 R m -Kvk 2 R m 

C = c—^-T = —^-T (6.22) 
o o 

It is noteworthy that, even though the dispersion of the oscillator modes is independent 
of N, the spectrum remembers, via the zero-modes, the finite radius of the nanotube. 
Furthermore, R is quantized by N; in the language of Luttinger liquids, this means that the 
Luttinger parameter is fixed by topology, similarly to the case of edge states in a fractional 
quantum Hall fluid |85|], and in contrast to quantum wires (where the Luttinger parameter 
can vary continuously). Because there is no inter-layer hopping, eft is tied to transverse, 
rather than to longitudinal, density fluctuations along the tube. 

Finally, let us briefly view the spin tube as a quantum spin ladder to see if it yields a gap- 
less state in the degenerate plateaus. A standard approach is to use a Lieb-Schultz-Mattis 
(LSM) argument, in which the spins are deformed slowly along the length pB|. Applying 
it to our tube, we find that a plateau is gapless if S — M is not an integer, where S and 
M are the total spin and magnetization, respectively, per layer (a layer being the N sites 
around the perimeter). Using S = N/2 and the magnetizations from Eqn. ( |6.3j ), we find 
that S — M is an integer in the macroscopically degenerate plateaus, so that the LSM argu- 
ment is insufficient in this case. A conclusive argument must take the geometric frustration 
into account, which is further evidence that our state is strongly correlated. 

In conclusion, we have studied the problem of monolayer adsorption on carbon nanotubes 
and identified several interesting filling fraction plateaus. Since the difference between the 
plateaus decreases slowly, as the inverse of the tube diameter, experimental measurement 
should be feasible for large enough tubes. We have identified the zig-zag tubes as ex- 
ceptional, in which the geometric frustration together with quantum fluctuations lead to 
conformal symmetry. The effective theory is free compactified boson, which has a quantized 
radius to first order in the hopping. The only other such theory in nature that we are aware 
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of are the chiral edge states in a quantum Hall fluid, where the radius is quantized by the 
bulk filling fraction. There are interesting questions related to the large hopping limit. 
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Chapter 7 



Summary 



It is widely believed by both theorists and experimentalists that the composite particle 
construct is required for a full understanding of the fractional quantum Hall effect (FQHE). 
In this thesis we take the point of view that a great deal of the FQHE physics can be 
understood by projecting to the lowest Landau level (LLL) at the outset, which allows us 
to develop a composite particle formalism at various filling fractions, v. Initially the com- 
posites are abstract operators building up the many-particle Fock space, but later analysis 
of physical operators, many-particle wavefunctions, and response functions reveals an in- 
terpretation in terms of vortices bound to an underlying particle. Our theory departs from 
most other theoretical work in this area, which uses singular flux attachment to map the 
problem into a Chern-Simons action. We treat both fermionic and bosonic statistics of the 
underlying particles. 

When the underlying particles (or particles, for short) are bosons, the appropriate formu- 
lation is in terms of composite fermions. The original Hamiltonian for the bosons is mapped 
exactly into a Hamiltonian for composite fermions with an infinite number of constraints. 
To preserve the constraints, we use a conserving approximation to derive an effective theory 
microscopically. In addition to a self-consistent theory of incompressible quantum liquids 
in the LLL, our approach provides a method to calculate the effective mass (gap) and the 
single particle spectrum of the composite fermions, which arise solely from the interactions 
between the particles. This calculation extends previous work in the special case v = 1 to 
arbitrary v. A recent proposal raises the intriguing prospect of observing a FQHE of bosons 



120 



in rotating atomic Bose-Einstein condensates. 

The complementary case is composite bosons. The ground state for a perturbation 
expansion is macroscopically degenerate, which precludes a microscopic derivation of an 
effective theory. We follow an alternate approach by constructing a phenomenological 
Landau-Ginzburg action based on a symmetry analysis. A crucial ingredient in the ac- 
tion is the internal structure of the composite particle, which, at v = 1/p, consists of p 
vortices bound to the particle. For p = even the particles are bosons while for p = odd 
they are fermions. Our model incorporates this structure through gauge potentials that 
couple to the internal degrees of freedom. The spectrum of the effective theory contains the 
so-called magnetoroton excitation, which seems to be the first analytic observation of this 
phenomenon. 

The next portion of the thesis is an examination of paired states of bosons or fermions in 
two dimensions. The class of pairings that we investigate includes non-zero relative angular 
momentum, which is a state that breaks both parity and time reversal invariance. Part of 
our method is based on Bardeen-Cooper-Schreiffer (BCS) theory and is largely independent 
of the FQHE. However, in the context of the FQHE, we use a mean field approximation 
that maps particles in a net magnetic field into composite particles in zero field. The BCS 
approach then applies to the composite particles. 

For bosonic particles, we consider p-wave pairing of spin-1/2 bosons and d-wave pairing 
of spinless bosons. In the FQHE, the former case is a singlet known as the "permanent" 
state and the latter is the "Haffnian" . Each state can be written as a trial wavefunction that 
is the unique ground state of its corresponding Hamiltonian. By analyzing the spectrum 
of the Hamiltonian directly or by applying the mean field BCS theory, we find that the 
permanent sits on the transition between the polarized Laughlin state (a Bose condensate 
of composite particles) and a Bose condensate with helical order. The spin order of the 
permanent is that of an anti-Skyrmion. Similarly, the Haffnian is on the transition between 
a Laughlin state and a strong coupling paired state. Of course, these conclusions can stand 
alone without reference to the FQHE because they can be derived solely within the BCS 
framework. 

For fermionic particles, we carry through a conserving approximation for the spin con- 
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ductivity. We show that the induced action for an external gauge field that couples to spin 
is a Chern-Simons (CS) term. In the d-wave case we obtain the non-abelian SU(2) CS term 
because the system is spin-rotationally invariant, whereas the p-wave case gives an abelian 
U(l) term since it has only U(l) symmetry. In both cases, the Hall spin conductivity is a 
topological invariant, which characterizes the winding of the order parameter in momentum 
space. This is a microscopic proof of quantization that was proposed in earlier works by 
other authors. 

The last part of the thesis deals with adsorption on carbon nanotubes, which is a one- 
dimensional problem. The hexagon centers serve as adsorption sites for hard core atoms, 
allowing the system to be treated as a lattice gas on a triangular lattice wrapped on a 
cylinder. This model is equivalent to a type of quantum spin tube. The wrapping introduces 
geometric frustration on top of the frustration of the triangular lattice, leading to interesting 
physics. In the spin language, we find magnetization plateaus in all tubes in the Ising 
limit, which is confirmed both analytically and numerically. However, the zig-zag tubes 
are exceptional and contain plateaus that are macroscopically degenerate. When quantum 
hopping is allowed, the special plateaus become gapless phases that are described by a c = 1 
conformal theory of a compactified boson. The theory is derived analytically and confirmed 
numerically. Perhaps the most remarkable feature is that the radius of compactification is 
quantized by the tube diameter. This brings us back to the FQHE, where the theory of 
edge states is also characterized by a (chiral) conformal boson compactified on a quantized 
radius. 
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Appendix: Non- Commutative 
Fourier Transform 



Our ultimate aim in this subsection is to introduce constructs that will allow us to take the 
thermodynamic limit, N —* oo, and to introduce the momentum k so as to take advantage 



of translational invariance. As discussed by several authors |L2, 15, [L8|], k is a good quantum 
number and z x k will turn out to be the dipole moment. The formalism in this appendix 
has been discussed in more detail by |l8|], and we include it here for completeness. 

Consider the special case of one attached vortex with B\ = —B2, as introduced in 
Sections 2.1.b and |2.1.c| . We first introduce real space wavefunctions by analogy to the 



usual matter field: 

c(z,rj) = ^u m (z)u n {rf) c mn 

ran 

c j (v,z) = ^2u n (r])u m (z) c ] nm . (7.1) 

ran 

Our convention is to use z for the left coordinate and r] for the right. Complex conjugation 
reflects the two opposite charges. There is only one magnetic length £b m this problem 
because the magnitude of the charges is equal. This is a rather singular limit of our two- 
particle construction in Section 2.1.b| ; in this limit the effective magnetic field is B = 



B\ + B2 = 0, and the pseudomomentum and translation operators, tt and K, are identical 
(eqn. Q). 

In the z, 77 basis, the densities become 

P R (V,V') = J d 2 z c j (r],z)c(z,rj') 

p L (z,z') = J d 2 V c\nJ)c{z,Ti) (7.2) 
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Thus integration has replaced summation over indices. It is convenient to introduce a binary 
operation * to represent integration over one set of coordinates, 

(a * b)(z,z') = / d 2 zi a(z,zi)b(zi,z') , (7-3) 



of two operators, a and b, which is just matrix multiplication. It is also convenient to define 
the *-commutator by 



[a* b] = a*b- b* a . (7.4) 



This allows us to write 



p R = ct * c 

p L = : c * c f : (7.5) 

Note the normal ordering in p L necessary to avoid sign ambiguities. 

Consider the plane wave, e* kr , projected to the LLL. Following the previous discussion 



of two particles in Section 2.1.b, eqn. fl2.22|), r = (Ri + R-2)/2 in zero effective field. The 



plane wave operator now becomes 

e ik-r = e ik-(R 1 +R 2 )/2 _ ^ 

Its representation in the z, r] coordinates is obtained by acting on the lowest weight eigen- 
function Vo,o- Recall that, for B\ = —E>2, ijjofi(z,v) = ^e^^ z ^~^ 2+ ^ zT ' , where the 
magnetic length has been set to unity. It is not difficult to show that ipQfi(z,r]) is identical 
to the delta function in the LLL. 

5(z,rj) =^2u m (z)u m (r]) . (7.7) 

m 

In the LLL, 5 acts as expected: S*a = a*5 = a. The constraint is thus p R = 5, or 

p R (rj,rf) = 6(rj,rf) . 

The differential representation of Ri 2 in the z, ij coordinates was constructed in Sec- 
tion ( |2.1.bD . According to the prescription, the plane wave acting on 5 yields the LLL 
representation 

Tic{z,rj) = 5{z,rj)e^ Iz+k ^-i\ k \ 2 . (7.8) 
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By either straightforward integration or by using the commutator of Ri^, eqn. ( |2.19| ), we 
find that the f k obey 



r k *r k , = f k+k , e2 lkAk , (7.9) 

where we have introduced the shorthand notation, Ak = — z x k and k A k' = k • Ak'. 
The phase factor is the area of a triangle formed by k and k' so that the phase counts the 
flux enclosed by the triangle. Hence we interpret fk as a magnetic translation in the plane 



by Ak IS]. The connection to plane waves extends to completeness and orthonormality 
properties, which defines a "noncommutative Fourier transform". In particular, 

Tr f k * f k ' = 27r<5(k + k') , (7.10) 
— T k (z,z / ) T -k(??,?? / ) = 5( z >v')5{v,z') ■ (7-11) 

The Tr stands for a trace defined by Tr a = J d 2 z a(z,z). We can now define the Fourier 
transform and its inverse: 

c(z,rj) = j , 2 s 3/2 Ck T k (z,r?) , (7.12) 
c k = (2vr) 1/2 Tr c*f_ k • (7.13) 

The extra \p2/n factors are not used in general; they are specific to the Fourier transform of 
the c's in order for the commutators to retain their conventional form 

[c k , 4,] ± = (27r) 2 <5(k - k') . (7.14) 

Therefore, we have constructed composite particles with momentum k and dipole moment 
Ak, which emerged from magnetic translations perpendicular to k. When the underlying 



particles are bosons, anticommuators are required in eqn. ( 7.14 ). This defines the composite 
fermion. The Fock space of fermions at v = 1 is trivial, consisting of exactly one function, 
so composite bosons are not useful in this case. 

Next we would like to apply the Fourier transform to the left and right densities in eqn. 
( |7.2| ), which become 
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The commutators of /5 q are familiar in the quantum Hall effect [f)5| , p0| , K|, defining an 
infinite Lie algebra known as Woo ■ The constraints appear particularly simple in momentum 
space: 

- 2vrp $(q) = , (7.17) 
which enforces a uniform vortex density, 7J. 
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